The module does not cope well with discontinuities.
plot(1/x)
plot(1/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))], (-2,1))