function P3()
% Plots the error in MyErf across [-7,7]
x = linspace(-7,7,100)';
err = abs(erf(x) - MyErf(x));
semilogy(x,err+eps)
xlabel('x','fontsize',14)
ylabel('| erf(x) - MyErf(x) |','fontsize',14)
shg