1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
|
//User defined hessian is of wrong dimensions
function y=f(x)
y=x(1)*x(2)+x(2)*x(3);
endfunction
x0=[1,1,1];
A=[];
b=[];
Aeq=[];
beq=[];
lb=[0 0.2,-%inf];
ub=[0.6 %inf,1];
intcon=[1,2];
// !--error 10000
// intfmincon_options: Expected type ["string"] for input argument GradObj at input #6, but got "function" instead.
// at line 56 of function Checktype called by :
// at line 359 of function intfmincon called by :
// [x,fval,exitflag,grad,hessian] =intfmincon(f,x0,intcon,A,b,Aeq,beq,lb,ub,nlc,options)
// at line 46 of exec file called by :
// exec intfmincon_hess2.sce
function [c,ceq]=nlc(x)
c=[x(1)^2-1,x(1)^2+x(2)^2-1,x(3)^2-1];
ceq=[x(1)^3-0.5,x(2)^2+x(3)^2-0.75];
endfunction
function y= fGrad(x)
y= [x(2),x(1)+x(3),x(2)];
endfunction
//Wrong dimensions of Hessian here
function y= lHess(x,obj,lambda)
y= obj*[0,1;1,0;0,1] + lambda(1)*[2,0;0,0;0,0] + lambda(2)*[2,0;0,2;0,0] +lambda(3)*[0,0;0,0;0,0] + lambda(4)*[6*x(1),0;0,0;0,0] + lambda(5)*[0,0;0,2;0,0];
endfunction
options=list("MaxIter", [1500], "CpuTime", [500], "GradObj", fGrad, "Hessian", lHess);
[x,fval,exitflag,grad,hessian] =intfmincon(f,x0,intcon,A,b,Aeq,beq,lb,ub,nlc,options)
|