double test_write_lsfit72(double a, double b, double c,
       double d, double e, double f, double g) 
{ 
  int i, j; 
  int m = 2; 
  double t, val=0.0; 
  double vpwr[252]; 
  double C[36] = { 
     1.0000000030592109e+00, 
     2.0000000000028639e+00, 3.0000000000048535e+00, 4.0000000000064890e+00, 
     5.0000000000081322e+00, 6.0000000000097700e+00, 7.0000000000113864e+00, 
     8.0000000000131504e+00, 8.9999999992626751e+00, 9.9999999999859792e+00, 
     1.5999999999099918e+01, 1.0999999999986276e+01, 1.6999999999963961e+01, 
     2.1999999999087770e+01, 1.1999999999984873e+01, 1.7999999999962245e+01, 
     2.2999999999950898e+01, 2.6999999999081854e+01, 1.2999999999983697e+01, 
     1.8999999999959751e+01, 2.2999999999952468e+01, 2.7999999999942290e+01, 
     3.0999999999070432e+01, 1.3999999999982474e+01, 1.9999999999958106e+01, 
     2.4999999999948088e+01, 2.8999999999940304e+01, 3.1999999999934371e+01, 
     3.3999999999064606e+01, 1.4999999999981210e+01, 2.0999999999956554e+01, 
     2.5999999999946443e+01, 2.9999999999937760e+01, 3.2999999999932065e+01, 
     3.4999999999927226e+01, 3.599999999905840e+01}; 
 
  int pwrs[252] = { 
    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, 
    2, 0, 0, 0, 0, 0, 0, 
    1, 1, 0, 0, 0, 0, 0, 
    0, 2, 0, 0, 0, 0, 0, 
    1, 0, 1, 0, 0, 0, 0, 
    0, 1, 1, 0, 0, 0, 0, 
    0, 0, 2, 0, 0, 0, 0, 
    1, 0, 0, 1, 0, 0, 0, 
    0, 1, 0, 1, 0, 0, 0, 
    0, 0, 1, 1, 0, 0, 0, 
    0, 0, 0, 2, 0, 0, 0, 
    1, 0, 0, 0, 1, 0, 0, 
    0, 1, 0, 0, 1, 0, 0, 
    0, 0, 1, 0, 1, 0, 0, 
    0, 0, 0, 1, 1, 0, 0, 
    0, 0, 0, 0, 2, 0, 0, 
    1, 0, 0, 0, 0, 1, 0, 
    0, 1, 0, 0, 0, 1, 0, 
    0, 0, 1, 0, 0, 1, 0, 
    0, 0, 0, 1, 0, 1, 0, 
    0, 0, 0, 0, 1, 1, 0, 
    0, 0, 0, 0, 0, 2, 0, 
    1, 0, 0, 0, 0, 0, 1, 
    0, 1, 0, 0, 0, 0, 1, 
    0, 0, 1, 0, 0, 0, 1, 
    0, 0, 0, 1, 0, 0, 1, 
    0, 0, 0, 0, 1, 0, 1, 
    0, 0, 0, 0, 0, 1, 1, 
    0, 0, 0, 0, 0, 0, 2}; 
 
  for(j=0; j<7; j++) 
    vpwr[j] = 1.0; 
  for(i=1; i<=2; 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<36; i++) 
  { 
    t = C[i]; 
    for(j=0; j<7; j++) 
    { 
      t = t * vpwr[pwrs[i*7+j]*7+j]; 
    } 
    val = val + t; 
  } 
  return val; 
}