function ud=p(x,mu,ssq,rho)

%function p=mynorm(x,mu,ssq,rho)
x=x-mu;
ud=exp( -1/(2*ssq*(1-rho^2)) * (x(1)^2-2*rho*x(1)*x(2)+x(2)^2) );