求大侠指点牛顿法如何解非线性规划

灰烬重生

目标函数：
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来解，哪位大侠知道，求帮助