%   Figure 8.22 conjugate gradient paths

    [X,Y] = meshgrid(-2:.2:2, -2:.2:2);
        Z = X .^2+(1.5*Y).^2;
        figure;
        contour(X,Y,Z)
        colormap(gray)
        axis([-2 2 -2 2])
        axis square