Wie stellt man die Bessel-Funktion in Tikz auf Overleaf dar?

Question

Ich habe versucht, vorhandenen Code zu verwenden, aber das scheint in Overleaf nicht zu funktionieren.

Außerdem wurde empfohlen, dies in auszuführen:, pdflatex -shell-escapeaber ich habe diese Option in den Compilern nicht gefunden.

Code:

\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}