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

close all; clc; clear; 

X = [ 1, 2; 2, 3 ]; 
U = [ 1, 1; 1, -1 ] / sqrt(2); 

A00 = 0.5 * [ 1, 1; 1, 1 ]; 
A11 = 0.5 * [ 1, -1; -1, 1 ]; 
A01 = 0.5 * [ 1, -1; 1, -1 ]; 
A10 = 0.5 * [ 1, -1; 1, -1 ]; 

% take the required innner products: 
%
t = conj(X) .* A00;
fprintf('Y(0,0)= %10.6f\n', sum( t(:) ) );

t = conj(X) .* A01;
fprintf('Y(0,1)= %10.6f\n', sum( t(:) ) );

t = conj(X) .* A10;
fprintf('Y(1,0)= %10.6f\n', sum( t(:) ) );

t = conj(X) .* A11;
fprintf('Y(1,1)= %10.6f\n', sum( t(:) ) );