c  カオス：ロギスティック方程式
      real t,x,dt,t0,x0
c  i：世代数、x：個体数
c  x0：初期値
c  nlim：時間の最大値
c  a：パラメータ

      open(unit=2,file='logistic.dat')
      write(6,*) ' a ?'
      read(5,*) a
      write(6,*) ' x0 ?'
      read(5,*) x0
      write(6,*) ' nlim ?'
      read(5,*) nlim
      nwrite=nlim/10
      if(nwrite.eq.0) nwrite=1
      i=0
      x=x0                        
      write(2,601) i,x
   10 continue
      i=i+1
      x=a*(1.0-x)*x
      write(2,601) i,x
      if(mod(i,nwrite).eq.0) then
        write(6,601) i,x
      endif
      if(i.gt.nlim) then
        stop
      endif
      go to 10                                                          
  601 format(1h ,i5,2x,1pe12.5)                                         
      end