solve the ghk equation using the above constants using the bisection method to an error tolerance of tol and an initial guess of x1bi and x2bi. in matlab

main.m
function result = ghk_eq(x)
    % Constants
    R = 8.314;
    T = 273.15 + 25;
    F = 96485;
    P_K = 10^-2;
    P_Na = 10^-2;
    P_Cl = 10^-2;
    C_Kin = 10^-2;
    C_Nain = 145;
    C_Cl = 110;
    V = -0.070;

    % GHK equation
    result = x - (P_K * C_Kin + P_Na * C_Nain * exp(-x * F / (R * T)) + P_Cl * C_Cl) / (P_K + P_Na * exp(-x * F / (R * T)) + P_Cl);
end
373 chars
17 lines

main.m
function result = bisection_method(tol, x1bi, x2bi)
    % Bisection method implementation
    while abs(x2bi - x1bi) > tol
        xbi = (x1bi + x2bi) / 2;
        if ghk_eq(xbi) == 0
            break;
        elseif ghk_eq(xbi) * ghk_eq(x1bi) < 0
            x2bi = xbi;
        else
            x1bi = xbi;
        end
    end
    result = xbi;
end
352 chars
15 lines

main.m
tol = 1e-6;
x1bi = -0.100;
x2bi = 0.100;

solution = bisection_method(tol, x1bi, x2bi);
disp(solution);
104 chars
7 lines