Function obj(x);
Begin
Return 2*x[0]^2-4*x[0]*x[1]+4*x[1]^2-6*x[0]-3*x[1];
End;
//非线性约束:
function constr(x);
begin
ne := array();
e := array() ;
//ne[0]:= x[0]+x[1]^2-5;//x[0]+x[1]^2<5非线性不等式约束
e[0]:=x[1]^2+x[0]^2-4;//x[1]^2=4非线性等式约束
return array(ne,e);
end