Revision 1110 trunk/src/df1b2-separable/df1b2f20.cpp

df1b2f20.cpp (revision 1110)
26 26
  const double pi =3.1415926535897932384626432;
27 27
  const double lpp =0.9189385332046727417803297;
28 28
  int n=7;
29
  const double c[9]={0.99999999999980993, 
30
    676.5203681218851, 
29
  const double c[9]={0.99999999999980993,
30
    676.5203681218851,
31 31
    -1259.1392167224028,
32
     771.32342877765313, 
33
    -176.61502916214059, 
32
     771.32342877765313,
33
    -176.61502916214059,
34 34
    12.507343278686905,
35
     -0.13857109526572012, 
36
    9.9843695780195716e-6, 
35
     -0.13857109526572012,
36
    9.9843695780195716e-6,
37 37
    1.5056327351493116e-7};
38 38
  z-=1.0;
39 39
  double x=c[0];
......
49 49
    xdot-=c[i]*square(zinv);
50 50
    x2dot+=2.0*c[i]*cube(zinv);
51 51
    x3dot-=6.0*c[i]*fourth(zinv);
52
  }    
52
  }
53 53
  double t=z+n+0.5;
54 54
  double tdot=zdot;
55 55
  //return lpp + (z+0.5)*log(t) -t + log(x);

Also available in: Unified diff