program main
      implicit none
      real*8 flux,energy
      real*8 params(2),test
      real*8 energies(16),expflux(16),error(16)
      integer i,j,k,flag,counter

      open (10, file='excess.dat',status='unknown')
      do i=1,16
      read(10,*)energies(i),expflux(i),error(i)
      end do
      close(10)
c      print*,'experimental data file has been read'

      do i=1,50
          params(1)=90.+5.*i

       do j=1,50
          params(2)=10.**(-6+0.04*j)

         counter=0
         do k=1,16
            energy=energies(k)
            call yieldsub(energy,params,flux)
            test=expflux(k)-flux
                if (test .ge. 0.) then
                   counter=counter+1
                else
                   counter=counter
                end if

         end do

             if (counter .lt. 16) then
                flag = 0
             else
                flag = 1
             end if

c            print*,k,energies(k),params(1),params(2),flag
            write(*,'(f10.2,e15.5,2i5)')params(1),params(2),
     .                               counter,flag



       end do
      end do

      end program