Doxygen Source Code Documentation
cdf_71.c File Reference
#include "cdflib.h"Go to the source code of this file.
Functions | |
| double | spmpar (int *i) |
Function Documentation
|
|
Definition at line 2 of file cdf_71.c. References i, ibeta(), ipmpar(), and spmpar(). Referenced by bratio(), cdfbet(), cdfbin(), cdfchi(), cdff(), cdfgam(), cdfnbn(), cdfnor(), cdfpoi(), cdft(), cumnor(), gaminv(), gratio(), psi(), spmpar(), and Xgamm().
00030 {
00031 static int K1 = 4;
00032 static int K2 = 8;
00033 static int K3 = 9;
00034 static int K4 = 10;
00035 static double spmpar,b,binv,bm1,one,w,z;
00036 static int emax,emin,ibeta,m;
00037 /*
00038 ..
00039 .. Executable Statements ..
00040 */
00041 if(*i > 1) goto S10;
00042 b = ipmpar(&K1);
00043 m = ipmpar(&K2);
00044 spmpar = pow(b,(double)(1-m));
00045 return spmpar;
00046 S10:
00047 if(*i > 2) goto S20;
00048 b = ipmpar(&K1);
00049 emin = ipmpar(&K3);
00050 one = 1.0;
00051 binv = one/b;
00052 w = pow(b,(double)(emin+2));
00053 spmpar = w*binv*binv*binv;
00054 return spmpar;
00055 S20:
00056 ibeta = ipmpar(&K1);
00057 m = ipmpar(&K2);
00058 emax = ipmpar(&K4);
00059 b = ibeta;
00060 bm1 = ibeta-1;
00061 one = 1.0;
00062 z = pow(b,(double)(m-1));
00063 w = ((z-one)*b+bm1)/(b*z);
00064 z = pow(b,(double)(emax-2));
00065 spmpar = w*z*b*b;
00066 return spmpar;
00067 } /* END */
|