惩罚函数最优化方法 MATLAB实现

%罚函数法  
%目标函数：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