function d = mhdist(X, Z)
% MHDIST1: Mahalanobis distance under estimated parameters
% for the Gaussian.

mu = mean(Z, 1);
sg = cov(Z);    

d = mhdist1(X, mu, sg);

return

