Sigma = [1, 0; 0, 1];

mu1 = [1, -1];

x1 = mvnrnd(mu1, Sigma, 200);

mu2 = [5.5, -4.5];

x2 = mvnrnd(mu2, Sigma, 200);

mu3 = [1, 4];

x3 = mvnrnd(mu3, Sigma, 200);

mu4 = [6, 4.5];

x4 = mvnrnd(mu4, Sigma, 200);

mu5 = [9, 0.0];

x5 = mvnrnd(mu5, Sigma, 200);

% obtain the 1000 data points to be clustered

X = [x1; x2; x3; x4; x5];

% Show the data point

plot(x1(:,1), x1(:,2), "r."); hold on;

plot(x2(:,1), x2(:,2), "b.");

plot(x3(:,1), x3(:,2), "k.");

plot(x4(:,1), x4(:,2), "g.");

plot(x5(:,1), x5(:,2), "m.");

save myX %將X儲存到檔案中，在其他檔案中load就可以了