3.3.N.sagews

from sympy import *

oneoverthree = S(1)/3
eightoverthree = S(8)/3
threeovertwo = S(3)/2

print("\n--- User input -----------------------")

# mass m and length l in Mahnken-task:
mM, lM = var("mM, lM")

ome = var("omega")

# mass of bar:
mS = threeovertwo * mM
xS = oneoverthree * lM

ThetaS = eightoverthree * mM * lM**2

print("\n--- a --------------------------------")

ThetaA = ThetaS + mS*xS**2
pprint(ThetaA)

print("\n--- b --------------------------------")

vS = omega*xS
pprint(vS)

print("\n--- c --------------------------------")
omega_b = var("omega_b")
b, F = var("b, F")
dom = omega_b - omega

eq1 = Eq( ThetaA*dom , -b*F )

print("\n--- d --------------------------------")

eq1 = Eq( ThetaS*dom , -(b-xS)*F )

print("\n--- e --------------------------------")

eq2 = Eq( mS*xS*dom , -F )

print("\n--- f --------------------------------")

sol=solve([eq1, eq2], [b, F])
sol=sol[0]
b=sol[0]

print("\n--- g --------------------------------")

e = var("e")
F=sol[1]
F = F.subs(omega_b, -e*omega)
pprint(F)

--- User input ----------------------- --- a -------------------------------- 2 17*lM *mM --------- 6 --- b -------------------------------- lM*omega -------- 3 --- c -------------------------------- --- d -------------------------------- --- e -------------------------------- --- f -------------------------------- --- g -------------------------------- lM*mM*(e*omega + omega) ----------------------- 2