% 
% 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.
% 
%-----

X1 = [ 0.2, 0.7; 0.3, 0.3; 0.4, 0.5; 0.6, 0.5; 0.1, 0.4 ];
n1 = size(X1,1); 
X2 = [ 0.4, 0.6; 0.6, 0.2; 0.7, 0.4; 0.8, 0.6; 0.7, 0.5; ]; 
n2 = size(X2,1); 

X = [ X1; X2 ]; 

% append a one as the last column of X:
%
X = [ X, ones(size(X,1),1) ]; 

XTX = X' * X;
disp( XTX ); 

XTy = X' * y;
disp( XTy );

y = [ 1*ones(n1,1); -1*ones(n2,1) ]; 

what = ( XTX ) \ ( XTy );
disp( what );