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

Diff of /gsl/cdf/paretoinv.c

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

revision 1.1.2.1 by bjg, Sun Jul 20 15:01:32 2003 UTC revision 1.1.2.2 by bjg, Sun Jul 20 18:51:44 2003 UTC
# Line 36  gsl_cdf_pareto_Pinv (const double P, con Line 36  gsl_cdf_pareto_Pinv (const double P, con
36        return b;        return b;
37      }      }
38    
39    x = b / pow(1-p, 1/a);    x = b * exp(-log1p(-P)/a);
40    
41    return x;    return x;
42  }  }
# Line 55  gsl_cdf_pareto_Q (const double Q, const Line 55  gsl_cdf_pareto_Q (const double Q, const
55        return b;        return b;
56      }      }
57    
58    x = b / pow(p, 1/a);    x = b * exp(-log(Q) / a);
59    
60    return x;    return x;
61  }  }

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