#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>

#define N 10000
double r[N], theta[N], x[N][N];

int main()
{
  int i,j;

  // Initialization
  for (i=0; i<N; i++) r[i] = rand();
  for (i=0; i<N; i++) theta[i] = rand();

  // Start time
  clock_t time0 = clock();

  // Main calculation
    for (j=0; j<N; j++) {
      double d = cos(theta[j]);
      for (i=0; i<N; i++) {
	x[i][j] = r[i] * d;
    }
  }

  // Finish time
  clock_t time1 = clock();

  // Output time
  double time = (double)(time1-time0)/CLOCKS_PER_SEC;
  printf("Time = %15.7f sec\n", time);

  return 0;
}