erzwungenp.m

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