/[gsl]/gsl/cdf/exponential.c
ViewVC logotype

Diff of /gsl/cdf/exponential.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1.2.1 by bjg, Mon Feb 17 17:28:32 2003 UTC revision 1.1.2.2 by bjg, Mon Feb 17 17:29:28 2003 UTC
# Line 32  gsl_cdf_exponential_P (const double x, c Line 32  gsl_cdf_exponential_P (const double x, c
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  }  }
# Line 47  gsl_cdf_exponential_Q (const double x, c Line 47  gsl_cdf_exponential_Q (const double x, c
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  }  }
# Line 60  gsl_cdf_exponential_Q (const double x, c Line 60  gsl_cdf_exponential_Q (const double x, c
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  }  }

Legend:
Removed from v.1.1.2.1  
changed lines
  Added in v.1.1.2.2

savannah-hackers-public@gnu.org
ViewVC Help
Powered by ViewVC 1.1.26