/******************************************************************************
*									      *
*  INTERPOLATE BETWEEN FRACTIONS					      *
*  05/19/15 (dkc)							      *
*									      *
******************************************************************************/
#include <math.h>
double interp(int N, int D, int N1, int D1, unsigned int I, double init,
	      double *out) {
unsigned int i;
double F,F1,inc;
F=(double)N/(double)D;
F1=(double)N1/(double)D1;
inc=(F1-F)/(double)I;
F=init;
for (i=0; i<I; i++) {
   out[i]=F;
   F=F+inc;
   }
return F1;
}