Although this is a solution in Python (not in Matlab)
It produced the following graphs (which make sense due to the increased aerodynamic friction in the lower layers of the atmoshpere.
Additionally the code below produces also the phase portrait which is:
the following code creates the above:
from scipy.integrate import odeint
import numpy as np