plot Ei(-inversegammaregularized(0.0001,-x/gamma(0.0001)))