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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
|
module m_gsl_sf_erf
implicit none
interface
! Complementary Error Function
! erfc(x) := 2/Sqrt[Pi] Integrate[Exp[-t^2], {t,x,Infinity}]
!
! exceptions: none
integer function gsl_sf_erfc_e (x, result)
use m_gsl_sf_result
implicit none
double precision x
type(gsl_sf_result) result
end function gsl_sf_erfc_e
double precision function gsl_sf_erfc (x)
implicit none
double precision x
end function gsl_sf_erfc
! Log Complementary Error Function
!
! exceptions: none
integer function gsl_sf_log_erfc_e (x, result)
use m_gsl_sf_result
implicit none
double precision x
type(gsl_sf_result) result
end function gsl_sf_log_erfc_e
double precision function gsl_sf_log_erfc (x)
implicit none
double precision x
end function gsl_sf_log_erfc
! Error Function
! erf(x) := 2/Sqrt[Pi] Integrate[Exp[-t^2], {t,0,x}]
!
! exceptions: none
integer function gsl_sf_erf_e (x, result)
use m_gsl_sf_result
implicit none
double precision x
type(gsl_sf_result) result
end function gsl_sf_erf_e
double precision function gsl_sf_erf (x)
implicit none
double precision x
end function gsl_sf_erf
! Probability functions:
! Z(x) : Abramowitz+Stegun 26.2.1
! Q(x) : Abramowitz+Stegun 26.2.3
!
! exceptions: none
integer function gsl_sf_erf_Z_e (x, result)
use m_gsl_sf_result
implicit none
double precision x
type(gsl_sf_result) result
end function gsl_sf_erf_Z_e
integer function gsl_sf_erf_Q_e (x, result)
use m_gsl_sf_result
implicit none
double precision x
type(gsl_sf_result) result
end function gsl_sf_erf_Q_e
double precision function gsl_sf_erf_Z (x)
implicit none
double precision x
end function gsl_sf_erf_Z
double precision function gsl_sf_erf_Q (x)
implicit none
double precision x
end function gsl_sf_erf_Q
#if 0
/* Does not exist in older versions of GSL*/
! Hazard function, also known as the inverse Mill's ratio.
!
! H(x) := Z(x)/Q(x)
! = Sqrt[2/Pi] Exp[-x^2 / 2] / Erfc[x/Sqrt[2]]
!
! exceptions: GSL_EUNDRFLW
integer function gsl_sf_hazard_e (x, result)
use m_gsl_sf_result
implicit none
double precision x
type(gsl_sf_result) result
end function gsl_sf_hazard_e
double precision function gsl_sf_hazard (x)
implicit none
double precision x
end function gsl_sf_hazard
#endif
end interface
end module m_gsl_sf_erf
|