Connect via MCP →

Enter Calculation

Formula

Advertisement

Results

First table value I_v(x) at x = initial value
1
51 rows generated
x I_v(x)
0 1
0.1 1.0025015629
0.2 1.0100250278
0.3 1.0226268794
0.4 1.0404017822
0.5 1.0634833707
0.6 1.0920453643
0.7 1.1263030183
0.8 1.1665149229
0.9 1.2129851657
1 1.2660658778
1.1 1.3261601837
1.2 1.3937255841
1.3 1.4692777979
1.4 1.5533950997
1.5 1.6467231898
1.6 1.7499806397
1.7 1.8639649621
1.8 1.9895593566
1.9 2.1277401941
2 2.2795853023
2.1 2.4462831294
2.2 2.6291428636
2.3 2.8296056006
2.4 3.049256658
2.5 3.2898391441
2.6 3.5532689042
2.7 3.8416509766
2.8 4.1572977035
2.9 4.5027486613
3 4.8807925859
3.1 5.2944914897
3.2 5.7472071872
3.3 6.2426304652
3.4 6.7848131604
3.5 7.3782034322
3.6 8.0276845471
3.7 8.7386175242
3.8 9.5168880261
3.9 10.3689579167
4 11.3019219521
4.1 12.323570116
4.2 13.4424561633
4.3 14.6679729918
4.4 16.0104355249
4.5 17.4811718556
4.6 19.0926234795
4.7 20.8584555266
4.8 22.7936779931
4.9 24.9147790758
5 27.2398718236

What this calculator does

This tool tabulates the modified Bessel function of the first kind, written \(I_v(x)\), for a fixed real order \(v\) over a sequence of \(x\) values. You provide the order, an initial \(x\), an increment (step), and how many rows to generate; the calculator builds the list \(x_i = \text{start} + i\cdot\text{step}\) and evaluates \(I_v(x_i)\) at each point, returning both a table and a graph. It is a pure-mathematics special-function tool and applies universally (no regional rules or units).

Family of monotonically increasing exponential-like curves for the modified Bessel function of the first kind for several orders
Modified Bessel functions I_v(x) for orders v = 0, 1, 2, 3 rise rapidly with x.

The formula

The modified Bessel function \(I_v(x)\) solves the modified Bessel equation $$x^2 y'' + x y' - (x^2 + v^2)y = 0.$$ It is computed here from its power series:

$$I_{v}(x) = \sum_{k=0}^{\infty} \frac{1}{k!\;\Gamma(k+v+1)}\left(\frac{x}{2}\right)^{v+2k}.$$

The factorial and the Gamma function let \(v\) be any real number. For numerical stability each term is evaluated in log space using a Lanczos approximation of \(\ln\Gamma\), then summed until the terms are negligible.

How to use it

Enter Order v (e.g. 0, 1, or 2.5), the Initial value of x, the Increment added to \(x\) for each row, and the Number of repetitions (rows). Click calculate to get a two-column table of \(x\) and \(I_v(x)\) plus a graph over the same range.

Worked example

With \(v = 0\), start = 0, step = 0.5, count = 5 you get \(x = 0, 0.5, 1, 1.5, 2\) and:

$$I_{0}(0) = 1,\quad I_{0}(0.5) \approx 1.0634834,\quad I_{0}(1) \approx 1.2660658,\quad I_{0}(1.5) \approx 1.6467232,\quad I_{0}(2) \approx 2.2795853.$$ These match standard reference tables.

FAQ

Can the order be negative or non-integer? Yes. For negative integer order the identity \(I_{-n}(x) = I_{n}(x)\) is used. Non-integer \(v\) is supported for \(x \ge 0\); for \(x < 0\) with non-integer \(v\) the value is complex, so NaN is returned.

Why does \(I_v(x)\) grow so fast? Unlike the oscillating ordinary Bessel \(J_v\), the modified function grows roughly like \(e^{x}/\sqrt{2\pi x}\) for large \(x\), so large \(x\) may overflow to infinity.

What is \(I_v(0)\)? \(I_{0}(0) = 1\), and \(I_{v}(0) = 0\) for \(v > 0\).

Last updated: