program essai

  use mathstuff
  implicit none

  integer                     :: foo
  real                        :: quux
  double precision            :: somme

  call init_random_seed()
  somme = 0.0

  do foo=1, 5
    quux = rand()
    somme = somme + quux
    print *, foo, quux, somme/foo
  enddo

end program