fndaw(xdaw,ydaw ) = aw0*0.5*( hx2*xdaw**2 + hy2*ydaw**2 )
        fneps(aw) = abs(aw)**2 / ( 2.0d0 *( 1.0d0+abs(aw)**2 ) )
        fnjj(eps) = cjj0 + eps* ( cjj1 + eps*( cjj2 + eps*cjj3 ) )