Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

12.9. Applicazioni di MATLAB® nel Calcolo Numerico
ycn ( i +1)= (2−h)/(2+h) * ycn ( i ) ;
t ( i +1)= t ( i )+h ;
yex ( i +1)=exp(− t ( i + 1 ) ) ;
end
plot ( t , yex , ’k’ , ’linewidth’ ,2 )
hold on %comando per s o v r a s c r i v e r e i g r a f i c i l ’ uno s u l l ’ a l t r o
plot ( t , ye ,’b’ , ’linewidth’ , 2 )
plot ( t , yi ,’r’ ,’linewidth’ , 2 )
plot ( t , ycn , ’g’ ,’linewidth’ , 2 )
legend (’soluzione esatta’ , ’Eulero esplicito’ , ’Eulero implicito’ , ’Crank-Nicolson’ , )
% la legenda , lo spessore d e l l e linee , o i l colore possono e s s e r e
% messi s i a direttamente dalla f i n e s t r a della figura
% s i a dal programma s t e s s o , come in questo caso ( solo in Octave s i
% deve operare direttamente dal programma per poter cambiare
% l e proprieta ’ della figura ) .
hold o f f
for i =1:n+1
f p r i n t f ( fid , ’\n%5.2f %12.6e %12.6e %12.6e %12.6e’ , t ( i ) , yex ( i ) , . . .
ye ( i ) , y i ( i ) , ycn ( i ) ) ;
end
f c l o s e ( f i d ) ;
Questo programma è specifico per l’equazione test assegnata. La sua esecuzioneci permette ci comprendere
meglio il concetto di stabilità dei metodi studiati per la soluzione di equazioni differenziali ordinarie.
203