# include <cstdlib>
# include <iostream>
# include <omp.h>
using namespace std;

int main ( int argc, char *argv[] )
{
  int i, n = 1000;
  double total, x[1000];

  for ( i = 0; i < n; i++ )
  {
    x[i] = ( double ) rand ( ) / ( double ) RAND_MAX;
  }

  total = 0.0;
# pragma omp parallel for
  for ( i = 0; i < n; i++ )
  {
    total = total + x[i];
  }
  cout << "Sum = " << total << "\n";
  return 0;
}