Beispiel - SAM - ETH Zürich

Berechney 1 aus (1.4.30) & (1.4.31)
Beispiel 1.4.32 (Störmer-Verlet-Verfahren für Pendelgleichung).
△
Numerische
Mathemtik
Listing 1.6: Störmer-Verlet-Verfahren für Pendelgleichung
1 f u n c t i o n sverletpend(y0,v0,T,N,filename)
2 % MATLAB function for Stoermer-Verlet simulation of movement of mathematical
3 % pendulum, Example 1.4.32
4 %
5 % y0,v0: initial values for angle and its temporal derivative
6 % T : final time
7 % N : number of timesteps
8
9 g = 9.8; l = 1;
0
1 % right hand side for (1.4.25) (Newton’s equations of motion)
2 f = @(y) -g/l*sin(y);
3
4 h = T/N; % uniform timestep
5 y_old = y0; y_new = h*v0+y0+ 0.5*h*h*f(y0); % initial step, see
Rem. 1.4.29
R. Hiptmair
rev 35327,
25. April
2011
1.4
p. 107