function [erg] = schussparam(t,y)

global rho0 m quer t0 var cl a vwx vwy z0

grav_c = 6.67390*10^(-11);
m_erde = 5.97223*10^24;
r_erde = 6.3708*10^6;

v = sqrt(y(2)^2 + y(4)^2);
 
% Beschreibung des cw-Wertes in Abhaengigkeit von Geschwindigkeit
if v <= 305
    cd = 0.2*(v/150)^0;
elseif v > 305 & v <= 365
    cd = 0.23*(v/305)^3;
elseif v > 365 & v <= 425
    cd = 0.4*(v/365)^0;
else
    cd = 0.4*(v/425)^(-0.6);
end

%Gravitation abhaengig von Hoehe
g = grav_c*m_erde/((r_erde + y(3))^2);
%Luftwiderstand abhaengig von Temperatur und Hoehe
rho = rho0*(273.15/t0)*exp(-(rho0/100)*g*y(3)/1000);

%Luftwiderstand
Fdx = rho*quer*cd*(y(2)+vwx)^2/(2*m);
Fdz = rho*quer*cd*y(4)^2/(2*m);
%Wind 
Fl = rho*quer*cl*((y(6)+vwy))^2/(2*m);

%Uebergabe an die ODE Funktion
erg = [y(2); -Fdx; y(4); -g + var*Fdz; y(6); -Fl];
   
       
       
       