program test_gamma
use slatec_module
print*,"4! = ", exp(ln_gamma (5.0))
end program test_gamma
