a = linspace(-5,5);
n = length(a);
for i = 1:n
    P(i) = integral(@stdnormdens,-Inf,a(i));
end
plot(a,P,'k-') ; grid
xlabel('z')
ylabel('Cumulative Probability')

