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