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