from sympy import *
from sympy import summation
import numpy as np
a, b, c = symbols('a b c')
a=2*a+1
print(a)
from sympy import *
left, right, pln, a = symbols('left right pln a')
left = a
right = -2 * a + 1
pln = left - right
print(pln) # 3*a - 1
print(pln / polys.polytools.LC(pln)) # a - 1/3
print(solvers.solve(pln, a)) # [1/3]