Numerical Integration in fortran with infinity as one of the limits

Question

I am asked to normalize a probability distribution P=A(x^2)(e^-x) within 0 to infinity by finding the value for A. I know the algorithms to ca