The Maple commands used to produce this plot are as follows:
z:=exp(-x^2-2*y^2); plotlist:=[implicitplot({z=exp(-1),z=exp(-2),z=exp(-3),z=exp(-4), z=exp(-5),z=exp(-6)},x=-2..2,y=-2..2,grid=[40,40], scaling='CONSTRAINED')]: plotlist:=[op(plotlist),plot(x^2,x=-2..2,-2..2,color='GREEN'), plot(-x^2,x=-2..2,-2..2,color='GREEN')]: for i from 2 to 5 do plotp:=plot(i*x^2,x=-2..2,-2..2,color='GREEN'): plotmp:=plot(-i*x^2,x=-2..2,-2..2,color='GREEN'): plotm:=plot(1/i*x^2,x=-2..2,-2..2,color='GREEN'): plotmm:=plot(-1/i*x^2,x=-2..2,-2..2,color='GREEN'): plotlist:=[op(plotlist),plotp,plotm,plotmp,plotmm]: od: display(plotlist);