clc;  close all; clear all;

X = [27 19; 11 46; 25 15; 36 27; 35 25; 10 43; 11 44; 36 24; 26 14; 26 14; 9 45; 33 23; 27 16; 10 47];

[n, m] = size(X);

k = 3;
M = [10 45; 25 20; 35 25];
W = [1/k 1/k 1/k];
p = zeros(n, k);
g = zeros(n, k);
for i = 1:n
   for j = 1:k
       C{j} = cov(X, 1);
       xm = X(i,:)-M(j,:);
       fr = 1 / (((2*pi)^(m/2))*sqrt(det(C{j})));
       ee = exp(-0.5 * xm * (C{j}^-1) * xm');
       p(i, j) = fr * ee;
   end

   for j = 1:k
       wp = W(j)*p(i,j);
       wpSum = sum(W(j)*p(i,:));
       g(i, j) = wp / wpSum;
   end

end