%
%  Matlab script to generate Sauer's table on page 17 of the textbook.
%
x = logspace(0,-12,13);
E1 = (1-cos(x))./(sin(x).^2);
E2 = 1./(1+cos(x));
fid=fopen('SauerP17Table.txt','w');
fprintf(fid,'----------------------------------------------------------\n');
fprintf(fid,'|        x         |        E1        |        E2        |\n');
fprintf(fid,'----------------------------------------------------------\n');
for k = 1:length(x)
  fprintf(fid,'| %1.14f | %1.14f | %1.14f |\n',x(k),E1(k),E2(k));
end
fprintf(fid,'----------------------------------------------------------\n');

