%
% Written by:
% --
% John L. Weatherwax                2009-04-21
%
% email: wax@alum.mit.edu
%
% Please send comments and especially bug reports to the
% above email address.
%
%-----

% EPage 152

close all; clc; clear all; 

addpath('../../BookCode');

a = [ [ 5, 4, 1, 1 ]; [ 4, 5, 1, 1 ]; [ 1, 1, 4, 2 ]; [ 1, 1, 2, 4 ] ];
b = [ 1 ; 2 ; 3 ; 4 ];
xdash = solvercg(a,b,4,1.e-6)
norm( b - a * xdash )