matlab编的二分法求根的程序

matlab编的二分法求根的程序

matlab编的二分法求根的程序

2010-03-10 16:44

function varargout=erfliu(varargin)

clear,clc

F=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)\n'); while((b-a)>tol)

m=a+(b-a)/2;

fprintf('%10.6f %10.6f %10.6f %10.6f\n',a,feval(F,a),b,feval(F,b)'); if feval(F,a)*feval(F,m)>=0

a=m;

else

b=m;

end

end

结果

a f(a) b f(b) 1.000000 -2.365884 3.000000 8.435520

1.000000 -2.365884 2.000000 0.362810

1.500000 -1.739980 2.000000 0.362810

1.750000 -0.873444 2.000000 0.362810

1.875000 -0.300718 2.000000 0.362810

1.875000 -0.300718 1.937500 0.019849

1.906250 -0.143255 1.937500 0.019849

1.921875 -0.062406 1.937500 0.019849

1.929688 -0.021454 1.937500 0.019849

1.933594 -0.000846 1.937500 0.019849

1.933594 -0.000846 1.935547 0.009491

1.933594 -0.000846 1.934570 0.004320

1.933594 -0.000846 1.934082 0.001736

1.933594 -0.000846 1.933838 0.000445

1.933716 -0.000201 1.933838 0.000445

1.933716 -0.000201 1.933777 0.000122

1.933746 -0.000039 1.933777 0.000122

1.933746 -0.000039 1.933762 0.000041

1.933746 -0.000039 1.933754 0.000001

1.933750 -0.000019 1.933754 0.000001

1.933752 -0.000009 1.933754 0.000001

1.933753 -0.000004 1.933754 0.000001

第1页下一页