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

#define N 100000

int main(){

  int i;
  double A[N+2];
  double B[N+2];
  double suma;

  for(i=0;i<N+2;i++) A[i] = (double)i/N;
  for(i=0;i<N+2;i++) B[i] = 1.0 - (double)i/N;

  double t1 = omp_get_wtime();
  for(i=0; i<N; i++){
    A[i] += A[i+2] + sin(B[i]);
  }
  t1 = omp_get_wtime() - t1;

  suma = 0.0;
  for(i=0;i<N+2;i++) suma+=A[i];
  printf("suma %lf, czas obliczen %lf\n", suma, t1);



  for(i=0;i<N+2;i++) A[i] = (double)i/N;
  for(i=0;i<N+2;i++) B[i] = 1.0 - (double)i/N;

  t1 = omp_get_wtime();

  // wersja równoległa
	#pragma omp parallel default(none) shared(A, B, i)
{
int i;
	for(i=0; i<N; i++){
		#pragma omp critical
		{
		A[i] += A[i+2] + sin(B[i]);
		}
}
}

  t1 = omp_get_wtime() - t1;

  suma = 0.0;
  for(i=0;i<N+2;i++) suma+=A[i];
  printf("suma %lf, czas obliczen rownoleglych %lf\n", suma, t1);


}