newtondiedai.m

来自「基于matlab实现的实质分析实例的计算 对基本的插值的计算 可运行」· M 代码 · 共 30 行

function x=newtondiedai(x0,n,e) %方程求根牛顿迭代(课本137页例5)

for k=1:n
    k
    
    
    if x0==-1
        '奇异'
        break
    else
        x1=x0-(x0-exp(-x0))/(1+x0);
    end
    if abs(x1-x0)<e
        x=x1;
        break
    else
        if k==n
            '迭代失败'
        else
            x0=x1;
        end
    end
    x1
    
end

x=x1;