matlab 二分法求根

二分法是一种求解非线性方程的数值方法，也称为折半法。在MATLAB中，可以通过编写函数来实现二分法求根。下面是一个示例函数：

function [x0,i] = dichotomy(a,b,eps)
a0 = a;
b0 = b;
x0 = (a0+b0)/2;
for i = 1:1000
    if abs((b-a)/(2^(i))) > eps
        if f(a0)*f(x0) < 0
            b0 = x0;
            x0 = (a0+b0)/2;
        elseif f(a0)*f(x0) > 0
            a0 = x0;
            x0 = (a0+b0)/2;
        elseif f(a0)*f(x0) == 0
            break
        end
    else
        break
    end
end

其中，a和b分别为区间的左右端点，eps为精度要求，f为待求解的非线性方程。函数返回近似根x0和二分次数i。