// 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 implicitem EC-Verfahren
vy = vy - (-g + 2*d*vy + sq(omega0)*(y - l0))*dt; // 1. Integration
y = y + vy*dt; // 2. Integration mit v_i+1
t = t + dt; // Zeitachse
...
}