function test_weka()
S1= [1,0.9;
    0.9,2];
S2 = [6,0.2;
     0.2,2];
M1 = [13,15];
M2 = [2,7];
        
X = [mvnrnd([13,15],S1,n);
      mvnrnd([2,15],S2,n)];
  
C = [zeros(n,1)+1;zeros(n,1)+2;
plot_distribution2(X,[1,2],C);