Remove `gamfun` and replace it with the analytical solution from scipy

[chris-ashe]

Currently within the profiles module the gamma function is calculated by gamfun within the maths_library.f90. This is a Stirling approximation approach. SciPy is able to calculate the full integral with higher accuracy