function [res]=gaussvalue(x,m,c)
t=(x-m)'*pinv(c)*(x-m);
t1=exp(-t/2);
t2=1/((2*pi)^(length(x)/2));
[E,V]=eig(c);
V=diag(V);
[P,Q]=find(V>0);
t3=1;
for i=1:1:length(P)
    t3=t3*V(P(i));
end
t3=1/(abs(t3)^(1/2));
res=t1*t2*t3;