第四周
3.
中的三个根。
,在求8] [0,
041.76938.7911.1-)(23=-+=x x x x f
function y=mj()
for x0=0:0.01:8
x1=x0^3-11.1*x0^2+38.79*x0-41.769;
if (abs(x1)<1.0e-8)
x0
end
end
4.分别用简单迭代法、埃特金法、牛顿法求解方程,并比较收敛性与收敛速度(ε分别取10-3、10-5、10-8)。
简单迭代法:
function y=jddd(x0)
x1=(20+10*x0-2*x0^2-x0^3)/20;
k=1;
while (abs(x1-x0)>=1.0e-3)
x0=x1;
x1=(20+10*x0-2*x0^2-x0^3)/20;k=k+1;
end
x1
k
埃特金法:
function y=etj(x0)
x1=(20-2*x0^2-x0^3)/10;
x2=(20-2*x1^2-x1^3)/10;
x3=x2-(x2-x1)^2/(x2-2*x1+x0);
k=1;
while (abs(x3-x0)>=1.0e-3)
x0=x3;
x1=(20-2*x0^2-x0^3)/10;
x2=(20-2*x1^2-x1^3)/10;
x3=x2-(x2-x1)^2/(x2-2*x1+x0);k=k+1;
end
2,020102)(023==-++=x x x x x f