function [dgl]=ma(t,x)
global a b C ha hb Jw Jz Jges mW mZ mges k r1 r2 s0 tspan ys rho B Z
%Berechnung der Voreinstellung
ha=0.018*(1+sin(100*t)/5); %Bandhöhe links [m]
hb=0.018*(1+cos(100*t)/5); %Bandhöhe rechts [m]
%DGL
dgl=[x(2);((1/2)*sqrt(2)*b*k*(sqrt(2*r1*ha-2*s0*r1-4*r1*x(1)+b*r1*sin(x(3)))-2*x(1)+(1/2)*sqrt(2)*b*k*(sqrt(2*r1*hb-2*s0*r1-4*r1*x(1)-b*r1*sin(x(3)))/mges)));x(4);((2*C*a^2*sin(x(3)))+((1/4)*sqrt(2)*b^2*k*sqrt((2*r1*ha)-(2*s0*r1)-(4*r1*x(1))+(b*r1*sin(x(3))))-((1/4)*sqrt(2)*b^2*k*sqrt((2*r1*hb)-(2*s0*r1)-(4*r1*x(1))-(b*r1*sin(x(3)/Jges))))))];