function dist=mahalan(X,mu,C)
% MAHALAN 付竅扼疇厚膠 芭府甫 拌魂茄促. 
%
% 儈過:
%  dist=mahalan(X,mu,C)
%
% 汲疙:
%  青紡 X客 氦磐 mu 弊府絆 青紡 C牢 痢甸 埃狼 付竅扼疇厚膠 芭府甫 拌魂茄促. 
%    dist = (x-mu)'*inv(C)*(x-mu),
%  咯扁輯 x綽 X狼 凱氦磐撈促. 
%
% 澇仿:
%  X [dim x num_data] 澇仿 單撈磐.
%  mu [dim x 1] 氦磐 mu.
%  C [dim x dim] 青紡 C.
%
% 免仿:
%  dist [1 x num_data] 付竅扼疇厚膠 芭府.
%


[dim, num_data] = size( X );

XC = X - repmat(mu,1,num_data);
dist= sum((XC'*inv( C ).*XC')',1);

return;