function dy=gravitacija(t,y,p)
%dy=gravitacija(t,y,p)
%vrne desno stran Newtonovega sistema enacb 
k=4.02*1e+5;%gravitacijska konstanta krat masa Zemlje
r=y(1:2);%trenutni polozaj
v=y(3:4);%trenutna hitrost
dy=[v;-r*k/norm(r)^3];%Newtonova enacba za gravitacijo