We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 8dea409 commit 435417dCopy full SHA for 435417d
src/besseli.jl
@@ -8,7 +8,7 @@
8
#
9
# Branch 1: x < 7.75
10
# besseli0 = [x/2]^2 P16([x/2]^2)
11
-# Branch 2: x >=
+# Branch 2: x >= 7.75
12
# sqrt(x) * exp(-x) * besseli0(x) = P22(1/x)
13
# where P16 and P22 are a 16 and 22 degree polynomial respectively.
14
@@ -23,8 +23,8 @@
23
24
25
# besseli1 = x / 2 * (1 + 1/2 * (x/2)^2 + (x/2)^4 * P13([x/2]^2)
26
27
-# sqrt(x) * exp(-x) * besseli0(x) = P22(1/x)
+# sqrt(x) * exp(-x) * besseli1(x) = P22(1/x)
28
# where P13 and P22 are a 16 and 22 degree polynomial respectively.
29
30
# Remez.jl is then used to approximate the polynomial coefficients of
0 commit comments