[x, y] = meshgrid(linspace(-2, 2, 101), linspace(-2, 2, 101));
z = (x.^2).*(y.^2) - x.*y + x + 2*y -2;
contour(x, y, z);