# Example 18.1
m0 = 250e3 kg
tt = 700 kN
ax = tt / m0
BEGIN
  a = 1 kN
  b = 1 kN
  fy = a + b - m0 * grav()
  M = 2 m * tt + 22 m * b - 5 m * a
END
solve(fy,M,a,b)


# Example 18.2
m0 = 45 kg
b = 0.15 m
mu = 0.4
ia = 4 kg m^2
MM = 50 N m
q = 20 deg
BEGIN
  alp = 1 rad/s^2
  ax = b * alp
   err = MM - b * m0 * grav() * (sin(q) + mu * cos(q)) - b * m0 * ax - ia * alp
END
solve(err, alp)


# Example 18.12
xbar = 0.3 m
mbar = 3 kg
xdisk = (0.6 + 0.2) m
mdisk = 2 kg
x = (xbar * mbar + xdisk * mdisk) / (mbar + mdisk)
ibar = (1/12) * mbar * (0.6 m)^2 + mbar * (0.2 m)^2
idisk = 0.5 * mdisk * (0.2 m)^2 + mdisk * (0.3 m)^2
ii = ibar + idisk



# Homework Problems:

TBD