  function P1(a)
% Looks at the quality of various Rqq(x) across the interval [-a,a];
x = linspace(-a,a,1000)';
for q = 1:10
   y = DiagPade(x,q);
   semilogy(x,abs(y-exp(x))+eps)
   title(sprintf('| exp(x) - R_{qq}(x) |, q = %1d',q),'fontsize',14)
   shg
   pause(1)
end
