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

clc; close all; 

var_scale = 1e-2;
cA = mvnrnd( [ -3, +2 ], var_scale*eye(2), 20 );
cB = mvnrnd( [ +1.5, +2 ], var_scale*eye(2), 20 );
cC = mvnrnd( [ +3, +1 ], var_scale*eye(2), 20 ); 
X = [ cA; cB; cC ];

N = size(X,1); 

scatter( X(:,1), X(:,2), 25, ones(N,1), 'filled' );

labels = BSAS( X, sqrt(2), 3 ); 
scatter( X(:,1), X(:,2), 25, labels, 'filled' );
if( 1 ) saveas( gcf, '../../WriteUp/Graphics/Chapter12/dup_fig_12_1_q_3', 'eps' ); end

labels = BSAS( X, sqrt(2), 2 ); 
scatter( X(:,1), X(:,2), 25, labels, 'filled' );
if( 1 ) saveas( gcf, '../../WriteUp/Graphics/Chapter12/dup_fig_12_1_q_2', 'eps' ); end