/usr/share/openturns/validation/dProb.txt is in openturns-validation 1.5-7build2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 | > restart;
> dbeta:=beta*(sqrt(1-2*ln(epsilon)/beta^2)-1);
// 2 ln(epsilon)\1/2 \
dbeta := beta ||1 - -------------| - 1|
|| 2 | |
\\ beta / /
> phi:=proc(x)
> 1/(2*Pi)^(N/2)*exp(-x*x/2)
> end:
> PHI:=proc(x)
> int(phi(t),t=-infinity..x)
> end:
> khi2:=proc(x,n)
> 1/GAMMA(n/2)*(1/2)^(n/2)*x^(n/2-1)*exp(-1/2*x)
> end:
> KHI2:=proc(x,n)
> int(khi2(t,n),t=0..x)
> end:
>
> sol1:=solve(phi(beta+dbeta)/phi(beta)=epsilon,dbeta);
> sol2:=solve(PHI(-beta-dbeta)/PHI(-beta)=epsilon,dbeta);
> sol3:=solve((1-KHI2(beta+dbeta,6))/(1-KHI2(beta,6))=epsilon,dbeta);
2 1/2
sol1 := -beta + (beta - 2 ln(epsilon)) ,
2 1/2
-beta - (beta - 2 ln(epsilon))
/
| 1/2
sol2 := -1/2 |2 beta
\
1/2 \
2 beta |
- 2 RootOf(-erf(_Z) + 1 + epsilon erf(---------) - epsilon)|
2 /
1/2
2
sol3 :=
> plot(subs(beta=5,n=2,[sol1[1],sol2,sol3]),epsilon=0.001..1);
> series(sol1[1]/sol2,epsilon);
Error, (in series/RootOf) unable to compute series
> with(stats);
Warning, these names have been redefined: anova, describe, fit,
importdata, random, statevalf, statplots, transform
[anova, describe, fit, importdata, random, statevalf, statplots,
transform]
>
|