erzwungen.m

function dydt = 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);
dydt = [dy1 dy2]';