Bisection script running endlessly

Question

Abdulla Al-Mohannadi 2015-9-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/240490-bisection-script-running-endlessly

回答： Walter Roberson 2015-9-7

I'm trying to develop a script to obtain the root of a function using the bisection method. This is what I have:

function [ TC ] = bisect(xl, xu, es, osf)
format long
xl = xl + 273.15;
xu = xu + 273.15
xr = (xl + xu)/2;
ea = inf; 
ii = 1; 
while ea > es;
 xr_vec = zeros(1,ii) 
    if fTa(xl, osf)*fTa(xr, osf) < 0 
        xu = xr
        ii = (ii-1) + 1; 
        xr = (xl + xu)/2;
        xr_vec(ii) = xr;
        if ii > 1; 
            xr = (xl + xu)/2;
        xr_vec(ii) = xr; 
        ea = (xr_vec(ii) - xr_vec(ii-1))/(xr_vec(ii))*100;
        end
    else xl = xr
    ii = (ii-1)+1;
    xr = (xl + xu)/2;
    xr_vec(ii) = xr;
    if ii > 1;
        xr = (xl + xu)/2;
        xr_vec(ii) = xr;
        ea = (xr_vec(ii) - xr_vec(ii-1))/(xr_vec(ii))*100;
    end
    end
end
TC = xr - 273.15;
end
function f = fTa(Ta, osf)
f = -139.34411 + (1.575701e5)/(Ta) - (6.642308e7)/(Ta^2) + (1.243800e10)/(Ta^3) - (8.621949e11)/(Ta^4) - log(osf);
end

Command window input: bisect(0, 0.4, 0.05, 8)

TC = temperature in Celsius
TA = temperature in Kelvin
osf = oxygen concentration (just another input variable)
xl = specified lower guess
xu = specified upper guess
xr = average of xl and xu
ea = the relative % approximate error
es = specified error boundary

I'm putting all of the xr values in a vector so that ea can be calculated easily. The fTa function at the bottom is being called in the bisect function.

Every time I run it, I keep getting an endless loop with recurring numbers, and I'm not sure where it's going wrong.

I appreciate any help, Thank you