Displaying 1 result from an estimated 1 matches for "loesung".
2007 Sep 03
2
The quadprog package
...erden
dvec
mu<-0 # (in Mio. €)
bvec<-c(1,mu,matrix(0,7,1)) # muss als Spaltenvektor eingegeben werden
bvec
mu_r<-c(19.7,33.0,0.0,49.7, 82.5, 39.0,11.8)
Amat<-matrix(c(matrix(1,1,7),7*mu_r,diag(1,7,7)),9,7,byrow=T)
# muss als Matrix angegeben werden, wie sie wirklich ist
Amat
meq<-2
loesung<-solve.QP(Dmat,dvec,Amat=t(Amat),bvec=bvec,meq=2)
loesung
# Überprüfen, ob System richtig gelöst wurde
loesung$solution %*% mu_r
sum(loesung$solution)
for (i in 1:7){
a<-loesung$solution[i]>=0
print(a)
}
Thanks in advance for your answers.
______________________________
Thomas Schwan...