%-----------------------------------------------------------------------
% Script File: ExpPlot
% Examine the function f(x)=((1+x/24)(1-x/12+x^2/384))^8
% as an approximation to exp(x) across [0,1]
close all
x = linspace(0,50,1000);
num = 1+x/24;
denom = 1 - x/12 + (x/384).*x;
quot = num./denom;
y = quot.^8;
plot(x,y',x,exp(x),'o')
  axis([0 50 0 1000000])
