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;