\documentclass[pstricks]{standalone}
\usepackage{pst-contourplot,multido,pst-plot}
\begin{document}
% équation donnée par Paul Bourke dans :
% http://paulbourke.net/papers/conrec/
\begin{pspicture}(-5,-5)(5,5)
\multido{\r=0.25+0.25,\n=0.0+0.1}{11}{%
\pstVerb{/isovalue \r\space def}%
\psContourPlot[algebraic,unit=2.5,a=0.02,function=1/((y^2+x^2-0.71)^2+4*y^2*(x-0.842)^2)-isovalue,linecolor={[hsb]{\n,1,1}}](-2,-2)(2,2)
}%
\end{pspicture}

\end{document}