function compParams = Mstep(X, p)
  [n,d] = size(X) ;
  k = size(p,2) ;
  compParams = [] ;
  for i=1:k
    prior =
    mu =
    Sigma = 
    compParams = [compParams, struct('prior',prior,'mean',mu,'cov',Sigma)];
  end;

    