function rdotdot = rzut(t,r)
   Fop = -0.5*cx*rho*S*[r(2)^2,r(4)^2];
   rdotdot = [r(2);
              -Fop(1)/m;
              r(4);
              -g-Fop(2)/m];
endfunction
//Parametry modelu
m = 0.058; r = 0.0335; S = %pi*r^2; cx = 0.45; rho = 1.225;
V0 = 170/3.6; alpha = 45; wp =