void copy_array2(double a[], int length, double b[])
{
  int i;
  for ( i = 0; i < length; i++ )
    *b++ = *a++;
}