% main per newton
close all
clear all
clc


% definizione della funzione nel file miafunz.m
% definizione della derivata della funzione nel file miadfunz.m

% APPROSSIMAZIONE DELLA PRIMA RADICE
disp('PRIMA RADICE');
x0 = 0.1;
toll = 1e-8;
nmax = 1000;

% richiamo bisezione
[xvect,xdif,fx,nit] = newton(x0,toll,nmax,'miafunz','miadfunz');

% stampo soluzione
disp(['Vettore soluzione: ']);
disp(xvect(end));
disp(['Valore dello zero: ']);
disp(miafunz(xvect(end)));

% numero di iterazioni eseguite
disp(['Numero di iterazione eseguite: ',num2str(nit)]);



% APPROSSIMAZIONE DELLA SECONDA RADICE
disp(' ');disp(' ');disp(' ');
disp('SECONDA RADICE');

x0 = 3;
toll = 1e-8;
nmax = 1000;

% richiamo bisezione
[xvect,xdif,fx,nit] = newton(x0,toll,nmax,'miafunz','miadfunz');

% stampo soluzione
disp(['Vettore soluzione: ']);
disp(xvect(end));
disp(['Valore dello zero: ']);
disp(miafunz(xvect(end)));

% numero di iterazioni eseguite
disp(['Numero di iterazione eseguite: ',num2str(nit)]);
