[x, y] = meshgrid(-.5:0.05:0.5, -.5:.05:.5);
xdot = 3*x - y - x.^3;
ydot = x-3*y- y.^3;
quiver(x, y, xdot, ydot)