# Example 10-1
3e-4 rad; deg

# Example 10-3
a = 0.06 rad/s^2
del_t = 8 s
om = a * del_t
r = 2.5 m
v = om * r; m/s
at = r * a
ar = v^2 / r
a_vec = [ar, at]
mag(a_vec)
PolarAngle(a_vec); deg

# Example 10-4
om = 7200 rpm; rad/s
r = 3 cm
v = om * r
x = 0.5e-6 m
del_t = v / x

# Example 10-5
r = 3 m
a = 1.6 rad/s
b = 1.2 rad/s^2
tt = 2 s
om = a + b * tt
v = om * r
a = [r * b, om^2 * r]

# Example 10-6
om = 20000 rpm
tt = 30 s
a = om / tt
q = 0.5 * a * tt^2; rev

# Example 10-7
ra = 30 cm
rb = 50 cm
fa = 50 N
fb = 50 N
q = 30 deg
tor = ra * fa - rb * fb * cos(q); N m

# Example 10-8
m1 = 5 kg
m2 = 7 kg
r = 2 m
ii = m1 * r^2 + m2 * r^2
x1 = 0.5 m
x2 = 4.5m
ii = m1 * x1^2 + m2 * x2^2

# Example 10-9
Ft = 15 N
m1 = 4 kg
r0 = 33 cm
om = 30 rad/s
tt = 3 s
tfric = 1.1 N m
tor = r0 * Ft - tfric; N m
a = om / tt
ii = tor / a

# Example 10-10
w = 15 N
m0 = w / grav()
r0 = 33 cm
tfric = 1.1 N m
ii = 0.385 kg m^2
a = (w * r0 - tfric) / (ii + m0 * r0^2)
aa = r0 * a

# Example 10-10
w = 15 N
m0 = w / grav()
r0 = 33 cm
tfric = 1.1 N m
ii = 0.385 kg m^2
a = (w * r0 - tfric) / (ii + m0 * r0^2)
aa = r0 * a

# Example 10-21
a = 0.5 * grav()
m1 = 1200 kg
sum_f = m1 * a
x1 = 1.5 m
x2 = 0.75 m
mu =0.3
BEGIN
fn1 = 5000 N
fn2 = 5000 N
mu = m1 * a / (fn1 + fn2)
eq1 = x1 * (fn1 - fn2) - x2 * mu * (fn1 + fn2)
eq2 = fn1 + fn2 - m1 * grav()
END
solve(eq1,eq2,fn1,fn2)

# Homework Problems:

TBD