function dy = erzwungenp(t, y, m, b, c, A, omega) % rechte Seite der DGL bei erzwungener Schwingung % Parameter: % m Masse % b Reibungskoeffizient % c Federkonstante % A Amplitude der äußeren Kraft % omega Kreisfrequenz der äußeren Kraft dy1 = y(2); dy2 = -b/m*y(2) - c/m*y(1) + A/m*cos(omega*t); dy = [dy1; dy2];