from sympy.plotting import plot plot(1/x) plot(1/x,(x, 0, 3)) plot(Heaviside(x)) plot(sqrt(x)) plot(-sqrt(sqrt(x)),sqrt(sqrt(x))) plot(LambertW(x)) plot(LambertW(x), sqrt(LambertW(x)), (x, -2, 1))