%% 注意输入的参数。函数fun以及给定的区间[a,b]，求出精确度为ep的零点x_star，k为迭代次数 ;
function [x_star,k] = bisect1(fun,a,b,ep)
%二分法求方程的近似值
if nargin<4 % nargin,nargout 是函数文件里的固有变量，意思是输入（输出）参数的个数
    ep=1e-5; % 如果bisect1只输入前三个参数，那么精确度默认为1e-5
end
fa=feval(fun,a);fb=feval(fun,b); %函数fun 在点a,b 处的函数值
if(fa*fb>0)  %如果a,b之间不存在零解
    disp('error');
    x_star=NaN;
    k=NaN;
else
    k=1; % 迭代次数初始化
    while abs((b-a)/2)>ep %区间大于精确度的话就迭代
        x=(a+b)/2;
        fx=feval(fun,x);
        if fa*fx<0 %[a,x]区间存在零解
            b=x;fb=fx; %缩小b的上限
        else
            a=x;fa=fx; %如果无解，加大a的下限
        end
        k=k+1; %每循环一次 迭代加一
    end
    x_star=(a+b)/2; %把最后的结果赋值给输出的x_star
end

fun=inline('x^3-x-1','x'); %定义函数
[x_star,k]=bisect1(fun,1,15,0.0000005)