%罰函數法??
%目標函數：fun=x1^2+2*x2^2;?
%約束條件：x1+x2=1?

syms?x1?x2?k?
fun0=x1^2+2*x2^2?
fun1=x1^2+2*x2^2+k*（x1+x2-1）^2?
f0=[diff（fun1x1）;diff（fun1x2）];?
[x1x2]=solve