Hola, buenas a todos!

Tengo un problemilla con este programa.

Tengo que hacer que de como resultado el numero pi pero no se que hago mal que las interaciones no me funcionan correctamente (algunos numeros se repiten) y hace que el proceso falle. Esta hecho con threads. Cualquier ayuda sera bienvenida!

# include <stdio.h>
# include <stdlib.h>
# include <string.h>
# include <unistd.h>
# include <pthread.h>
double factores [1000]; 
int i,N;
double step,x, totalsum, pi;
 
void *calcularsum (void *argument){
	x = (i+0.5)*step;
	factores[i] = 4.0/(1.0+x*x);
}
 
int main (int argc, char *argv[]) {
	if (argc != 2){
		printf("Se ha de introducir un valor. Ejemplo: ./pi 10\n");
}else{
	N=atoi(argv[1]);
	step = 1.0/N;
	pthread_t h[N];
 
	for (i=0; i<N; i++) { 
		pthread_create(&h[i] , NULL , calcularsum , NULL);
 
	}
	for (i=0; i<N; i++){
		pthread_join(h[i], NULL);
		totalsum = factores[i] + totalsum;
		printf("%f\n",totalsum);
	}
	pi = totalsum * step;
	printf("%.10f\n", pi);
}
	return 0;
}
 

Al hacer el "printf("%f\n",totalsum); " para comprobar que numeros van haciendo los diferentes threads, me encuentro con que alguna vez (random) la hace bien, pero lo mas normal es que el primer número salga como 0.000000 y/o algunos otros se repitan.

El resultado de ejecucion seria:
./pi 10
0.000000
3.911980
7.676686
7.676686
11.003090
14.074107
16.886058
19.446058
21.768264
23.870761
2.3870760525

En cambio otra vez, es un número mas aproximado a pi:
./pi 10
3.990025
7.902005
11.666711
15.230186
18.556589
21.627606
24.439557
26.999557
29.321763
31.424260
3.1424259850