Skip to content

Instantly share code, notes, and snippets.



Created Sep 22, 2012
What would you like to do?
subroutine Fm(maxm, t, F)
integer, intent(in) :: maxm
real(dp), intent(in) :: t
real(dp), intent(out) :: F(0:)
integer :: m
do m = 0, maxm
F(m) = dgami(m+0.5_dp, t) / (2*t**(m+0.5_dp))
end do
end subroutine
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment