%牛顿法应用:求x^3-x-1=0的根:  迭代公式 x(k+1)=x(k)-f(x(k))/df(x(k))
clc
f=inline('x^3-x-1');
df=inline('3*x^2-1');
c=0.5*10^-6;
x(1)=input('请输入迭代初值=');
x(2)=x(1)-f(x(1))/df(x(1));
fprintf('k=%.0f, x=%.8f\n',1,x(1))
fprintf('k=%.0f, x=%.8f\n',2,x(2))
k=2;
while abs(x(k)-x(k-1))>=c
    k=k+1;
    x(k)=x(k-1)-f(x(k-1))/df(x(k-1));
    fprintf('k=%.0f, x=%.8f\n',k,x(k))
end
% 请输入迭代初值=0.6
% k=1, x=0.60000000
% k=2, x=17.90000000
% k=3, x=11.94680233
% k=4, x=7.98552035
% k=5, x=5.35690931
% k=6, x=3.62499603
% k=7, x=2.50558919
% k=8, x=1.82012942
% k=9, x=1.46104411
% k=10, x=1.33932322
% k=11, x=1.32491287
% k=12, x=1.32471799
在线计算器% k=13, x=1.32471796

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。