program main

      include 'omp_lib.h'

      integer n
      parameter ( n = 10000000 )

      integer i
      double precision x(n), y(n), s

      s = 123.456

      do i = 1, n
        x(i) = mod ( i, 17 )
        y(i) = mod ( i, 31 )
      end do

c$omp parallel

c$omp do
      do i = 1, n
        y(i) = y(i) + s * x(i)
      end do
c$omp end do

c$omp end parallel

      stop
      end