eps = .5;
tol = 1e-9;
A = [1,1+eps; 1-eps , 1];
xtrue = [1.1;-1];
g = A*xtrue;
xcgn = CGNE1([0;0],A,g,1e-9,2);
[xcg,xicg,ricg] = CG1([0;0],A'*A,A'*g,tol,2);
[xsd,xisd,risd] = SD1([0;0],A,g,tol,10);

figure(1); clf ; 
subplot(1,2,1);hold on;
plot(xisd(1,:),xisd(2,:),'-b+');
plot(xicg(1,:),xicg(2,:),'--k+');
title('solution steps');

subplot(1,2,2);hold on;
plot(risd(1,:),risd(2,:),'-b+');
plot(ricg(1,:),ricg(2,:),'--k+');
title('gradient directions steps');