D=diag([1 1 10^(-2)]);
b=[1 1 0]';
x=[1.05 0.95 10]';
W=x;
for j=1:30,
r=b-D*x;
a=r'*r/(r'*D*r);
x=x+a*r;
W=[W x];
end
k=W(1,:);l=W(2,:);m=W(3,:);

% Plots
subplot(2,2,1), plot3(k,l,m,'--',k,l,m,'or'),...
   set(gca,'box','on'),...
   xlabel('x'), ylabel('y'),zlabel('z'),...
   grid, title('3D View')
subplot(2,2,2), plot(k,l,'--',k,l,'or'),...
   xlabel('x'), ylabel('y'),grid, title('X-Y Plane')
subplot(2,2,3), plot(k,m,'--',k,m,'or'),...
   xlabel('x'), ylabel('z'),grid, title('X-Z Plane')
subplot(2,2,4), plot(l,m,'--',l,m,'or'),...
   xlabel('y'), ylabel('z'),grid, title('Y-Z Plane')

