# Example 5-1
MKS(N)
m1 = 10 kg
fn = m1 * grav()
mus = 0.4
muk = 0.3
#
fa = 0 N
ftr = 0 N
#
fa = 10 N
f_static_max = mus * fn
f_k = muk * fn
ftr = fa
#
fa = 20 N
f_static_max = mus * fn
f_k = muk * fn
ftr = fa
#
fa = 38 N
f_static_max = mus * fn
f_k = muk * fn
ftr = fa
#
fa = 40 N
f_static_max = mus * fn
f_k = muk * fn
ftr = f_k
a = (fa - ftr) / m1

# Example 5-3
MKS(N)
m1 = 10 kg
fg = m1 * grav()
fp = PolarVec(40 N, 30 deg)
fpx = component(fp,1)
fpy = component(fp,2)
fn = fg - fpy
mus = 0.4
muk = 0.3
f_static_max = mus * fn
f_k = muk * fn
ffr = f_k
ax = (fpx - ffr) / m1

# Example 5-5
MKS(N)
ma = 5 kg
fn = ma * grav()
mu = 0.2
ffr = mu * fn
mb = 2 kg
fb = mb * grav()
a = (fb - ffr) / (ma + mb)
Ft = ffr + ma * a

# Example 5-6
MKS(N)
q = 30 deg
mu = tan(q)

# Example 5-7
MKS(N)
ma = 10 kg
q = 37 deg
mus = 0.4
mb_min = ma * sin(q) - mus * ma * cos(q)
mb_max = ma * sin(q) + mus * ma * cos(q)
muk = 0.3
mb = 10 kg
a = (mb * grav() - ma * grav() * sin(q) - muk * ma * grav() * cos(q)) / (ma + mb)

# Example 5-8
m1 = 150 g
r = 0.6 m
om = 2 rps
v = om * r
a = v^2 / r

# Example 5-9
r = 384000 km
tt = 27.3 day
v = 2 * pi() * r / tt
a = v^2 / r
a_rel = a / grav()

# Example 5-10
om = 50000 rpm
r = 6 cm
v = om * r
a = v^2 / r
a_rel = a / grav()

# Example 5-11
MKS(N)
m1 = 0.15 kg
r = 0.6 m
tt = 0.5 s
v = 2 * pi() * r / tt
a = v^2 / r
f = m1 * a

# Example 5-12
MKS(N)
m1 = 0.15 kg
r = 1.1 m
v1 = sqrt(grav() * r)
v2 = 2 * v1
ft2 = m1 * v2^2 / r + m1 * grav()

# Example 5-15
MKS(N)
v = 14 m / s
r = 50 m
q = atan(v^2 / r / grav()); deg

# Example 5-16
v1 = 0
v2 = 35 m/s
tt = 11 s
r = 500 m
a_tan = (v2 - v1) / tt
ar = v2^2 / r

# Homework Problems:

TBD