¿Cómo trazar la función de Bessel en tikz al dorso?

Question

Intenté usar algún código existente, pero parece que no funciona al dorso.

También se recomendó ejecutar esto en: pdflatex -shell-escapepero no encontré esta opción en los compiladores.

código:

\documentclass[12pt]{article}
\usepackage{pgfplots}
\usepackage{tikz}

\begin{document}
  \begin{tikzpicture}
    \begin{axis}[width=\textwidth, height=0.5*\textwidth, xlabel=$x$]
    \addplot+[id=parable,domain=0:20, samples=500, mark=none, width=2pt]
    gnuplot{besj0(x)} node[pin=95:{$J_0(x)$}]{};
    \addplot+[id=parable,domain=0:20, samples=500, mark=none, width=2pt, color=red]
    gnuplot{besj1(x)} node[pin=130:{$J_1(x)$}]{};
    \addplot+[id=parable2,domain=0:20, samples=500, mark=none, width=2pt, color=black]
    gnuplot{2*1/x*besj1(x)-besj0(x)} node[pin=-140:{$J_2(x)$}]{};
   \end{axis}
  \end{tikzpicture}
\end{document}