matlab编的二分法求根的程序.doc

《matlab编的二分法求根的程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、matlab编的二分法求根的程序2010-03-1016:44functionvarargout=erfliu(varargin)clear,clcF=inline('x^2-4*sin(x)');a=1;b=3;erff(F,a,b)function  erff(F,a,b)  tol=0.0000005;  fprintf('   a         f(a)         b        f(b)');while((b-a)>tol)m=a+(b-a)/2;fprintf('%10.6f%10.6f%10.6

2、f%10.6f',a,feval(F,a),b,feval(F,b)');iffeval(F,a)*feval(F,m)>=0    a=m;else    b=m;endend结果       a               f(a)               b          f(b)1.000000  -2.365884   3.000000   8.4355201.000000  -2.365884   2.000000   0.3628101.500000  -1.739980   2.000000  

3、 0.3628101.750000  -0.873444   2.000000   0.3628101.875000  -0.300718   2.000000   0.3628101.875000  -0.300718   1.937500   0.0198491.906250  -0.143255   1.937500   0.0198491.921875  -0.062406   1.937500   0.0198491.929688  -0.021454   1.937500   0.0198491.933594 

4、 -0.000846   1.937500   0.0198491.933594  -0.000846   1.935547   0.0094911.933594  -0.000846   1.934570   0.0043201.933594  -0.000846   1.934082   0.0017361.933594  -0.000846   1.933838   0.0004451.933716  -0.000201   1.933838   0.0004451.933716  -0.000201   1.933

5、777   0.0001221.933746  -0.000039   1.933777   0.0001221.933746  -0.000039   1.933762   0.0000411.933746  -0.000039   1.933754   0.0000011.933750  -0.000019   1.933754   0.0000011.933752  -0.000009   1.933754   0.0000011.933753  -0.000004   1.933754   0.000001