%Section 2.2 computer problem #8 

format long

n = 15;     %number of iterations

f = @(x) (x-sin(x))/x^3;
g = @(x) 1/6 - x^2/120 + x^4/5040;

disp('naive estimates of f are...')

for i=0:n
    f(10^(-i))
end

disp('accurate estimates of f near 0 are...')

for i=0:n
    g(10^(-i))
end
