牛顿法如何解非线性规划

[灰烬重生]

目标函数：
function [f,df,d2f]=nwfun(x);
f=x(1)^4+25*x(2)^4+x(3)^4+x(1)^2*x(2)^2+4*x(1)^2*x(3)^2+16*x(2)^2*x(3)^2;
约束条件：
function [g,ceq] = nyueshu(x)
g(1) = x(1)*x(1)+x(2)*x(2)+x(3)*x(3)-100
g(2) = 60 - x(1)*x(1) + 10*x(3)*x(3)
ceq(1) = x(1) + x(2)*x(2) + x(3) - 80
ceq(2) = x(1)^3 + x(2)*x(2) + x(3) - 80
那主程序如何编译，要用牛顿迭代法来解，不是用fmincon来解，哪位大侠知道，求帮助

[laosam280]

这个你的约束都是等式约束吗？还是包含了不等式和等式？
有约束问题最好用拉格朗日乘子法和惩罚函数法求解
牛顿法用于求解无约束问题

[灰烬重生]

是含有约束的，我们老师是说用拉格朗日牛顿法来解，应该就是你说的拉格朗日乘子法，能给我具体说说怎么运用的么，有代码最好了，谢谢你了