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

Diff of /gsl/cdf/t.c

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

revision 1.1.2.6 by jstover, Fri Feb 7 16:20:37 2003 UTC revision 1.1.2.7 by jstover, Mon Feb 10 04:02:45 2003 UTC
# Line 210  int gsl_cdf_t_e ( double t, double nu, g Line 210  int gsl_cdf_t_e ( double t, double nu, g
210              {              {
211                q = M_PI_2 * (atan(y) + ck * y / b);                q = M_PI_2 * (atan(y) + ck * y / b);
212              }              }
213            if ( tail == GSL_CDF_UPPER )            if ( ((tail == GSL_CDF_UPPER ) && (t > 0.0))||((tail==GSL_CDF_LOWER)&&(t<=0.0)))
214              {              {
215                result->val = (1.0 - q)/2;                result->val = (1.0 - q)/2;
216              }              }
# Line 228  int gsl_cdf_t_e ( double t, double nu, g Line 228  int gsl_cdf_t_e ( double t, double nu, g
228             */             */
229            u = cornish_fisher ( t, nu );            u = cornish_fisher ( t, nu );
230            q = gsl_cdf_gauss ( u, GSL_CDF_UPPER );            q = gsl_cdf_gauss ( u, GSL_CDF_UPPER );
231            if ( tail == GSL_CDF_UPPER )            if ( ((tail == GSL_CDF_UPPER )&&(t>0.0))||((tail==GSL_CDF_LOWER)&&(t<=0.0)))
232              {              {
233                result->val = q;                result->val = q;
234              }              }
# Line 282  int gsl_cdf_t_e ( double t, double nu, g Line 282  int gsl_cdf_t_e ( double t, double nu, g
282        diff = (lg2.val) - (lg1.val);        diff = (lg2.val) - (lg1.val);
283        p *= pow(y,nu) * exp (diff) / sqrt (M_PI);        p *= pow(y,nu) * exp (diff) / sqrt (M_PI);
284    
285        if ( tail == GSL_CDF_UPPER )        if ( ((tail == GSL_CDF_UPPER )&&(t>0.0))||((tail==GSL_CDF_LOWER)&&(t<=0.0)))
286          {          {
287            result->val = p;            result->val = p;
288          }          }

Legend:
Removed from v.1.1.2.6  
changed lines
  Added in v.1.1.2.7

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