function y=M(x)
global lumda sigma l m;
c=confun(x);
sum=0;
for i=1:l
    sum=sum+c(i)^2;
end
sum=sigma/2*sum;
for i=1:l
    sum=sum-lumda(i)*c(i);
end
for i=(l+1):m
    if lumda(i)-sigma*c(i)>=0;
        temp=(lumda(i)-sigma*c(i))^2-lumda(i)^2;
    else
        temp=-lumda(i)^2;
    end
    sum=sum+1/(2*sigma)*temp;
end
y=fun(x)+sum;