惩罚函数法（内点法、外点法）求解约束优化问题最优值 matlab

1、 用外点法求下列问题的最优解

方法一：外点牛顿法：
clc
m=zeros(1,50);a=zeros(1,50);b=zeros(1,50);f0=zeros(1,50);%a b为最优点坐标，f0为最优点函数值，f1 f2最优点梯度。 m=zeros