r1 = sqrt(rand(90,1)); % radius
t1 = 2*pi*rand(90,1); % angle
data1 = [r1.*cos(t1), r1.*sin(t1)]; % points

r2 = sqrt(3*rand(10,1)+1); % radius
t2 = 2*pi*rand(10,1); % angle
data2 = [r2.*cos(t2), r2.*sin(t2)]; % points

data3 = [data1;data2];
theclass = cat(1,repmat('A',[90,1]),repmat('B',[10,1]));

idx = randperm(size(data3,1));
data3 = data3(idx,:);
theclass = theclass(idx,:);

Y = [A;A;B;D;B;A;C;A;A;A;D, ... ];