double precision, parameter :: pi=3.14159265359d0
double precision, parameter :: dsqrt_pi=1.77245385091d0