写一段用二分法求二元一次方程的根matlab程序

在MATLAB中，我们可以使用二分搜索法（Binary Search）来近似求解一元二次方程，因为对于一般形式的二次方程ax^2 + bx + c = 0，我们通常无法找到精确的解析解。这里是一个简单的示例程序，用于寻找实数解：

function roots = bisection_method(a, b, c)
    % 检查系数是否满足二次方程条件
    if a == 0
        error('The equation is not quadratic (a must be non-zero).');
    end

    % 初始化搜索区间
    x_lower = -Inf;
    x_upper = Inf;

    % 设置初始精度（例如1e-6）
    tolerance = 1e-6;
    
    % 迭代直到找到解或达到精度限制
    while x_upper - x_lower > tolerance
        mid_point = (x_lower + x_upper) / 2; % 中点
        f_mid = a * mid_point^2 + b * mid_point + c; % 计算函数值

        % 如果f(x_lower) * f(mid_point) < 0，那么根在mid_point和x_lower之间
        if f_mid * c < 0
            x_upper = mid_point;
        % 否则，根在mid_point和x_upper之间
        else
            x_lower = mid_point;
        end
    end

    % 返回近似的根
    roots = [x_lower, x_upper];
end

运行此函数，你可以传入二次方程的系数a、b和c，比如`roots = bisection_method(1, -3, 2)`。注意，这个算法返回的是两个可能的根，因为二次方程可能有两个实根。