diff options
Diffstat (limited to 'demos/qpipoptmat.dem.sce')
-rw-r--r-- | demos/qpipoptmat.dem.sce | 36 |
1 files changed, 16 insertions, 20 deletions
diff --git a/demos/qpipoptmat.dem.sce b/demos/qpipoptmat.dem.sce index 61263a8..bbaa42c 100644 --- a/demos/qpipoptmat.dem.sce +++ b/demos/qpipoptmat.dem.sce @@ -3,26 +3,6 @@ mode(1) // Demo of qpipoptmat.sci // -//Find x in R^6 such that: -halt() // Press return to continue - -Aeq= [1,-1,1,0,3,1; --1,0,-3,-4,5,6; -2,5,3,0,1,0]; -beq=[1; 2; 3]; -A= [0,1,0,1,2,-1; --1,0,2,1,1,0]; -b = [-1; 2.5]; -lb=[-1000; -10000; 0; -1000; -1000; -1000]; -ub=[10000; 100; 1.5; 100; 100; 1000]; -x0 = repmat(0,6,1); -param = list("MaxIter", 300, "CpuTime", 100); -//and minimize 0.5*x'*Q*x + p'*x with -f=[1; 2; 3; 4; 5; 6]; H=eye(6,6); -[xopt,fopt,exitflag,output,lambda]=qpipoptmat(H,f,A,b,Aeq,beq,lb,ub,[],param) -clear H f A b Aeq beq lb ub; -halt() // Press return to continue - //Find the value of x that minimize following function // f(x) = 0.5*x1^2 + x2^2 - x1*x2 - 2*x1 - 6*x2 // Subject to: @@ -37,6 +17,22 @@ b = [2; 2; 3]; lb = [0; 0]; ub = [%inf; %inf]; [xopt,fopt,exitflag,output,lambda] = qpipoptmat(H,f,A,b,[],[],lb,ub) +// Press ENTER to continue halt() // Press return to continue +//Find x in R^6 such that: +Aeq= [1,-1,1,0,3,1; +-1,0,-3,-4,5,6; +2,5,3,0,1,0]; +beq=[1; 2; 3]; +A= [0,1,0,1,2,-1; +-1,0,2,1,1,0]; +b = [-1; 2.5]; +lb=[-1000; -10000; 0; -1000; -1000; -1000]; +ub=[10000; 100; 1.5; 100; 100; 1000]; +x0 = repmat(0,6,1); +param = list("MaxIter", 300, "CpuTime", 100); +//and minimize 0.5*x'*Q*x + p'*x with +f=[1; 2; 3; 4; 5; 6]; H=eye(6,6); +[xopt,fopt,exitflag,output,lambda]=qpipoptmat(H,f,A,b,Aeq,beq,lb,ub,[],param) //========= E N D === O F === D E M O =========// |