please solve in matlab

Felix Baumgartner ascended to 39 km in a stratospheric balloon and made a free-fall jump rushing toward earth at supersonic speeds before parachuting to the ground. As he fell, his drag coefficient changed primarily because the air density changed. Recall from Chap. 1 that the terminal velocity, Uterminal (m/s), of a free-falling object can be computed as gm DierminalVc where g gravitational acceleration (ms2), m = mass (kg), and Cd = drag coefficient(kg/m). The drag coefficient can be computed as 0.5pAC where = the fluid density (kg/m3), A = the projected area (m2), and Cd-a dimensionless drag coefficient. Note that the gravitational acceleration, g (m/s2), can be related to elevation by g = 9.80641 2-0.00303 9734z where z - elevation above the earth's surface (km) and the density of air, p (kg/m3), at various elevations can be tabulated as Felix Baumgartner ascended to 39 km in a stratospheric balloon and made a free-fall jump rushing toward earth at supersonic speeds before parachuting to the ground. As he fell, his drag coefficient changed primarily because the air density changed. Recall from Chap. 1 that the terminal velocity, Uterminal (m/s), of a free-falling object can be computed as gm DierminalVc where g gravitational acceleration (ms2), m = mass (kg), and Cd = drag coefficient(kg/m). The drag coefficient can be computed as 0.5pAC where = the fluid density (kg/m3), A = the projected area (m2), and Cd-a dimensionless drag coefficient. Note that the gravitational acceleration, g (m/s2), can be related to elevation by g = 9.80641 2-0.00303 9734z where z - elevation above the earth's surface (km) and the density of air, p (kg/m3), at various elevations can be tabulated as