program binwrite
  implicit none
  integer, parameter :: n = 20
  integer            :: i
  real               :: x(n)

  do i=1,n 
     call random_number(x(i))
  enddo

  ! ascii output
  open(unit=1,file='myascii.dat')
  write (1,*) n
  do i=1,n
     write (1,*) i, x(i)
     write (*,*) i, x(i)
  enddo
  close(1)

  ! binary output
  open(unit=1,file='mybin.dat',form='unformatted')
  write (1) n
  write (1) x
  close(1)


end