Beispiel - SAM - ETH Zürich

7
8 i f (nargin > 3),
p r i n t(’-depsc2’, s p r i n t f(’%s_orbit.eps’,filename)); end
9
0 % Tracking of energies
1 E_kin = 0.5*(y_p.^2);
2 E_pot = -g/l*cos(y_sv);
3 E_pot = E_pot - min(E_pot) + min(E_kin);
4 E_tot = E_kin + E_pot;
5
6 f i g u r e(’name’,’Pendulum: energy’);
7 p l o t(tg,E_kin,’b-’,tg,E_pot,’c-’,tg,E_tot,’r-’);
8 x l a b e l(’{\bf time t}’,’fontsize’,14);
9 y l a b e l(’{\bf energy}’,’fontsize’,14);
0 legend(’kinetic energy’,’potential energy’,’total
energy’,’location’,’southeast’);
1 t i t l e(’Energies for {\bf Stoermer-Verlet} discrete evolution’);
2 i f (nargin > 3),
p r i n t(’-depsc2’, s p r i n t f(’%s_EnSV.eps’,filename)); end
Numerische
Mathemtik
R. Hiptmair
rev 35327,
25. April
2011
1.4
p. 110