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

L = [ 1 2 3; 4 5 6; 7 8 9 ];

b1 = [ 1 1 1 ].'; b2 = [ -2 1 1 ].'; b3 = [ 0 1 -1 ].';
% d1 = b1/sqrt( ( b1.' * b1 ) ); d2 = b2/sqrt( ( b2.' * b2 ) ); d3 = b3/sqrt( ( b3.' * b3 ) ); 
d1 = b1; d2 = b2; d3 = b3; 

Ld1 = L*d1
Ld2 = L*d2
Ld3 = L*d3

pause

( d1.' )*Ld1
( d1.' )*Ld2
( d1.' )*Ld3

pause

( d2.' )*Ld1
( d2.' )*Ld2
( d2.' )*Ld3

pause

( d3.' )*Ld1
( d3.' )*Ld2
( d3.' )*Ld3

pause

d1 * d1.'

% Check:
(15/3)*( d1 * d1.' ) + (3/6)*( d1 * d2.' ) - (3/6)*( d1 * d3.' ) + ...
(9/6) *( d2 * d1.' ) - (9/6)*( d3 * d1.' )