//平均値補完
//array:  補完前の配列
//array2: 補完後の配列
//num:    arrayの配列長
//num2:   array2の配列長
void mean_imputation(double *array, double *array2, int num, int num2)
{
  double space, d;
  int i, j, k=0;
  
  if(num < 2 )
    {
      printf("num は 2 以上でければ補完できません\n");
      exit(1);
    }

  for(i=0,j=0; i<num-1; i++,k=0)
    {
      space = (int)((num2-1.0-(double)j)/(num-1.0-(double)i));
      //補完すべき個数を計算
      d = (array[i+1]-array[i])/space;
      while(j<=(i+1)*space)
	{
	  array2[j] = array[i]+d*(double)k;
	  j++;
	  k++;
	}
      j--;
    }
}