program demo_lamsf use M_datapac, only : lamsf implicit none ! call lamsf(x,y) end program demo_lamsf