C
      IMPLICIT NONE
      REAL*8 delta_x, prob(1:100), x
      INTEGER ibin, n, nbin, nmax
      DATA prob/100*0.0D0/
      nbin = 100                            ! ���̐�
      delta_x = 1.0D0/nbin
      n = 0                                 ! ������
      x = 0.6D0
      nmax = 10000
      DO WHILE(n .LT. nmax)
         x = 4*x*(1-x)
         n = n + 1
         ibin = int(x/delta_x) + 1
         prob(ibin) = prob(ibin) + 1
      END DO
      DO ibin = 1, nbin
         IF(prob(ibin) .GT. 0) 
     +      WRITE(*,'(I12,F13.6)') ibin, prob(ibin)/nmax
      END DO
      END