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