如何使用多个解决方案解决多项式，包括python中的复杂解？

I'm trying to solve the following polynomial in python. Some of the constants (2.11 and .125567481) are dynamic and will change in the program. I need a way to calculate the solutions to this equation in python. I tried using Sympy but it was executing forever and never came up with a solution. Wolfram Alpha gives 4 solutions, 3 of them complex, and 1 real. I only need the real one. Is there any way to solve this equation in python efficiently? It's really easy to do in mathematica, but I can't find a way to do it in python3 that is efficient.I'm trying to solve the following polynomial in