
subplot(211);
y = -5:0.01:5;
p = (1.0/sqrt(4*pi))*(exp(-(y.^2)/4));
plot(y,p);
xlabel('y');
ylabel('p(y)');
axis([-5,5,-0.2,1.0]);

subplot(212);
y = 1.6;
x = -5:0.01:5;
p = (1.0/sqrt(pi))*(exp(-((x-y/2).^2)));
plot(x,p);
xlabel('x_1');
ylabel('p(x_1 | y=1.6)');
axis([-5,5,-0.2,1.0]);


