clc;
close all;
clear all;

sigma1 = 0.5*0.5;
sigma2 = 2*sigma1;
sigma3 = 0.5*sigma1;

N = 1000;

c{1}.x = sqrt(sigma1)*randn(N, 1);
c{1}.y = sqrt(sigma1)*randn(N, 1);

r1 = randn(N, 1);
r2 = randn(N, 1);
c{2}.x = sqrt(sigma2)*r1.*r1.*r2 -2;
c{2}.y = sqrt(sigma2)*r2.*r2.*r2 +5;

plot(c{1}.x, c{1}.y, 'b*', c{2}.x, c{2}.y, 'r*');
xlabel('Feature 1');
ylabel('Feature 2');
legend('Class A', 'Class B');



