1.1.A.sagews

from sympy import *

# SI units:
m  = Symbol("m",  positive=True)
s  = Symbol("s",  positive=True)
kg = Symbol("kg", positive=True)

# Derived units:
N = kg*m/s/s

# Position vector:
(r_x, r_y, r_z)    = (1*m, 0*m, 2*m)

# forces
(F1_x, F1_y, F1_z) = (1*N, 0*N, 1*N)
(F2_x, F2_y, F2_z) = (0*N, S(3)/2*N, 0*N)

# Using matrices:
r = Matrix([r_x,r_y,r_z])

F1 = Matrix([F1_x, F1_y, F1_z])
F2 = Matrix([F2_x, F2_y, F2_z])

# Resultant:
F = F1 + F2
# M = r x F:
M = r.cross(F)

print("\nF, M, |M|:")
pprint(["M / Nm", M/N/m])
pprint(["F / N", F/N])
pprint(["|M| / Nm", (M/N/m).norm()])

F, M, |M|: [M / Nm, [-3 ]] [ ] [ 1 ] [ ] [3/2] [F / N, [ 1 ]] [ ] [3/2] [ ] [ 1 ] [|M| / Nm, 7/2]