program demo_expsf use M_datapac, only : expsf implicit none ! call expsf(x,y) end program demo_expsf