double test_write_lsfit71(double a, double b, double c,
       double d, double e, double f, double g) 
{ 
  int i, j; 
  int m = 1; 
  double t, val=0.0; 
  double vpwr[56]; 
  double C[8] = { 
     1.0000000000000009e+00, 
     2.0000000000025109e+00, 3.0000000000047011e+00, 4.0000000000065050e+00, 
     5.0000000000081615e+00, 6.0000000000097398e+00, 7.0000000000113580e+00, 
     8.000000000012985e+00}; 
 
  int pwrs[56] = { 
    0, 0, 0, 0, 0, 0, 0, 
    1, 0, 0, 0, 0, 0, 0, 
    0, 1, 0, 0, 0, 0, 0, 
    0, 0, 1, 0, 0, 0, 0, 
    0, 0, 0, 1, 0, 0, 0, 
    0, 0, 0, 0, 1, 0, 0, 
    0, 0, 0, 0, 0, 1, 0, 
    0, 0, 0, 0, 0, 0, 1}; 
 
  for(j=0; j<7; j++) 
    vpwr[j] = 1.0; 
  for(i=1; i<=1; i++) 
  { 
    vpwr[i*7+0] = vpwr[(i-1)*7+0] * a; 
    vpwr[i*7+1] = vpwr[(i-1)*7+1] * b; 
    vpwr[i*7+2] = vpwr[(i-1)*7+2] * c; 
    vpwr[i*7+3] = vpwr[(i-1)*7+3] * d; 
    vpwr[i*7+4] = vpwr[(i-1)*7+4] * e; 
    vpwr[i*7+5] = vpwr[(i-1)*7+5] * f; 
    vpwr[i*7+6] = vpwr[(i-1)*7+6] * g; 
  } 
  for(i=0; i<8; i++) 
  { 
    t = C[i]; 
    for(j=0; j<7; j++) 
    { 
      t = t * vpwr[pwrs[i*7+j]*7+j]; 
    } 
    val = val + t; 
  } 
  return val; 
}