implicit none
  integer               :: j, n
  real,dimension(10000) :: x
  real                  :: t1,t2
  . . . .
  call CPU_TIME( t1 )
  !do j=1,10000
  !  x(j)=log(23.)*j
  !end do
  x(:) = log(23.)*[ (j,j=1,10000) ]
  call CPU_TIME( t2 )
  . . . .
  read *, n ; print *, x(n)
  print *,t2 - t1
  . . . .