function Y = gencluster(m,c,n)
% generate cluster of n vectors with mean m and variance c

X = randn(n,2);
M = ones(n,1);
T = [[X] [M]];

% y1 = a*x1 + 0*x2 +1*m1 ; y2 = b*x1 + d*x2 + 1*m2;

a = sqrt(c(1,1));
b = c(1,2)/a;
d = sqrt(c(2,2) - b^2); 

W = [[a;0.0] [b;d] ;[m]];
Y = (T*W);