Error al crear distribución normal en PGFPlots

Error al crear distribución normal en PGFPlots

Hola, quiero crear una figura distribuida Normal. con sigma = 10 y mu =185;

Tengo el siguiente código basado enDistribución normal bivariada

\usepackage{tikz}
\usepackage{pgfplots}
\begin{document}    
\begin{tikzpicture}[
        declare function={mu=185;}, %
        declare function={sigma=10;}, %
        declare function={normal(\m,\s)=1/(2*\s*sqrt(pi))*exp(-(x-\m)^2/(2*\s^2));}] %
        \begin{axis}[
        width=15cm,
        enlargelimits=false,
        grid=major,
        domain=-1:4,
        y domain=-1:4,
        samples=201,
        xlabel=$x_1$,
        ylabel=$P$,
        ]
        \addplot [domain=150:220,samples=201, samples y=0, thick, smooth] (x,4,{normal(mu,sigma)});
        \end{axis}
    \end{tikzpicture}
\end{document}

El caso es que cada vez aparece el siguiente error.

El párrafo finalizado antes de que \pgfflt@readlowlevelfloat se completara
\begin{axis} en la línea de entrada 480 terminada en \end{document}. \end{document}
Falta \endgroup insertado \end{document}

¿Qué hago mal? Lo estúpido es que cuando copio el original funciona bien y la mía menos grande no funciona…

Espero que alguien pueda ayudarme

Gracias

información relacionada