Как решить ошибку при расчете?

При выполнении кода:

p = sym.Symbol('p')#z3.BitVec('p', 128)
q = sym.Symbol('q')#z3.BitVec('q', 128)
x = sym.Symbol('x')
y = sym.Symbol('y')

p1 = 1.732 # sym.Symbol('p1')#z3.BitVec('p', 128)
q1 = 1 # sym.Symbol('q1')#z3.BitVec('q', 128)
p2 = 1.732 # sym.Symbol('p1')#z3.BitVec('p', 128)
q2 = 1 # sym.Symbol('q1')#z3.BitVec('q', 128)
x1 = 0
y1 = 0
x2 = 2
y2 = 0
eq = sym.solve({calc.equation12_minus(x, y, p, q, f, x1, y1),
                calc.equation14_minus(x, y, p, q, f, x1, y1, p1, q1),
                calc.equation12_plus(x, y, p, q, f, x2, y2),
                calc.equation14_plus(x, y, p, q, f, x2, y2, p2, q2)},
                {p, q, x, y}, -1)

Появляется ошибка:

File "/*/CalculationOfQuarrySides/venv/lib/python3.7/site-packages/sympy/functions/elementary/complexes.py", line 58, in eval
    elif arg.is_imaginary or (S.ImaginaryUnit*arg).is_extended_real:
TypeError: unsupported operand type(s) for *: 'ImaginaryUnit' and 'FiniteSet'

Как можно это исправить?

В маткаде все работает.