newdim.m

这是牛顿法求解非线性方程组哦那个

function [[x;y],iter,err]=newdim(F,JF,P,delta,epsilon,max1)
x=P(1);y=P(2);
Y=feval('F',P);
for k=1:max1
    J=feval('JF',P);
    Q=P-(J\Y')';
    Z=feval('F',Q);
    err=norm(Q-P);
    relerr=err/(norm(Q)+eps);
    P=Q;
    Y=Z;
    iter=k;
    if (err<delta)|(relerr<delta)|(abs(Y)<epsilon)
        break
    end
end