# Example 16.1
m0 = 1200 kg
t = 10 s
ix = (720 N/s) * t^2 / 2
iy = 2160 N * t - (360 N / s) * t^2 / 2
iz = 0 N s
vx = ix / m0
vy = iy / m0
vz = iz / m0
v1 = [vx, vy, vz]
t = 20 s
v = [36, 8] m/s
fav = m0 * (v - [vx, vy]) / 10 s

# Example 16.2
t = 30 
v = (200 * t - 3 * t^2) / 150
f = 150 * v / 30


# Example 16.3
IPS(deg)
v = 1.9 in / 0.001 s
m0 = 1.62 ozm
d = 1.68 in
q = 21 deg
e = [cos(q), sin(q)]
f = m0 * v * e / (0.0006 s)


# Example 16.5
ma = 4 kg
va = 10 m/s
mb = 4 kg
vb = -5 m/s
vp = (ma * va + mb * vb) / (ma + mb)
BEGIN
   vap = 1 m/s
   vbp = 1 m/s
   e1 = (ma * va + mb * vb) - (ma * vap + mb * vbp)
   e2  = (vbp - vap) / (va - vb) - 0.8
END
solve(e1,e2,vap,vbp)


# Example 16.6
ma = 18e3 kg
mb = 6.6e3 kg
va = [0.2, 0.03, -0.02] m/s
vb = [0, 0, 0] m/s
v = (ma * va + mb * vb) / (ma + mb)
BEGIN
   vap = 1 m/s
   vbp = 1 m/s
   e1 = (ma * va[1] + mb * vb[1]) - (ma * vap + mb * vbp)
   e2  = (vbp - vap) / (va[1] - vb[1]) - 0.95
END
solve(e1,e2,vap,vbp)

# Example 16.8
vp = 7000 m/s
rp = 10000 km
re = 6370 km
va = 2 * grav() * re^2 / (rp * vp) - vp
ra = rp * vp / va


# Example 16.9
rho = 1000 kg/m^3
v = 300 m/s
a = 0.01 m^2
f = -rho * v^2 * a


# Example 16.12
dmcdt = 13.5 kg/s
dmfdt = 0.130 kg/s
ve = 480 m/s
vi = 120 m/s
TT = (dmcdt + dmfdt) * ve - dmcdt * vi


# Homework Problems:

TBD