implicit real*8(a-h,o-z)
         dimension xp(350000),yp(350000),nbin(10)
c
         pi=4.*datan(1.d0)
         do 1 k=1,343340
         read(24,*) xp(k),yp(k),z,a,b,c,d
      if (xp(k).gt.0..and.yp(k).gt.0.) then
      theta=datan(yp(k)/xp(k))
      else if (yp(k).gt.0.) then
      theta=pi-datan(dabs(yp(k)/xp(k)))
      else if (xp(k).gt.0.) then
      theta=2.*pi-datan(dabs(yp(k)/xp(k)))
      else
      theta=pi+datan(dabs(yp(k)/xp(k)))
      end if
      r=dsqrt(xp(k)**2.+yp(k)**2.)
      if(r.gt.0.2) then
      nth=dint(theta/2.*pi))*10+1
      nbin(nth)=nbin(nth)+1
      end if
1     continue
      do 2 j=1,10
      write(6,*) j,nbin(j)
2     continue
      stop
      end if