program test
      implicit none
      integer :: a(250,250),center(2),radius
      real(8) :: rd,distance
      integer :: i,j,k
        
      call random_seed()
      a(:,:) = 0
      do i = 1,100
          print*, '第',i,'个圆'
          call random_number(rd)
          center(1) = int(250*rd) ! [0,250]
          call random_number(rd)
          center(2) = int(250*rd) ! [0,250]
          call random_number(rd)
          radius = int(5*rd)+5    ! [5,10]
          print*, '圆心(随机): ', '(',center,')'
          print*, '半径(随机): ', radius
           
          print*, '半径内的点(值为1) ... 略'
          do j = 1,250
              do k = 1,250
                  distance = sqrt( (j-center(1)*1.0)**2 + (k-center(2)*1.0)**2 )
                  if(int(distance) <= radius) then
                      a(j,k) = 1
                  end if

                  if(a(j,k) == 1) then
                      !print*, '(',j,k,')'
                  end if
              end do
          end do
      end do
      
      stop
      end program