数值分析 实验四 方程求根程序设计及验证

《数值分析 实验四 方程求根程序设计及验证》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验四方程求根程序设计与验证二分法程序如下：functionx=nabisect(fname,a,b,e)%fnameΪÄÚǶº¯ÊýµÄ±í´ïʽf(x)%a,bÇø¼äµÄ¶Ëµã%eΪ¾«¶Èifnargin<4,e=1e-4;endfa=subs(fname,a);fb=subs(fname,b);iffa*fb>0,error('º¯ÊýÔÚÁ½¶ËµãÖµ±ØÐëÒìºÅ');endx=(a+b)/2while(b-a)>(2*e),fx=subs(fname,x);iffa*fx<0,b=x;fb=fx;elsea=x;fa=fx;endx=(a+b)/2;endsym

2、sx;实验结果：>>f=x^3-3*x-1f=x^3-3*x-1>>x=nabisect(f,1,3)x=2x=1.8793x=nabisect(f,-1,0)x=-0.5000x=-0.3474>>x=nabisect(f,-2,-1)x=-1.5000x=-1.5320牛顿法程序如下：functionx=nanewton(fname,dfname,x0,e,N)ifnargin<5,N=500;endifnargin<4,e=1e-4;endx=x0;x0=x+2*e;k=0;whileabs(x0-x)>e&k

3、(dfname,x0);disp(x)endifk==N,warning('已达函数上线');end实验结果为：x=nanewton(f,df,2.0,0.0001,10)1.88891.87951.8794x=1.8794讨论分析：