Beispiel - SAM - ETH Zürich

6 E_tot = E_kin + E_pot;
7
8 f i g u r e(’name’,’Pendulum: energy’);
9 p l o t(tg,E_kin,’b-’,...
0 tg,E_pot,’c-’,...
1 tg,E_tot,’r-’);
2 x l a b e l(’{\bf time t}’,’fontsize’,14);
3 y l a b e l(’{\bf energy}’,’fontsize’,14);
4 legend(’kinetic energy’,’potential energy’,’total
energy’,’location’,’southwest’);
5 t i t l e(’Energies for {\bf implicit midpoint} discrete evolution’);
6
7 i f (nargin > 3)
8 p r i n t(’-depsc2’, s p r i n t f(’%s_EnImid.eps’,filename));
9 end
Numerische
Mathemtik
R. Hiptmair
rev 35327,
25. April
2011
1.4
p. 95