I should add that the recurrence relation for the integral only terminates if α is a positive integer. In the general case you need to evaluate the integral numerically. There are a large number of subroutine libraries for this class of problem, as the algorithm of choice is problem dependent.
See https://fortran-lang.discourse.group/t/numerical-integration-library/ for some clues