do i = 1, n
    do j = 1, n

      x = ( ( n - j ) * xmin + ( j - 1 ) * xmax ) / ( n - 1 )
      y = ( ( n - i ) * ymin + ( i - 1 ) * ymax ) / ( n - 1 )
      pixel(i,j) = 0

      x1 = x
      y1 = y

      do k = 1, 1000

        x2 = x1 * x1 - y1 * y1 + x
        y2 = 2 * x1 * y1 + y

        if ( x2 < -2.0 .or. 2.0 < x2 .or. y2 < -2.0 .or. 2.0 < y2 ) then
          pixel(i,j) = 1
          exit
        end if

        x1 = x2
        y1 = y2

      end do
    end do
  end do