itload('../mixef_test.it');

figure(1);
hold off;
plot(Smp(1,:), Smp(2,:),'.');
hold on
xv=-5:0.1:5;
X = ones(size(xv))'*xv;
Y = xv'*ones(size(xv));

contour(X,Y,PPdf)