function classes = classify_GMM(X, mulist, sigmalist, weightlist);

% 單撈磐 籠欽 X俊輯 阿 氦磐俊 措茄 海撈令狼 盒幅甫 薦青 
% 撈 盒幅甫 烹竅咯 GMM葛膽撈 瘤瀝竅綽 努貳膠 吝俊輯 竅唱肺 撈 氦磐甫 且寸
%
% 澇仿 牢磊:
% X: d*numpoints 單撈磐 青紡
% mulist : 阿 努貳膠俊 措茄 MUS甫 啊瘤絆 樂(lè)綽 府膠飄 
% sigmalist : 阿 努貳膠俊 措茄 SIGMAS甫 啊瘤絆 樂(lè)綽 府膠飄
% weightlist : 阿 努貳膠俊 措茄 WEIGHTS甫 啊瘤絆 樂(lè)綽 府膠飄 
% 
% 免仿 牢磊 :
% classes : 阿 單撈磐 痢俊 且寸等 努貳膠甫 啊瘤絆 樂(lè)綽 凱 氦磐 

  [d,numpoints] = size(X);
  
  %if (is_list(mulist) & is_list(sigmalist) & is_list(weightlist)) == 0
  %   error('classify_GMM(X, mulist, sigmalist, weightlist) : lists should be lists!');
  %end    
      
  numclasses = length(mulist);
  
  %if (length(weightlist) == length(sigmalist) & length(sigmalist) == numclasses)
  %   error('List lengths inconsistent!');
  %end    
  
  bestloglike = loglikeGMM(X,mulist{:,1},sigmalist{:,1},weightlist{:,1});
  classes = ones(1,numpoints);
  for c = 2:numclasses
    newloglike = loglikeGMM(X,mulist{:,c},sigmalist{:,c},weightlist{:,c});
    better_idx = find(newloglike > bestloglike);
    bestloglike(better_idx) = newloglike(better_idx);
    classes(better_idx) = c;
  end