| ... | ... |
@@ -797,3 +797,3 @@ |
| 797 | 797 |
} while(nu>std::pow(xi,delta-1.0)*std::exp(-xi)); |
| 798 |
return theta*(xi |
|
| 798 |
return theta*(xi+gamma(int(std::floor(k)))); |
|
| 799 | 799 |
} |
| ... | ... |
@@ -821,3 +821,3 @@ |
| 821 | 821 |
{
|
| 822 |
return exponential(gamma(k,1.0/x_min)); |
|
| 822 |
return exponential(gamma(k,1.0/x_min))+x_min; |
|
| 823 | 823 |
} |
0 comments (0 inline)