function dy = erzwungen(t, y) % rechte Seite der DGL bei erzwungener Schwingung % Parameter m = 1.0; b = 0.3; c = 1.0; A = 1.0; omega = 1.0; dy1 = y(2); dy2 = -b/m*y(2) - c/m*y(1) + A/m*cos(omega*t); dy = [dy1; dy2];