/******************************************************************************
*									      *
*  INTERPOLATE BETWEEN FRACTIONS					      *
*  04/02/10 (dkc)							      *
*									      *
******************************************************************************/
#include <math.h>
double interp1(double ND, double ND1, unsigned int I, double init,
	      double *out) {
unsigned int i;
double F,inc;
inc=(ND1-ND)/(double)I;
F=init;
for (i=0; i<I; i++) {
   out[i]=F;
   F=F+inc;
   }
return ND1;
}