Bisection Method2
Bisection Method2
disp('Bisection Method');
y1 = feval(f, low);
y2 = feval(f, high);
i = 0;
if y1 * y2 > 0
m = 'Error'
return
end
i = i + 1;
m = (high + low)/2;
y3 = feval(f, m);
if y3 == 0
return
end
if y1 * y3 > 0
low = m;
y1 = y3;
else
high = m;
end
end
w = feval(f, m);