# Calculate a definite integral containing an ExpIntegral

I need to calculate the following integral

Integrate[Exp[-(y/bn)]  (-y)^c ExpIntegralE[d, p y], {y, 0, z/x}]


Where $$z<0,z\in \mathbb{R}$$, $$x\in\mathbb{R},x>0$$,$$c>0,c\in\mathbb{Z}$$,$$d\leq0,d\in\mathbb{Z}$$, $$p\in\mathbb{R},p<0$$. I tried to replace the Exponential Integral with functions from http://functions.wolfram.com/GammaBetaErf/ExpIntegralE/introductions/ExpIntegrals/ ,but mathematica still can't calculate the integral. Any ideas of how or where to look for a workaround?

• Probably you should try numerical integration... You are trying to integrate a value that contains an integration itself. I think this is too much for Mathematica – user59583 Apr 23 '19 at 10:02