1.2.I.sagews

from sympy.physics.units import *
from sympy import *

# Given symbols:
F = 4*newton
a, b, c = 4*m, 5*m, 6*m

tex = True
tex = False

pprint("\n--- a ------------------------------------")
pprint("\nrP/m, rQ/m:")
rP = Matrix([0, a, c])
rQ = Matrix([a, b, 0])

if tex:
    pprint(latex(rP/m))
    pprint(latex(rQ/m))
else:
    pprint(rP/m)
    pprint(rQ/m)

pprint("\n--- b ------------------------------------")
pprint("\ndPO/m, dPQ/m:")
dPO = -rP
dPQ = rQ - rP
if tex:
    pprint(latex(dPO/m))
    pprint(latex(dPQ/m))
else:
    pprint(dPO/m)
    pprint(dPQ/m)

pprint("\ne, f:")
e = dPO/dPO.norm()
f = dPQ/dPQ.norm()
if tex:
    pprint(latex(e))
    pprint(latex(N(e,2)))
    pprint(latex(f))
    pprint(latex(N(f,2)))
else:
    pprint(N(e,2))
    pprint(N(f,2))

pprint("\n--- c ------------------------------------")
pprint("\nF.f/N:")
F = F*f
if tex:
    pprint(latex(F/newton))
    pprint(latex(N(F/newton,2)))
else:
    pprint(F/newton)
    pprint(N(F/newton,2))

pprint("\n--- d ------------------------------------")
pprint("\nFe/N:")
Fe = F.dot(e)*e
if tex:
    pprint(latex(N(Fe/newton,2)))
else:
    pprint(N(Fe/newton,2))

pprint("\n--- e ------------------------------------")
pprint("\n||Fe||/N:")
Fe = Fe.norm()
if tex:
    pprint(latex(N(Fe/newton,2)))
else:
    pprint(N(Fe/newton,2))

--- a ------------------------------------ rP/m, rQ/m: [0] [ ] [4] [ ] [6] [4] [ ] [5] [ ] [0] --- b ------------------------------------ dPO/m, dPQ/m: [0 ] [ ] [-4] [ ] [-6] [4 ] [ ] [1 ] [ ] [-6] e, f: [ 0 ] [ ] [-0.55] [ ] [-0.83] [0.55 ] [ ] [0.14 ] [ ] [-0.82] --- c ------------------------------------ F.f/N: [ ____ ] [ 16*\/ 53 ] [ --------- ] [ 53 ] [ ] [ ____ ] [ 4*\/ 53 ] [ -------- ] [ 53 ] [ ] [ ____ ] [-24*\/ 53 ] [-----------] [ 53 ] [2.2 ] [ ] [0.55] [ ] [-3.3] --- d ------------------------------------ Fe/N: [ 0 ] [ ] [-1.4] [ ] [-2.0] --- e ------------------------------------ ||Fe||/N: 2.4