function neu()
global m a b c g my_G A B C D E F G H I J K L M;
m=750; a=3.9; b=3.6; c=2.1; g=9.81; my_G=0.07;
alpha0=[1.047197551,0];
tspan=[0 1];
[t,alphaout]=ode45(@bewegungsgleichung,tspan,alpha0);
gamma=-((pi/2)-alphaout(:,1));
alpha_punkt=alphaout(:,2);
figure;plot(t,gamma), grid on, title('Winkel gamma')
figure;plot(t,alpha_punkt),grid on, title('Winkelgeschwindigkeit')
function [alphapunkt]=bewegungsgleichung(t,alphaout)
beta=(asin((b-(a*(sin(alphaout(1)))))/c));
A=(((a/2)*sin(alphaout(1)))+(a*cos(alphaout(1))));
B=(((a/2)*cos(alphaout(1)))-(tan(beta)*a*sin(alphaout(1)))-((((tan(beta))^2)*(a^2)*cos(alphaout(1)))/(c*cos(beta)))-((a^2)*((cos(alphaout(1)))^2)));
C=((a/2)*cos(alphaout(1)));
D=((a/2)*sin(alphaout(1)));
E=((m*my_G)/(cos(beta)+(my_G*sin(beta))));
F=(((m*A)/my_G)+((E*cos(beta)*C)/my_G)+((E*A*cos(beta))/(my_G^2)));
G=(((m*B)/my_G)+((E*cos(beta)*B)/(my_G^2))-((E*D*cos(beta))/my_G));
H=((E*g*cos(beta))/my_G);
I=((m*(a^2))/12);
J=(((a/2)*sin(beta)*cos(alphaout(1)))+((a/2)*cos(beta)*sin(alphaout(1))));
K=(-I-(J*E*(C+(A/my_G)))-((a/2)*sin(alphaout(1))*my_G*F)+((a/2)*cos(alphaout(1))*F));
L=((J*E*((B/my_G)-D))+((a/2)*sin(alphaout(1))*my_G*G)-((a/2)*cos(alphaout(1))*G))
M=((E*g)+((a/2)*sin(alphaout(1))*my_G*H)-((a/2)*cos(alphaout(1))*F));
alphapunkt=[alphaout(2);(((L*((alphaout(2))^2))+M)/K)]
end

end