/* Variablendeklaration */
var vy; // Geschw. in y-Richtung
var y; // Ortskoordinate
var t; // Zeit
var v0 = 1; // Startgeschwindigkeit 1 m/s
...
function setup()
{
vy = vy0; // Initialisierung mit Anfangswerten
y = 0;
...
}
...
function draw()
{
...
// Lösung der DGl. 2. Ordnung mittels explicitem EC-Verfahren
var vy_ = vy; // v_i merken
vy = vy - (-g + 2*d*vy + sq(omega0)*(y - l0))*dt; // 1. Integration mit v_i
y = y + vy_*dt; // 2. Integration mit v_i
t = t + dt; // Zeitachse
...
}