{
  int i;
  double av=0.0;
  *n=0;
  *m=0;
  for(i=0;i<N;i++)
     av=av+a[i].s;
  av=av/N;           
  for(i=0;i<N;i++)
     if(av<=a[i].s) 
        { 
         b[*n]=a[i];   
         *n=*n+1;
        }    
     else   
       { 
         c[*m]=a[i];   
         *m=*m+1;
        }                
  return av;                       
}