计算方法中方程的近似解法中二分法matlab实现

计算方法

方程的近似解法

二分法

fun.m中代码如下：

function fun(a,b,e)
%f是自定义的函数
%a为隔根区间左端点，b为隔根区间右端点，e为绝对误差限
if nargin==2
    e=1.0e-6;
elseif nargin<2
    input('变量输入错误！');
    return;
end
if a>=b
    input('隔根区间输入错误！');
    return;
end
a1=a;
b1=b;
c1=(a1+b1)/2;
n=0;  %迭代计数器,初值为0
while  (b-a)/(2^(n)) >= 1/2*e
    c1
    if f(c1)==0
        c1
    elseif f(a1)*f(c1)>0
        a1=c1;
        c1=(a1+b1)/2;
        n=n+1;
    elseif f(b1)*f(c1)>0
        b1=c1;
        c1=(a1+b1)/2;
        n=n+1;
    end
end
n
 

f.m文件为自己需要编写的函数

例：用二分法求方程 x^3+4x^2-10=0在区间[1，2]内的根，要求绝对误差不超过1/2*10^(-2)

则f.m文件代码如下：

fun