32 |
{ |
{ |
33 |
if (x < 0) |
if (x < 0) |
34 |
{ |
{ |
35 |
return 0 ; |
return 0; |
36 |
} |
} |
37 |
else |
else |
38 |
{ |
{ |
39 |
double P = - expm1(-x/mu); |
double P = -expm1 (-x / mu); |
40 |
|
|
41 |
return P; |
return P; |
42 |
} |
} |
43 |
} |
} |
47 |
{ |
{ |
48 |
if (x < 0) |
if (x < 0) |
49 |
{ |
{ |
50 |
return 1 ; |
return 1; |
51 |
} |
} |
52 |
else |
else |
53 |
{ |
{ |
54 |
double Q = exp (-x/mu); |
double Q = exp (-x / mu); |
55 |
|
|
56 |
return Q; |
return Q; |
57 |
} |
} |
58 |
} |
} |
60 |
double |
double |
61 |
gsl_cdf_exponential_P_inv (const double P, const double mu) |
gsl_cdf_exponential_P_inv (const double P, const double mu) |
62 |
{ |
{ |
63 |
double x = -mu * log1p(-P); |
double x = -mu * log1p (-P); |
64 |
|
|
65 |
return x; |
return x; |
66 |
} |
} |
67 |
|
|
68 |
double |
double |
69 |
gsl_cdf_exponential_Q_inv (const double Q, const double mu) |
gsl_cdf_exponential_Q_inv (const double Q, const double mu) |
70 |
{ |
{ |
71 |
double x = -mu * log(Q); |
double x = -mu * log (Q); |
72 |
|
|
73 |
return x; |
return x; |
74 |
} |
} |