program demo_gamcdf use M_datapac, only : gamcdf implicit none ! call gamcdf(x,y) end program demo_gamcdf