double histograma(int *red,int *percolantes,double *histo,int n)
{
  // *red es la red ya clasificada
  // *percolantes contiene "1" si esa posicion corresponde a una
  // etiqueta percolante.
  // *histo es la distribucion de fragmantos sin normalizar 
  // n es el lado de la red

  int    i,j,*masas;
  double s;

  s = 0.0;

  masas = (int*) malloc(n*n*sizeof(int));

  for(i=0;i<n*n;i++) *(masas+i) = 0;

  for(i=0;i<n*n;i++)
    {
      j = *(red+i);
      if ((*(percolantes+j))==0) (*(masas+j))++;
    }

  for(i=1;i<n*n;i++)
    {
      j = *(masas+i);
      *(histo+j) = (*(histo+j))+1.0;
      if (j) s = s + 1.0;
    }

  *histo = 0;

  free(masas);

  return s;
}


int regresion(double *tr,double *histo,double suma,int n)
{
 // *tr es un puntero de dos valores [R^2, q0 , tau]
 // *histo es la distribución de fragmentos sin normalizar
 // suma es la suma total de los valores de *histo
 // n es el lado de la red

 int    i,m;
 double r,r2,sumX=0,sumX2=0,sumY=0,sumY2=0,sumXY=0,a,b,*x,*y;

 x = (double*) malloc(n*n*sizeof(double));
 y = (double*) malloc(n*n*sizeof(double));

 m = 0;

 sumX = 0.0;
 sumX2= 0.0;
 sumY = 0.0;
 sumY2= 0.0;
 sumXY= 0.0;

 for (i=SMIN;i<SMAX;i++)
   {
     r = (*(histo+i))/suma;

     if (r!=0.0) 
       {
         *(x+m) = log(i);
         *(y+m) = log(r);

         sumX  = sumX + *(x+m);
         sumX2 = sumX2 + (*(x+m))*(*(x+m));
         sumY  = sumY + (*(y+m));
         sumY2 = sumY2 + (*(y+m))*(*(y+m));
         sumXY = sumXY + (*(x+m))*(*(y+m));

         m++;
       }
   }

 // y = a+b.x

 b = (m*sumXY-sumX*sumY)/(m*sumX2-sumX*sumX); // pendiente
 a = (sumY - b*sumX)/m;                       // ordenada al origen

 r2 = (m*sumXY-sumX*sumY)*(m*sumXY-sumX*sumY);
 r2 = r2/((m*sumX2-sumX*sumX)*(m*sumY2-sumY*sumY));

 *(tr+0) = r2;
 *(tr+1) =  a;
 *(tr+2) = -b;

 free(x);
 free(y);

 return 1;
}