Skip to content

Instantly share code, notes, and snippets.

@certik

certik/fm.f90

Created September 22, 2012 19:56
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save certik/3767621 to your computer and use it in GitHub Desktop.
Save certik/3767621 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
fm.f90
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