From: Maikel on
Hi,

I'm trying to get a contour plot of the following function

[q,p] = meshgrid(-1:0.1:1);
H = q.*p.^2;
contour(q,p,H);

The result is shown here:
http://www.maikelbosschaert.nl/math/contour_H.jpg

The problem is that the zero level set is wrong.
I'm missing the right line when H=0: (q,p)=(q,0) with q \in R.

I think this because of the interpolation method being used to create the contour plot data.

Is there a way to resolve this problem?

Thanks in advance,

Maikel