function [logp] = evalgauss(X,param)

[n,m] = size(X);
Sinv = inv(param.cov);

Xt = (X-ones(n,1)*param.mean); 
logp = -sum( (Xt*Sinv).*Xt, 2 )/2 - 0.5*log(det(param.cov))-0.5*m*log(2*pi);

