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

set terminal svg
set output "BesselK_plot.svg"
set title "Modified Bessel function of second kind"
set xlabel "x"
set ylabel "K_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
plot "besselk.dat" using 1:2 title "K_0(x)" with lines linewidth 2, \
     "besselk.dat" using 1:3 title "K_1(x)" with lines linewidth 2, \
     "besselk.dat" using 1:4 title "K_2(x)" with lines linewidth 2

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

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