function C=plot_normal(y,A)
d=length(y);
if d ~= 2
    return
end
range=5;
x1=-range:0.02:range;
x2=x1;
for i=1:length(x1)
    for j=1:length(x2)
    x=[x1(i) x2(j)]';
    C(j,i)=1/((2*pi).^(d/2)*sqrt(det(inv(A))));
    C(j,i)=C(j,i)*exp(-(x-y)'*A*(x-y)*(0.5));
    end
end
mesh(x1,x2,C);
return
 