clc
clear
syms x
f=input('请输入函数f(x):');%x^3/3-x
x0=input('请输入迭代初值:');%初值大小
fai=x-f/diff(f,x); %求牛顿迭代函数
k=0 ; %迭代次数
err=1.000;
phi=input('请输入容许误差:');%允许误差大小
disp('迭代次数 求得根 误差');
while(err>phi)
x1=subs(fai,x0);
err=abs(x1-x0);
x0=x1;
k=k+1;
end
fprintf('%d %1.10f %1.10f\n',k,x1,err);
matlab牛顿迭代法代码
最新推荐文章于 2024-05-22 11:32:00 发布