double complex fndaw
        double precision xdaw, ydaw
        double precision fneps
        double complex aw
        double precision fnjj, eps
        double precision cjj0, cjj1, cjj2, cjj3
      parameter( cjj0 = 9.9999023630596d-1,
     %           cjj1 = -4.995632267751d-1,
     %           cjj2 = -2.541359939218d-1,
     %           cjj3 = 7.6140808316962d-2 )