module para
  implicit none
  integer,parameter :: imax = 10
  real(8),parameter :: alph = 3.4d0,      &
                     & z0   = 0.070588d0, &
                     & eps  = 1.d15
  integer :: i
end module para
!*********************************************************************
program rogi
  use para
  implicit none
  real(8) :: z1,z2
  
  z1 = z0
  write(*,*) 'z(',i,')','z(',i+1,')'
  do i = 0,imax
    z2 = alph*z1*(1.d0-z1)
    write(*,*) i,z1,i+1,z2
    if(dabs(z2-z1) > eps) then
      write(*,*) 'overflow error!!'
    end if
    z1 = z2
  end do
  
  stop
end program rogi