% X-hr of 10/1/08, Jon Brown's task.  code by Zaven and Alex
% needs taylor_error function

xarr=0:0.0001:2;
n = 20;
e = zeros(1,n);
for i=1:n, e(i)=taylor_error(i,xarr); end
%e
plot(1:n,e, '+');
semilogy(1:n,e, '+');

% 3d plot of pt-wise error vs N
xarr=0:0.01:2;
n = 20;
e = zeros(numel(xarr),n);
for i=1:n, [out,a,b]=taylor_error(i,xarr); e(:,i) = b; end
mesh(xarr, 1:n, log10(e'));
xlabel('x'); ylabel('N');