Revision 1245 trunk/src/linad99/combv.cpp
combv.cpp (revision 1245)  

240  240 
static dvariable gammlnguts(const prevariable _z) 
241  241 
{ 
242  242 
double z = value(_z); 
243 
double zdot=1.0; 

244 
const double lpi =1.1447298858494001741434272; 

245 
const double pi =3.1415926535897932384626432; 

243 
//double zdot=1.0;


244 
//const double lpi =1.1447298858494001741434272;


245 
//const double pi =3.1415926535897932384626432;


246  246 
const double lpp =0.9189385332046727417803297; 
247  247 
int n=7; 
248  248 
const double c[9]={0.99999999999980993, 
...  ...  
265  265 
xdot=c[i]*square(zinv); 
266  266 
} 
267  267 
double t=z+n+0.5; 
268 
double tdot=zdot; 

269  268 
//return lpp + (z+0.5)*log(t) t + log(x); 
270  269 
double ans= lpp + (z+0.5)*log(t) t + log(x); 
270 
//double tdot=zdot; 

271  271 
//double ansdot=zdot*log(t) + (z+0.5)/t*tdot tdot +xdot/x; 
272  272 
// since tdot=1.0 
273  273 
// since zdot=1.0 
Also available in: Unified diff