real function flxnumu(e,pol,emu)

c
c gives numu flux from muon decay at 1 km 
c in 1/GeV/m^2/muon
c
c

      real e, pol, emu
      real val
      real x 

c      pol = 1. 
c      emu = 4. 

      gam = emu/0.10566
      x = e/emu 
      
      val = (1./emu/1.e6)
      val = val*gam*gam*x*x*2./4./3.14159265 
      val = val*((3-2*x)+(1-2*x)*pol)

      flxnumu = val

      return 
      end