float tablei(nsample,array,tab)

long nsample;
float *array,*tab;

{
	register loc1,loc2;
        float frac = ((float)(nsample)/(*tab)) * *(tab+1);
	if(frac >= *(tab+1))  {
		loc1 = *(tab+1);
		loc2 = loc1;
		frac = 0;
		}
	else    {
		loc1 = frac;
		loc2 = loc1+1;
		frac = frac - (float)loc1;
		}
	return(*(array+loc1) + frac * (*(array+loc2) - *(array+loc1)));
}