MATLAB.rar_Factorial_interpolation

%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