#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <unistd.h>
#include <sys/types.h>

// #include "rand.h"

#define MAX_FREC 50
#define MAX_ITER 1000

int vector[MAX_FREC];

void inicializar_vector (void)
{
	int i;
	for (i=0;i<MAX_FREC;i++) {
		vector[i] = 0;
	}
}

int main (void) 
{
	int i,j;
	double aleat,z,v;

	inicializar_vector();
	
	srand(getpid());
	
	for (i=0;i<MAX_ITER;i++) {
		aleat = (rand()%1000) / 999.0;
		// aleat = rand(i);
		z = (pow(aleat,0.135)-(pow((1-aleat),0.135)))/0.1975;

		if ((z >= -4) || (z <= 4)) {
			v = 8.0 / MAX_FREC; // intervalo [-4,4]  -->  longitud 8
			vector[abs(z/v)+1]++;
			// printf("índice: %d\n",abs(z/v)+1);
			// printf("%f\n",v);
		}
		
		printf("i: %d  aleat: %f  z: %f\n",i,aleat,z);
		
	}
	for (j=0;j<MAX_FREC;j++) {
		printf("indice: %d  valor: %d\n",j,vector[j]);
	}

	return 0;
}