mian.m

来自「罚函数法matlab程序.希望对大家有用.」· M 代码 · 共 37 行

%主函数main.m
clc
clear
a=input('请输入根的上限：');
b=input('请输入根的下限：');
e=input('请输入控制精度：');

count=0;
if f1(a)*f1(b)<0
    while abs(b-a)>e
          count=count+1;
          x=1/2*(a+b);
          if f1(x)*f1(a)>0
          a=x; 
          elseif f1(x)*f1(a)<0
          b=x;
          elseif f1(x)==0
             return 
          end
          
   end
elseif f1(a)*f1(b)>0
   fprintf('无法用试位法！');
   return;
elseif f1(a)==0
    x=a;
elseif f1(b)==0
    x=b;
end
fprintf('2分法的解：\n');
x
fprintf('2分法迭代次数：\n');
count