newton.m

用newton迭代法求矩阵,可以求任意阶的矩阵,非常好用

function[pl,err,k,y]=newton(f,de,p0,delta,max1)
%f是非线性函数
%df是f的微商
%p0是初始值
%delta是给定允许误差
%max1是迭代的最大次数
%p1是牛顿迭代法求得的方程的近似值
%err是p0的误差估计
%k是迭代次数
%y=f(p1)
p0,feval('f',p0)
for k=1:max1
    p1=p0-feval('f',p0)/feval('df',p0);
    err=abs(p1-p0);
    p0=p1;
    p1,err,k,y=feval('f',p1)
    if(err<delta)|(y==0)
        break;
    end
    p1,err,k,y=feval('f',p1)
end