还剩2页未读,继续阅读
文本内容:
实验一求方程根的数值方法牛顿法的算法程序clear;x=2;fori=1:10x=x-x^3-3*x-1/3*x^2-3end结果x=
1.8889x=
1.8795x=
1.8794x=
1.8794x=
1.8794x=
1.8794x=
1.8794x=
1.8794x=
1.8794x=
1.8794故方程的根为x=
1.
879.双点弦算法的程序#includestdio.hvoidmain{floatx1=
1.0x2=
2.0x;inti;fori=1;i=5;i++{x=x1*x2*x2*x2-3*x2-1-x2*x1*x1*x1-3*x1-1/x2*x2*x2-3*x2-1-x1*x1*x1-3*x1-1;printf%f\nx;x1=x2;x2=x;}getch();}结果:故方程的根为x=
1.879Aitken法算法程序symsfxgxx;gx=sinx*2;fx=
0.5*x-sinx;dispkxx1x2x=
1.5;k=0;ffx=subsfxxx;whileabsffx
0.0001;u=subsgxxx;v=subsgxxu;disp[num2strknum2strxnum2strunum2strv]x=v-v-u^2/v-2*u+x;k=k+1;ffx=subsfxxx;enddisp[num2strknum2strxnum2strunum2strv]fprintf所求的解是x=%f迭代步数是%d\nxk结果kxx1x
201.
51.
9951.
822711.
86721.
91281.
884221.
89521.
89571.
895431.
89551.
89571.8954所求的解是x=
1.895494迭代步数是3。