c4_1.m

这是我从一本基于MATLAB实现的数值分析算法的书中自带的光盘程序

%C4_1.m
%二分法求解非线性方程的根
function C4_1
a=input('a=');
b=input('b=');
k=0;
epsilon=1.0e-4;
it_limit=20
while (b-a)/2>epsilon
      k=k+1
      if k>=it_limit  break;  end
      y_a=fun_f(a);
      x0=(b+a)/2
      y_x0=fun_f(x0)
      y_b=fun_f(b);
      if y_x0==0 
         fprintf('方程的根x = ');
      else if y_a*y_x0<0
             b=x0;
           else  a=x0;
          end
      end
end
fprintf('方程的根x = ');
x=(a+b)/2