schwing2d.m

function dy = schwing2d(t, y)
% rechte Seite der DGL bei zweidimensionaler Schwingerkette

% Parameter
m1 = 1.0;
m2 = 1.0;
c1 = 4.0;
c2 = 0.6;
c3 = 4.0;

dy1 = y(3);
dy2 = y(4);
dy3 = (1/m1)*(-(c1+c2)*y(1) + c2*y(2));
dy4 = (1/m2)*(c2*y(1) - (c2+c3)*y(2));

dy = [dy1; dy2; dy3; dy4];