Perror.m

x=[0:0.01:100];
for i=100:10000
Pe(i)=0.01*sqrt(2/pi)*sum(exp(-0.5*x(i:10000).^2));
end