سلام دکترجان. میتونین از روش مونت کارلو برای محاسبه این انتگرال استفاده کنید.

سلام. تشکر فراوان از راهنماییتون. برای محاسبه تابع گشتاور یک توزیع در یه مقاله لازمش دارم.

Hello

I think when b = a+1, you can easily integrate it. Or you can use Monte Carlo approximation.

Thanks

Hello. Thank you very much for your reply

let h(u,b)=u^a*(exp(gamma(b,u)).

Then the integral is the expected value of h(U) with U distributed as standard exponential r.v.

Simulate a large number, say, 100000 standard exponential random variables U(I), i=1,...,100000. Then a very close approximation to the integral is mean of {h(Ui), i=1,...,100000}.

Thank you very much for your reply