Erro ao criar distribuição normal em PGFPlots

Question

Olá, quero criar uma figura distribuída normal. com sigma = 10 e mu =185;

Eu tenho o seguinte código baseado emDistribuição 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}

O problema é que sempre recebo o seguinte erro.

O parágrafo terminou antes de \pgfflt@readlowlevelfloat ser concluído
\begin{axis} na linha de entrada 480 finalizado por \end{document}. \end{document}
Faltando \endgroup inserido \end{document}

O que estou errado? O mais estúpido é que quando copio o original ele funciona bem e a minha versão menor não funciona…

Espero que alguém possa me ajudar

Obrigado