%funzione analizzata f(x)=4*x^2 - 30*x + 5/6
clear all
clc
%faccio sceglie due valori x1 ed x2 arbitrariamente
%ovviamente il programma è ancora da rifinire poiché il programma non funziona con qualsiasi intervallo
x1 = input('scegli il primo valore dell''intervallo ');
x2 = input('scegli il secondo valore dell''intervallo ');
%calcolo la f(x) dei rispettivi valori
y1 = 4*(x1)^2 - 30*x1 + 5/6;
y2 = 4*(x2)^2 - 30*x2 + 5/6;
%inizializzo la variabile di controllo (y3) solo per evitare che ci siano problemi per entrare nel ciclo while
y3=100;
%variabile di controllo ulteriore
control=0;
%entrerà nel ciclo while ogni volta che y3 è diverso da 0 e che il prodotto delle due y sia negativo, e finché la variabile di controllo non superi le 1000 iterazioni
%nel momento in cui una delle tre condizioni risulti falsa (y3=0, il prodotto y1*y2 >=0 o che control>=1000) esce dal ciclo
while ((y3~=0) & (y1*y2)<0) & (control<1000)
%calcolo il valore medio e la sua ordinata
x3 = (x1 + x2)/2;
y3 = 4*(x3)^2 - 30*x3 + 5/6;
%il controllo sull'intervallo se y3 è positivo il suo valore diverrà il nuovo estremo destro dell'intervallo (cioè y2)
%altrimenti se negativo il valore di y3 sarà assegnato ad y1 (estremo sinistro)
if (y3>0) & (x3>0)
x2=x3;
y2=y3;
elseif (y3<0) & (x3>0)
x1=x3;
y1=y3;
end
if (y3>0) & (x3<0)
x1=x3;
y1=y3;
elseif (y3<0) & (x3<0)
x2=x3;
y2=y3;
end
control=control + 1;
end
%stampa a video
if x1==x2
fprintf('L''equazione analizzata è f(x)=4x^2 - 30x + 5/6');
fprintf('Il valor medio è approssimato al valore ');
disp(x3);
else
fprintf('Non è stato trovato alcun punto d''intersezione');
end
