科学计算案例_牛顿法解三次方程
08-25
%牛顿法应用:求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