colormap hsv;

x=linspace(-1,1.5);
y=x;

[X,Y]=meshgrid(x,y);
Z = (X-1).^2+10*(X.^2-Y).^2;

contourf(X,Y,Z,150);
axis([-1 1.5 -1 1.5]);
