schwingNd.m

function dy = schwingNd(t, y, M, C)
% rechte Seite der DGL bei n-dimensionaler Schwingerkette
N = length(y)/2;
x = y(1:N);
v = y((N+1):(2*N));
dx = v;
dv = -inv(M)*C*x;
dy = [dx; dv];