function PHR()
global lumda sigma  l m;
x=[2 1];            %初始點
l=1;                %等式
m=1;                %不等式
lumda=1;            %乘子向量
sigma=10;            %初始罰因子
g=2;                %放大系數
eps=1e-32;          %控制誤差
sita=0.5;           %常數
temp=terminate(x);  %終止函數
while 1
    x=fminsearch(@M, x);       %對M求最優解
    fai=terminate(x);
    if fai<eps
        break;
    end
    if temp<eps
        break;
    end
    if fai/temp>=sita
        sigma=g*sigma;
    end
    temp=fai;
    c=confun(x);
    for i=1:l                                   %乘子迭代
        lumda(i)=lumda(i)-sigma*c(i);
    end
    for i=(l+1):m
        if (lumda(i)-sigma*c(i))>=0
            lumda(i)=lumda(i)-sigma*c(i);
        else
            lumda(i)=0;
        end
    end
end
opt=x;
disp('The optimum solution is : ');
disp(opt);
disp('The optimum value is: ');
disp(fun(opt));