real function evtsnumu(e,pol,emu)

c
c gives numu events  from muon decay at 2540 km 
c in 1/GeV/muon  in 1 kt detector for 1.e20 muons
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)

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

      evtsnumu = val*e*4780.

      return 
      end