matlab二分法求函数的根

function root=fun1(fhandle, a, b)
%二分法求方程的根
fa=feval(fhandle, a);
fb=feval(fhandle, b);
root=(a+b)/2;
fc=feval(fhandle, root);
while abs(fc)>0.00001
    if fa*fc>0
        a=root;
    else
        b=root;
    end
    root=(a+b)/2;
    fc=feval(fhandle, root);
end

这里使用matlab中类似c语言函数指针的方法，将代码保存为matlab函数，在命令行输入函数名称后输入定义域区间即可实现求解