function f = multiGauss(x1,x2,mx,Cx)

for i=1:length(x1),
    for j=1:length(x2),
        f(i,j)=(1/(2*pi*det(Cx)^1/2))*exp((-1/2)*([x1(i) x2(j)]-mx')*inv(Cx)*([x1(i);x2(j)]-mx));
    end
end