%% Biseccion

close all
clc
clear

% f = @(x) x.^2 - 3;
% f = @(x) x.^3 - x - 1; 
f = @(x) cos(5*x) - x.^3 + 2; 

f(1)
f(2)

a = 1;
b = 2;
n_biseccion = 0;
x_biseccion = [];

while b-a > 1e-3
    plot(linspace(a-0.1,b+0.1), f(linspace(a-0.1,b+0.1)), 'Linewidth', 2)
    hold on
    plot(linspace(1,2), 0*linspace(1,2), '--k')
    plot(a,f(a),'ok', 'MarkerSize', 6, 'MarkerFaceColor', 'k')
    plot(b,f(b),'ok', 'MarkerSize', 6, 'MarkerFaceColor', 'k')
    pause
    close
    x = (a + b)/2;
    if f(x)*f(b) > 0
        b = x;
    else
        a = x;
    end
    n_biseccion = n_biseccion + 1;
    x_biseccion = [x_biseccion; x];
end

n_biseccion
x_biseccion