1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
#! /usr/bin/env python
import openturns as ot
nMin = 2
nMax = 4
pMin = 0.1
pMax = 0.9
nP = 2
alphaMin = 0.1
alphaMax = 0.9
nAlpha = 2
for n in range(nMin, nMax):
for iP in range(nP + 1):
p = pMin + (pMax - pMin) * iP / nP
for iAlpha in range(nAlpha + 1):
alpha = alphaMin + (alphaMax - alphaMin) * iAlpha / nAlpha
kF = ot.DistFunc.kFactor(2 * n, p, alpha)
print("kFactor(%d, %.12g, %.12g)=%.6g" % (2 * n, p, alpha, kF))
mMin = 1
mMax = 3
for n in range(nMin, nMax):
for m in range(mMin, mMax):
for iP in range(nP + 1):
p = pMin + (pMax - pMin) * iP / nP
for iAlpha in range(nAlpha + 1):
alpha = alphaMin + (alphaMax - alphaMin) * iAlpha / nAlpha
kF = ot.DistFunc.kFactorPooled(2 * n, m, p, alpha)
print(
"kFactorPooled(%d, %.12g, %.12g, %.12g)=%.6g"
% (2 * n, m, p, alpha, kF)
)
|
