I have a series of functional dependencies which look something like this.
y = 0.1222*x**0.8628
I would like to use sympy to get out x as a function of y.
Searching the internet I saw that solve can be used for that.
So I tried
x = Symbol('x', real=True, positive=True)
y = Symbol('y', real=True, positive=True)
solve(y - 0.1222*x**0.8628, x)
However, this doesn't work, because after 20 minutes the kernel is still busy. I noted, that if I exchange with 0.8628 to a simpler power, like 0.8 or even 0.86 then the computation works quite quickly.
I would like some help on speeding this up, or maybe there is a much better approach?