Missing contours in Beale's function

I want to plot the Beale’s function and see the contours close to the minimum.
I have not been able to display those contours
Here is my example:

documentclass[crop]{standalone}
usepackage[pdftex]{graphicx}
usepackage{pgfplots,tikz}
usepackage{tikz-3dplot}
pgfplotsset{compat=newest}

    begin{document}
    begin{tikzpicture}
      begin{axis}[]
        
        %addplot3[surf,shader=interp, opacity=1.0] 
        %{(1.5-x+x*y)^2 + (2.25 - x + x*y^2)^2 + (2.25 -x + x*y^2)^2 + (2.65 -x + x*y^3)^2};
        addplot3[contour gnuplot={number=100, labels=false},thick, samples=10] 
        {(1.5-x+x*y)^2 + (2.25 - x + x*y^2)^2  + (2.65 -x + x*y^3)^2};
      end{axis}
    end{tikzpicture}
end{document}

Here is the Figure that I am getting:

Here is a link with a guide for the contours I would like to see