Plot of the modified Bessel functions of 1st kind ${\displaystyle I_{0}(x)}$ in red, ${\displaystyle I_{1}(x)}$ in green and ${\displaystyle I_{2}(x)}$ in blue. The plot was produced running Gnuplot on the following code:

set terminal svg
set output "BesselI_plot.svg"
set title "Modified Bessel function of first kind"
set xlabel "x"
set ylabel "I_n(x)"
set xrange [0:4]
set yrange [0:4]
set xzeroaxis linetype -1 linewidth 0.5
set yzeroaxis linetype -1 linewidth 0.5
set xtics axis
set ytics axis
set key off
plot "besseli.dat" using 1:2 title "I_0(x)" with lines linewidth 2, \
     "besseli.dat" using 1:3 title "I_1(x)" with lines linewidth 2, \
     "besseli.dat" using 1:4 title "I_2(x)" with lines linewidth 2

the file "besseli.dat" contains the values of the Bessel functions and can be produced with the following Matlab commands:

t=0:0.001:4;
I = [t; besseli(0,t); besseli(1,t); besseli(2,t)];
I = I';
save -ascii "besseli.dat" I;