real function evtsnue(e,pol,emu)

c
c gives mumu 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*12./4./3.14159265 
      val = val*((1-x)+(1-x)*pol)

c 0.34e-38*1e-4*e*6.e32*1.e20/2900./2900.

      write(6,*) val


      evtsnue = val*e*2425.

      return 
      end