function xdot = vdpol(t,x)
%xdot = [-9.8/(1+x(2)/6370000)^2 - 0.5*exp(-x(2)/8000).*(x(1)*abs(x(1)))/70;x(1)];

xdot = [-9.8 - 0.5*exp(-x(2)/8000).*(x(1)*abs(x(1)))/70;x(1)];

%xdot = [-9.8 - 0.5*(x(1).*abs(x(1)))/70;x(1)];

%xdot = [-9.8;x(1)];


%xdot = [3*exp(-x(2));x(1)];
