chapter 02 problems

prob02.sage

print("Problem 3, 12, 18, 23, 44, 46, 66")

print("Problem 3:")
print("(a) distance traveled:", abs(10-2) + abs(8-10) + abs(11-8), " [m]")

print("(b) magnitude of the displacement:", abs(11-2), " [m]")

# (2.1)
print("(c) displacement:", 11-2, " [m].")


print("Problem 12:")
l = 1.1*units.length.meter
v = 18*(units.length.meter/units.time.second)

time = l/v
print(time)

print("Problem 18:")
a = 1.40 * (units.length.meter/(units.time.second)^2)
t = var("t")
# sol=solve(a*(t*units.time.second) == 2.00*(units.length.meter/units.time.second), [t])
sol=solve(a*t == 2.00*(units.length.meter/units.time.second), [t])
s= sol[0]
s = 1.0*s.rhs()
print(s)

a = var("a")
sol = solve(a*(0.800*units.time.second) == -2.00*(units.length.meter/units.time.second), [a])
s = sol[0]
s = 1.0*s.rhs()
print(s)

print("Problem 23:")
a = 1.35*(units.length.meter/(units.time.second)^2)
v = 80.0*(units.length.kilometer/units.time.hour)
t = var("t")
sol = solve( a*t == v, [t])
s = sol[0]
s = 1.0*s.rhs()
print(s.convert(units.time.second), s.convert(units.time.minute))

da = -1.65*(units.length.meter/(units.time.second)^2)
t = var("t")
sol = solve( v + da*t == 0, [t])
s = sol[0]
s = 1.0*s.rhs()
print(s.convert(units.time.second), s.convert(units.time.minute))

et = 8.30*units.time.second
eda = var("eda")
sol = solve( v + eda*et == 0, [eda])
s = sol[0]
s = 1.0*s.rhs()
print(s.convert((units.length.meter/(units.time.second)^2)))


print("Problem 44:")
g = 9.80*(units.length.meter/(units.time.second)^2)
vi = 1.40*(units.length.meter/units.time.second)
t = 1.8*units.time.second
print("(a): ", g, vi, t)

print("(b): ") 
h = var("h")

sol = solve( h == vi*t + g/2*t^2, [h])
s = sol[0]
s = 1.0*s.rhs()
print(s, n(s/units.length.meter, digits=3))


print("Problem 46:")
vi = 4.00*(units.length.meter/units.time.second)
hi = 1.80*units.length.meter

t = var("t")
sol = solve( 0 == hi + vi*t + (-g)*t^2/2, [t])
t0 = sol[1]
t0 = t0.rhs()
print(t0.convert())
print(n( t0/units.time.second, digits=3))

t = var("t")
sol = solve( 0 == vi + (-g)*t, [t])
s = sol[0].rhs()
print(1.0*s)
print( vi*s + (-g)*s^2/2 )

vf = vi + (-g)*t0 
vf = 1.0*vf
print(vf)
print(n(vf/(units.length.meter/units.time.second), digits=3))


print("Problem 66:")

reset()

print("prob 04_06\n")
decel = 196*units.length.meter/(units.time.second)**Integer(2)
mass = Integer(2100)units.mass.kilogram
f = var("f")
s =solve ( f == massdecel, f )
s = s[Integer(0)].rhs()
s = s.convert(units.force.newton)
print("the force is: ", s )

print("prob 04_13\n")
weight_on_moon = 250units.force.newton
g_moon = 1.625units.length.meter/(units.time.second)^2
m = var("m")
sol = solve( m*g_moon == weight_on_moon, m)
s =sol[0].rhs()
s = s.convert(units.mass.kilogram)
print("the mass is: ", s)

mass = s/units.mass.kilogram
print("i.e., ", mass.n(digits=3))

print(massunits.mass.kilogram g_moon)

print("prob 04_20\n")
mass = 60.0units.mass.kilogram
g = 9.80units.length.meter/(units.time.second)^2
f = mass*g
print("(a): ", f.convert(units.force.newton))

a = 1.50units.length.meter/(units.time.second)^2
f2 = mass(g+a)
print("(b): ", f2.convert(units.force.newton))

print("prob 04_22\n")
weight = 55*units.force.newton
mass_total = weight/g
print("(a): ", mass_total.convert())

t1, t2 = var("t1 t2")
eq1 = t1==weight
eq2 = t2==t1+ (0.500*units.mass.kilogram)*g
sol = solve([eq1,eq2], [t1,t2])
t1 = t1.subs(sol)
print("(b): ", t1)
t2 = t2.subs(sol)
t2 = t2.convert(units.force.newton)
print("(c): ", t2, " i.e., ", n(t2/units.force.newton, digits=3))

print("(d): ")

print("prob 04_35\n")
print("(a): ")

theta = 35units.angles.degree
f = 60.0units.force.newton
F = var("F")
eqn = f == F*cos( theta.convert(units.angles.radian)/units.angles.radian)
sol = solve( eqn, F )
F = F.subs(sol)/units.force.newton
print("(b): ", F.n(digits=3), " newton")

print("prob 04_41\n")
F_total = var("F_total")

theta = 20units.angles.degree
F_total = 2 200*units.force.newton * cos(theta.convert(units.angles.radian)/units.angles.radian)
F_total = F_total.convert(units.force.newton)
print("the total force is: ", F_total, " i.e., ", n(F_total/units.force.newton, digits=3) )