09 Mar '16 14:10>
Originally posted by humyOf course the easiest way of doing this is to make a change of variables to a standard integral. Up to a factor of h your integral is I(m, h + s), where:
Oh damn! You are right. And that was a misprint in my OP and which I propagated throughout my earlier workings causing me (and you ) confusion. Terribly sorry about that. My fault. That should have been;
∫[0, ∞] h (θ^(m-1)) / ( e^( θ*(h + s) )) dθ
I will now systematically iterate through all my files to correct that edit error wherever I have propagated it -that could take me hours.
I(m, a) = ∫[0, ∞] θ^(m-1) exp( -aθ ) dθ
Change variables so that x = aθ, then we have:
I(m, a) = ∫[0, ∞] (x/a)^(m-1) exp(-x) dx/a
= (1/a)^m ∫[0, ∞] dx x^(m-1) exp(-x)
The definition of the gamma function is:
Γ(t) = ∫[0, ∞] dx x^(t-1) exp(-x)
So:
I(m, a) = Γ(m)/a^m ∀a, m
So your integral is hΓ(m)/(h + s)^m. Γ(m) = (m - 1)! for positive integer m (which we can look up on Wikipedia or prove using one of the methods above). This gives your result of h (m - 1)! / (h + s)^m.