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