demo_causf Program

     program demo_causf
     use M_datapac, only : causf
     implicit none
     ! call causf(x,y)
     end program demo_causf