Re: Ordine di convergenza di metodi risolutivi di ODE a *k* passi

Messaggioda amatrix » 22/07/2019, 20:19

Ho provato a inserire anche i seguenti metodi di ordine 4 (AM e AB) e ricontrollato svariate volte i coefficienti ma dal grafico convergono come O(h)...

Codice:
%AB4
for i=5:ts
u(1)=y0;
u(2)=uex(1);
u(3)=uex(2);
u(4)=uex(3);
u(i)=u(i-1)*(1-h*55/24)+u(i-2)*(h*59/24)-u(i-3)*h*37/24+u(i-4)*h*9/24 ;
end


e

Codice:
%AM3
for i=4:ts
b(1)=y0;
b(2)=uex(1);
b(3)=uex(2);
b(i)=b(i-1)*(1-19*h/24)/(1+h*9/24)-b(i-2)*(-h*5/24)/(1+h*9/24)+b(i-3)*(-h/24)/(1+h*9/24);
end


E' possibile che definendoli attraverso il ciclo for incombo automaticamente in degli errori o che magari questi dipendano dalle condizioni iniziali benché stia inserendo come dati iniziali direttamente i valori della funzione?
Immagine
amatrix
Starting Member
Starting Member
 
Messaggio: 6 di 6
Iscritto il: 11/07/2019, 16:02

Precedente

Torna a Analisi Numerica e Ricerca Operativa

Chi c’è in linea

Visitano il forum: Nessuno e 2 ospiti