integer function CPDSA(n,z,cdn)
IMPLICIT DOUBLE PRECISION (A-B,D-H,O-Y)
EPS=1.0D-15
PI=3.141592653589793D0
CA0=EXP(-.25D0*Z*Z)
VA0=0.5D0*(1.0D0-N)
end function