%編寫目標函數
function y=objfun(x)
y=exp(x(1))*(3*x(1)^2+2*x(2)^2+3*x(1)*x(2)+2*x(2)+1);

%編寫返回約束值的函數
function [c,ceq]=confun(x)
%非線性不等式約束
c=[1.5+x(1)*x(2)-x(1)-x(2);-x(1)*x(2)-10];
%線性等式約束
ceq=[];

%初始點
x0=[-1,1]; 
%設置下界
lb=[0,0]; 
%無上界
ub=[]; 
%采用標準算法
options=optimset ('largescale','off');
[x,fval]=fmincon('objfun',x0,[],[],[],[],lb,ub,'confun',options)
[c,ceq]=confun(x)