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

#define MAX 3000

void histograma(double *y,double *x,int n,double a,double b,int M);

int main()
{

  int i,m;
  double a,b;
  double x[MAX],y[2*MAX];

  m=10;
  a=1.0;
  b=6.0;

  x[0]=3.35;
  x[1]=7.0;
  x[2]=4.75;
  x[3]=4.75;
  x[5]=1.1;

  histograma(y,x,5,a,b,m);

  for(i=0;i<m;i++)
    {
      printf("%f\t%f\n",y[i],y[i+m]);
    }

  return 1;

}


void histograma(double *y,double *x,int n,double a,double b,int m)
{

  /* 
     Histograma del vector de datos x[0]...x[n-1] 
     
     M  = cantidad de columnas del histograma
     a,b= limites del del histograma 
     
     y[0]...y[m-1]  valores de las columnas (normalizados)
     y[m]...y[2m-1] valores de las marcas de clase

  */

  int    i,j;
  double h,hh,s;

  s=1.0/(double)n;
  h=(b-a)/(double)m;
  hh=h/2.0;

  for(i=0;i<m;i++)
    {
      *(y+i)=0.0;
      *(y+m+i)=(double)i*h+a+hh;
    }

  for(i=0;i<n;i++)
    {
      j=(int)floor((*(x+i)-a)/h);

      if (j<0) j=0;
      if (j>m) j=m-1;

      y[j]=y[j]+s;
    }
}