Beispiel - SAM - ETH Zürich

%f,\\alpha(0)=%f,p(0)=%f’,g,l,y0,v0));
9 i f (nargin > 3),
p r i n t(’-depsc2’, s p r i n t f(’%s_alpha.eps’,filename)); end
0
1 % Plot of velocity vs. time
2 tg = h*(0:N);
3 f i g u r e(’name’,’Pendulum velocity’);
4 p l o t(t,y(:,2),’g-’,tg,y_p,’r-+’);
5 x l a b e l(’{\bf time t}’,’fontsize’,14);
6 y l a b e l(’{\bf velocity p}’,’fontsize’,14);
7 t i t l e( s p r i n t f(’Pendulum g = %f, l =
%f,\\alpha(0)=%f,p(0)=%f’,g,l,y0,v0));
8 i f (nargin > 3), p r i n t(’-depsc2’, s p r i n t f(’%s_p.eps’,filename));
end
9
0 % PLot of trajectory in phase space
1 f i g u r e(’name’,’Pendulum trajectory’);
2 ph = p l o t(y(:,1),y(:,2),’g--’,y_sv,y_p,’r-+’);
3 set(ph(1),’linewidth’,2);
4 x l a b e l(’{\bf angle \alpha}’,’fontsize’,14);
5 y l a b e l(’{\bf velocity p}’,’fontsize’,14);
6 t i t l e( s p r i n t f(’Pendulum g = %f, l =
%f,\\alpha(0)=%f,p(0)=%f’,g,l,y0,v0));
Numerische
Mathemtik
R. Hiptmair
rev 35327,
25. April
2011
1.4
p. 109