function du=oszitotal_dgl(t,u)            %definieren der hom. DGL
global m d rho                                                             %globale Variablen der DGL

du=[u(2);(-d*(u(1)-u(3))-d*(u(1)-u(23))-u(2)*rho)/m;...
    u(4);(-d*(u(3)-u(5))-d*(u(3)-u(1))-u(4)*rho)/m;...
    u(6);(-d*(u(5)-u(7))-d*(u(5)-u(3))-u(6)*rho)/m;...
    u(8);(-d*(u(7)-u(9))-d*(u(7)-u(5))-u(8)*rho)/m;...
    u(10);(-d*(u(9)-u(11))-d*(u(9)-u(7))-u(10)*rho)/m;...
    u(12);(-d*(u(11)-u(13))-d*(u(11)-u(9))-u(12)*rho)/m;...
    u(14);(-d*(u(13)-u(15))-d*(u(13)-u(11))-u(14)*rho)/m;...
    u(16);(-d*(u(15)-u(17))-d*(u(15)-u(13))-u(16)*rho)/m;...
    u(18);(-d*(u(17)-u(19))-d*(u(17)-u(15))-u(18)*rho)/m;...
    u(20);(-d*(u(19)-u(21))-d*(u(19)-u(17))-u(20)*rho)/m;...
    u(22);(-d*(u(21)-u(23))-d*(u(21)-u(19))-u(22)*rho)/m;...
    u(24);(-d*(u(23)-u(1))-d*(u(23)-u(21))-u(24)*rho)/m];               %hom. DGL in du-Vektor