program test_rand

    implicit none
    integer :: i, seed(2), M, n
    real( kind = 8 ) :: y

    write(*,*) 'input two integer numbers'
    read(*,*) seed(1), seed(2)
    call random_seed( size = M )
    call random_seed( put = seed(1:M) )
    do i = 1, 10
         call random_number( y )
         n = mod( int(100*y), 6 ) + 1
         write(*,10)  y, n
10     format(2x,'random number = ',d16.6, ' and integer = ', I2)
    end do

    end program