%%%ex1208.m 拉各朗日乘子法求最優化解
x=zeros(1,2);
%用syms表示出轉化后的無約束函數
syms x y lama
f=x+y+lama*(x^2+y^2-2);
%分別求函數關于x、y、lama的偏導
dx=diff(f,x);
dy=diff(f,y);
dlama=diff(f,lama);
%令偏導為零求解x、y
xx=solve(dx,x); %將x表示為lama函數
yy=solve(dy,y);  %將y表示為lama函數
ff=subs(dlama,{x,y},{xx,yy}); %代入dlama得關于lama的一元函數
lamao=solve(ff); %求解得lamao
xo=subs(xx,lama,lamao) %求得取極值處的xo
yo=subs(yy,lama,lamao)  %取極值處的yo
fo=subs(f,{x,y,lama},{xo,yo,lamao}) %極值點函數值