Beispiel - SAM - ETH Zürich

8 x l a b e l(’{\bf l}’,’fontsize’,14);
9 y l a b e l(’{\bf p}’,’fontsize’,14);
0 t i t l e( s p r i n t f(’Phase flow for Zeeman model (\\alpha =
%d,\\beta=%d)’,...
1 alpha,beta));
Numerische
Mathemtik
2
3 % Compute evolution of l (length) and p (potential), see Rem. 1.2.4
4 tspan = [0 100]; % Duration of simulation
5 [t,y] = ode45(odefun,tspan,[1;0],odeset(’abstol’,1E-12));
6
7 % Plot trajectory of solution
8 p l o t(1,0,’k*’,’markersize’,10);
9 p l o t(y(:,1),y(:,2),’b-’);
0 hold off;
1 p r i n t(’-depsc2’, s p r i n t f(’%s1.eps’,filename));
2
3 % Plot time-dependent solution
4 f i g u r e(’name’,’heartbeat’);
5 p l o t(t,y(:,1),’r-’,t,y(:,2),’b-’);
6 t i t l e( s p r i n t f(’heartbeat according to Zeeman model (\\alpha =
%d,\\beta=%d)’,alpha,beta));
7 x l a b e l(’{\bf time t}’,’fontsize’,14);
8 y l a b e l(’{\bf l/p}’,’fontsize’,14);
R. Hiptmair
rev 35327,
25. April
2011
1.2
p. 36