%% 
clear all
li = linspace(-10, 10, 1001);
[x,y] = meshgrid(li, li);
contour(x, y, y.^3 - x.*y + x, 'linewidth', 4)
axis equal;
