erzwungenp.m

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];