#include <stdio.h>
#include <omp.h>
#include <malloc.h>
#define N 10

int factorial(int x)
 {
  if (x == 0) return(1);
  return(x * factorial(x-1));
 }

void main(int argc, char **argv )
 {
  int i;
  int *a,*b,*c;

  a = malloc(N * sizeof(int));
  b = malloc(N * sizeof(int));
  c = malloc(N * sizeof(int));
 
  for(i=1;i<=N;i++) a[i-1] = factorial(i);

  #pragma omp parallel private(i)
   {
    printf("thread %d working on array b\n",omp_get_thread_num());
    for (i=0;i<N;i++) b[i] = a[i] * 2;
    for (i=0;i<N;i++) b[i] = b[i] + 1;
   }

  #pragma omp parallel private(i)
   {
    printf("thread %d working on array c\n",omp_get_thread_num());
    for (i=0;i<N;i++) c[i] = a[i] * 3;
    for (i=0;i<N;i++) c[i] = c[i] + 2;
   }

   for(i=0;i<N;i++)
    printf("a[%d]=%d b[%d]=%d c[%d]=%d\n",i,a[i],i,b[i],i,c[i]);

 }