equal
deleted
inserted
replaced
764 /// |
764 /// |
765 double gamma(double k,double theta=1.0) |
765 double gamma(double k,double theta=1.0) |
766 { |
766 { |
767 double xi,nu; |
767 double xi,nu; |
768 const double delta = k-std::floor(k); |
768 const double delta = k-std::floor(k); |
769 const double v0=M_E/(M_E-delta); |
769 const double v0=E/(E-delta); |
770 do { |
770 do { |
771 double V0=1.0-real<double>(); |
771 double V0=1.0-real<double>(); |
772 double V1=1.0-real<double>(); |
772 double V1=1.0-real<double>(); |
773 double V2=1.0-real<double>(); |
773 double V2=1.0-real<double>(); |
774 if(V2<=v0) |
774 if(V2<=v0) |