$$\begin{gathered} \mathrm{restart}\: \\ \mathrm{with}\left(\mathrm{ThermophysicalData}\right) \end{gathered}$$

$\left[\mathrm{Atmosphere}\,\mathrm{Chemicals}\,\mathrm{CoolProp}\,\mathrm{PHTChart}\,\mathrm{Property}\,\mathrm{PsychrometricChart}\,\mathrm{TemperatureEntropyChart}\right]$

$\mathrm{Atmosphere}\left(3⟦\mathrm{km}⟧\,\mathrm{all}\,\mathrm{altitudetype}=\mathrm{geopotential}\right)$

$268.6500⟦K⟧,70108.52215⟦\mathrm{Pa}⟧,0.9091219586⟦\frac{\mathrm{kg}}{m^3}⟧,328.5779005⟦\frac{m}{s}⟧,0.00001693718731⟦\mathrm{Pa}s⟧$

$\mathrm{Atmosphere}\left(3⟦\mathrm{km}⟧\,\mathrm{all}\,\mathrm{altitudetype}=\mathrm{geometric}\right)$

$268.6591984⟦K⟧,70121.13970⟦\mathrm{Pa}⟧,0.9092544424⟦\frac{\mathrm{kg}}{m^3}⟧,328.5835256⟦\frac{m}{s}⟧,0.00001693764616⟦\mathrm{Pa}s⟧$

$\mathrm{Atmosphere}\left(3000\,\mathrm{all}\,\mathrm{altitudetype}=\mathrm{geometric}\right)$

$268.6591984,70121.13970,0.9092544424,328.5835256,0.00001693764616$

$\mathrm{plotOpts}≔\mathrm{size}=\left[400\,400\right],\mathrm{background}=''LightGrey'',\mathrm{axis}=\left[\mathrm{gridlines}=\left[10\,\mathrm{color}=\mathrm{white}\right]\right],\mathrm{titlefont}=\left[\mathrm{Arial}\,16\right],\mathrm{axesfont}=\left[\mathrm{Arial}\right],\mathrm{labelfont}=\left[\mathrm{Arial}\,12\right],\mathrm{axes}=\mathrm{boxed},\mathrm{labeldirections}=\left[\mathrm{horizontal}\,\mathrm{vertical}\right]\:$

$\mathrm{de}≔m \frac{{\ⅆ}^2}{\ⅆt^2}y\left(t\right)=-m\cdot g\+\frac{1}{2}A\cdot \mathrm{Cd}\cdot \mathrm{\ρ}\left(t\right)\cdot \mathrm{vel}{\left(t\right)}^2,\mathrm{vel}\left(t\right)=-\frac{\ⅆ}{\ⅆt}y\left(t\right)\:$

m is mass

g is acceleration due to gravity

Cd is the drag coefficient

ρ(t) is the air density

$\mathrm{\ρ}≔t\to \mathrm{Atmosphere}:-\mathrm{Property}\left(''density''\,y\left(t\right)\right)\:$

$A≔1\:$$m≔85\:$$\mathrm{Cd}≔0.58\:$$g≔9.8\:$

$\mathrm{sol}≔\mathrm{dsolve}\left(\left\{\mathrm{de}\,\mathrm{vel}\left(0\right)=0\,y\left(0\right)=22000\,\mathrm{D}\left(y\right)\left(0\right)=0\right\}\,\mathrm{numeric}\,\mathrm{known} \= \left\{\mathrm{Atmosphere}:-\mathrm{Property}\right\}\, \mathrm{range}\=0..250\right)$

$\mathrm{sol}≔\mathbf{proc}\left(\mathrm{x\_rkf45\_dae}\right)...\mathbf{end\; proc}$

$\mathrm{plots}:-\mathrm{odeplot}\left(\mathrm{sol}\,\left[t\,y\left(t\right)\right]\,t=0..250\,\mathrm{labels}=\left[''time\; (s)''\,''Altitude\; (m)''\right]\,\mathrm{plotOpts}\right)$

$\mathrm{plots}:-\mathrm{odeplot}\left(\mathrm{sol}\,\left[t\,\mathrm{vel}\left(t\right)\right]\,t=0..250\,\mathrm{plotOpts}\right)$

$\mathrm{plots}:-\mathrm{odeplot}\left(\mathrm{sol}\,\left[\mathrm{vel}\left(t\right)\,y\left(t\right)\right]\,t=0..250\,\mathrm{plotOpts}\right)$