Example solving systems of equations with SymPy in Python for example 3 of plasticity lecture 1

import sympy as sym
sym.init_printing()

w_A, w_C, F = sym.symbols('w_A w_C F')
EA = sym.symbols('EA', positive=True) 

w_B = w_A - (w_A - w_C) / 3 * 2
eq1 = sym.Eq(EA * 4 * 2 * w_B + EA * 2 * w_C*6 - F * 3,0)
eq2 = sym.Eq(EA * 2 * w_A + EA * 2 * w_B + EA * 2 * w_C - F, 0)
sol = sym.solve([eq1, eq2], (w_A, w_C))
display(sol)

w_B.subs(sol)