用matlab解非線性方程組

    使用matlab的fsolve函數，首先需要寫一個m檔案來描述非線性方程組，比如

function eq=nxxf(x)

global number;

number=number+1

eq(1) = x(1)*normcdf(x(3),0,1)-80120231*normcdf(x(4),0,1)*exp(-0.02)-595042771.7;

eq(2) = (log(x(1)/80120231)+0.02+0.5*(x(2)^2))/x(2)-x(3);

eq(3) = x(3)-x(2)-x(4);

eq(4) = (normcdf(x(3),0,1)*x(1)*x(2))/595042771.7-0.0475;

其中每個方程都是移項後一邊為0後，另外一邊的式子。然後在matlab命令視窗輸入global number;number=0;y=fsolve('nxxf',[1,1,1,1])number其中fsolve函數第一個參數就是非線性方程組的函數，第二個是初始點，將輸出結果和迭代步數