syms?x1?x2?M;??????????????%M為罰因子。?
m（1）=1;?????????????????????
c=5;???????????????????????%c為遞增系數(shù)。賦初值。?
a（1）=-3;?
b（1）=-3;???????????????????????????????
f=x1+x2+M*（x1^2+x2^2-2）^2;%f=x1^3+x2^2+M*（x1+x2-1）^2;%外點罰函數(shù)?
f0（1）=1;?
?
%求偏導(dǎo)、Hessian元素?
fx1=diff（f‘x1‘）;?%求x1偏導(dǎo)?
fx2=diff（f‘x2‘）;%求x2偏導(dǎo)?
fx1x1=diff（fx1‘x1‘）;%x1二階導(dǎo)?
fx1x2=diff（fx1‘x2‘）;?
fx2x1=diff（fx2‘x1‘）;?
fx2x2=diff（fx2‘x2‘）;?
?
%外點法M迭代循環(huán)?
for?k=1:100???????????????????????????????????????????????????????????
????x1=a（k）;x2=b（k）;M=m（k）;?
????%牛頓法求最優(yōu)值?
????for?n=1:100?????????????????????????????????????????????????????
????????f1=subs（fx1）;?%求解梯度值和Hessian矩陣?
????????f2=subs（fx2）;?
????????f11=subs（fx1x1）;?
????????f12=subs（fx1x2）;?
????????f21=subs（fx2x1）;?
???