# Problems with plotting planar functions with Tikz

I tried to plot this graph using tikz and pgfplots. The code compiled just fine, but I got a peculiar result that I can't explain. My code.

\documentclass{article}

\usepackage{amsmath}
\usepackage{tikz}
\usepackage{pgfplots}

\begin{document}

\begin{tikzpicture}
\begin{axis}[xlabel=$x$,ylabel=$y$,
xmin=-6,xmax=6,ymin=-6,ymax=6, axis lines=center, axis equal]

\end{axis}
\end{tikzpicture}

\end{document}

This is the riddle I was presented with.

The shape my code produced (a line, strangely enough) perplexes me. I truly have no idea what I did wrong. Is there any internal process at fault? Can you please modify my code in such a way that the desired plot will actually show up? I'd be grateful if you (kind users who already mastered tikz and pgfplots) would help me. Thank you in advance.

I have seen this graph when I saw an introductory article on the small application GrafEq. The ideal behind is to compute the LHS and RHS using interval arithmetic. We cannot (at least it is pretty difficult to) do interval arithmetic in TeX. But we have the following for free

\documentclass[tikz]{standalone}

3 mul exch 3 mul exch % x y
2 copy                % x y x y
dup mul exch dup mul  % x y y² x²
add sin               % x y sin(x²+y²)
3 1 roll              % sin(x²+y²) x y
mul cos               % sin(x²+y²) cos(xy)
le                    % < or not <
{0 0 0}{1 1 1}ifelse  % black or white
}
3 mul exch 3 mul exch % x y
2 copy                % x y x y
dup mul exch dup mul  % x y y² x²
add sin               % x y sin(x²+y²)
3 1 roll              % sin(x²+y²) x y
mul cos               % sin(x²+y²) cos(xy)
sub dup mul           % (L-R)²
.01 le                % <.01 or not <.01
{0 0 0}{1 1 1}ifelse  % black or white
}

\begin{document}
\end{document}

# OR GNUPLOT

\documentclass[tikz]{standalone}
\usepackage{pgfplots}
\begin{document}
\begin{tikzpicture}
\begin{axis}[view={0}{90}]