+ Problem 1. The Basic Equations.

- If we let positive velocity be in the downward direction, then the two forces acting on the sky diver are:

force due to gravity, g, is F_g = m g;
force due to air resistance is F_r = - k v;

where m g = is the weight of the sky diver (and chute) and k is the coefficient of friction.

m (dv/dt) = F_g + F_r

that is

dv/dt = g - (k/m) v.

- First, we solve the initial value problem with arbitrary k.

Input := 

k = .
g = 32;
m = 155/32;
de1 = v'[t] == g - k/m v[t];
v0 = 0;
Input := 

ans = Flatten[DSolve[{de1, v[0]==v0}, v[t], t]]
Output =

         155         155
{v[t] -> --- - ---------------}
          k     (32 k t)/155
               E             k
Input := 

v[t] /. ans
Output =

155         155
--- - ---------------
 k     (32 k t)/155
      E             k
Input := 

LimVel = %[[1]]
Output =

155
---
 k

- Now, we can determine the physical drag coefficients for each of the limiting velocities. This then yields the three different (differential) equations.

Input := 

temp = Flatten[Solve[LimVel==278.8, k]];
k1 = k /. temp
Output =

0.555954
Input := 

temp = Flatten[Solve[LimVel==180.4, k]];
k2 = k /. temp
Output =

0.859202
Input := 

temp = Flatten[Solve[LimVel==21.32, k]];
k3 = k /. temp
Output =

7.27017
Input := 

de1 = v'[t] == g - k1 v[t]/m;
de2 = v'[t] == g - k2 v[t]/m;
de3 = v'[t] == g - k3 v[t]/m;