prtab.m

          x = cumprod(1.0e-4*ones(5, 1)); 
          a = [x, exp(x), exp(-x), (exp(x)- exp(-x))./(2*x)]; 
          fprintf('%15s %15s %15s %15s\n', 'x', 'exp(x)', 'exp(-x)', 'driv');
          fprintf('%15.4e %15.4e %15.4e %15.4e\n', a');