phi = solve(eqn,x)[0].rhs() # x en fonction de u : fonction phi(u)=x phi_inv = solve(eqn,u)[0].rhs() # u en fonction de x : inverse de phi : phi_inv(x)=u ...