% main per bisezione


% definizione della funzione nel file miafunz.m

% APPROSSIMAZIONE DELLA PRIMA RADICE
a = 0.13;
b = 0.14;
nmax = 1000;
toll = 1e-8;

% richiamo bisezione
[xvect,xdif,fx,nit] = bisezione(a,b,toll,nmax,'miafunz');

% stampo soluzione
xvect(end)

% numero di iterazioni eseguite
nit

% disegniamo errore
figure(1);
semilogy(abs(xvect(1:end-1)-xvect(end)));

% APPROSSIMAZIONE DELLA SECONDA RADICE
a = 1.5;
b = 1.6;
nmax = 1000;
toll = 1e-8;

% richiamo bisezione
[xvect,xdif,fx,nit] = bisezione(a,b,toll,nmax,'miafunz');

% stampo soluzione
xvect(end)

% numero di iterazioni eseguite
nit

% disegniamo errore
figure(2);
semilogy(abs(xvect(1:end-1)-xvect(end)));
