Doxygen Source Code Documentation
Main Page Alphabetical List Data Structures File List Data Fields Globals Search
cdf_45.c
Go to the documentation of this file.00001 #include "cdflib.h"
00002 double dln1px(double *a)
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039 {
00040 static double p1 = -.129418923021993e+01;
00041 static double p2 = .405303492862024e+00;
00042 static double p3 = -.178874546012214e-01;
00043 static double q1 = -.162752256355323e+01;
00044 static double q2 = .747811014037616e+00;
00045 static double q3 = -.845104217945565e-01;
00046 static double dln1px,t,t2,w,x;
00047
00048
00049
00050
00051 if(fabs(*a) > 0.375e0) goto S10;
00052 t = *a/(*a+2.0e0);
00053 t2 = t*t;
00054 w = (((p3*t2+p2)*t2+p1)*t2+1.0e0)/(((q3*t2+q2)*t2+q1)*t2+1.0e0);
00055 dln1px = 2.0e0*t*w;
00056 return dln1px;
00057 S10:
00058 x = 1.e0+*a;
00059 dln1px = log(x);
00060 return dln1px;
00061 }