
      do ie=1,100
         e = float(ie)
         print *,flux(12,e),flux(-12,e),flux(14,e),flux(-14,e)
      enddo
      end

