newtondiedai_exercise4_16.m

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

function x=newtondiedai(x0,n,e) %方程求根牛顿迭代(课本154页16题)

for k=1:n
    k
    
    
    if 3*x0^2+4*x0+10==0
        '奇异'
        break
    else
        x1=x0-(x0^3+2*x0^2+10*x0-20)/(3*x0^2+4*x0+10);
    end
    if abs(x1-x0)<e
        x=x1;
        break
    else
        if k==n
            '迭代失败'
        else
            x0=x1;
        end
    end
    x1
    
end

x=x1;