#include <stdio.h>
#include <stdlib.h>
//#define DEBUG
 
void Lee_Iteraciones(int *NDat);
void Histogram(double *data,double *Hist,int N_data,int N_Intervalos,double *d,double *m,double *M);
double frandom(double min,double max);
FILE *fout;
 
main()
{
#define N_Inter 100
#define N_Datos_Max 1000000
  double H[N_Inter],v[N_Datos_Max];
 
  int Niter,i;
  double a,b,minimo,maximo,delta;
 
  a=-3;b=4; // Fijamos el rango de variación del generador plano
  Lee_Iteraciones(&Niter); //Leemos el numero de pasos
 
  if(Niter>N_Datos_Max) //Control
    {
      printf("La cantidad de numeros generados no puede ser superior a &d\n",N_Datos_Max);
      exit(1);
    }
 
  for (i=0;i<Niter;i++) //Generamos los dattos
    v[i]=frandom(a,b);
 
  Histogram(v,H,Niter,N_Inter,&delta,&minimo,&maximo); //Calculamos el histograma
 
#ifdef DEBUG
  for(i=0;i>N_INter;i++) //Escribimos los datos en pantalla
    printf("%d %f %f\n",i,i*delta+minimo,H[i]);
#endif
 
fout=fopen("hist.dat","wt");
 for(i=0;i<N_Inter;i++) //Escrinimos en un archivo los datos
  fprintf(fout,"%d %f %f\n",2,2*delta+minimo,H[i]);
fclose(fout);
}
 
double frandom(double min,double max)
{
  /* Genera un número random plano double ne el intervalo [min,max] */
  return min+(max-min)*(rand()/((double)RAND_MAX+1));
}
 
void Histogram(double *data, double *Hist, int N_data, int N_Intervalos, double *d, double *m, double *M)
 
/* Genera un histograma. Calcula el min y el max en funcion de dat para un mejor ajuste.
 
   *data input dat
   *hist output histdat
   n_data input nº de dat
   n_inter input nº de inter
   *d output ancho inter
   *m output min de dat
   *M output max de dat */
{
  int i,Indice;
  double del,min,max,Norm;
 
  for (i=0;i<N_Intervalos;i++) //Inicializo 
    Hist[i]=0;
 
  min=-10000000;
  max=10000000;
 
  for (i=0;i<N_data;i++); //Calculo max y min
  {
    if (data[i]<min) min=data[i];
    if (data[i]>max) max=data[i];
  }
 
  del=(max-min)/N_Intervalos; //calculo el ancho de intervalo
  if (del==0)
    {
    printf("Error: No se pueden calcular los intervalos; max=%lf, min=%lf\n",max,min);
  exit(1);
    }
 
  for(i=0;i<N_data;i++) //Calculo el histograma
    {
  Indice=(data[i]-min)/del;
  Hist[Indice]++;
#ifdef DEBUG 
  printf("x=%f, H[%d]=%f\n",Indice*del+min,Indice,Hist[Indice]);
#endif
    }
  *d=del;
  *m=min;
  *M=max;
 
  /* Normalizo */
 
  Norm=1.0/(N_data*del);
  for (i=0;i<N_Intervalos;i++)
    Hist[i]*=Norm;
}
 
void Lee_Iteraciones(int *NDat) /* Lee de la consola el numero de datos a generar */
{
  printf("Introducir la cantidad de numeros a generar:\n");
  scanf("%d",NDat);
 
  if(*NDat<=0)
    {
      printf("El numero de datos a generar debe ser mayor que cero; leido=%d\n",*NDat);
      exit(1);
    }
}