如 此 循 环 往 复 , 可 得 一 列 逼 近 方 程 精 确 解 的 点
,其一般表达式为:
该公式所表述的求解方法称为 迭代法或切线法。
程序:
function y=f(x)%定义原函数
y=x^3-x-1;
end
function y1=f1(x0)%求导函数在 x0 点的值
syms x;
t=di(f(x),x);
y1=subs(t,x,x0);
end
function newton_iteration(x0,tol)%输入初始迭代点 x0 及精度 tol
x1=x0-f(x0)/f1(x0);k=1;%调用 f 函数和 f1 函数
while abs(x1-x0)>=tol
x0=x1;x1=x0-f(x0)/f1(x0);k=k+1;
end
fprintf('满足精度要求的数值为 x(%d)=%1.16g\n',k,x1);
fprintf('迭代次数为 k=%d\n',k);
end
3