clear all; close all;

y = 5*linspace(-1,1);

numer =  1*(normcdf(1-y)-normcdf(-y))...
        -1*(normcdf(-y)-normcdf(-1-y));
denom = normcdf(y+1)-normcdf(y-1);

% plot(y,th);
% plot(y,denom);
plot(y,numer./denom);
grid on;
set(gca,'FontSize',16);
xlabel('Noisy measurement, y');
ylabel('{\Theta}(y)');
title('EE 278B, PS#2, Problem 8b');