#include <stdio.h>
#include <omp.h>

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

  //Determine how many threads set by OMP_NUM_THREADS environment variable
  printf("%d threads possible\n",omp_get_max_threads());
  #pragma omp parallel private(i)
  for (i=0;i<4;i++) 
    {
     #pragma omp single
     printf("Currently %d threads in use\n",omp_get_num_threads());
     printf("thread %d working on i=%d\n",omp_get_thread_num(),i);
     a = i;
     b = 2 * i;
     a += b;
    }
 }