% Plots (x-1)^6 near x=1 with increasing refined scale
%
%
%
% Two ways to evaluate:
% 1. (x-1).^6
% 2. x.^6 - 6*x.^5 + 15*x.^4 - 20*x.^3 + 15*x.^2 -6*x + 1

function ZoomInPoly(delta)

if ~exist('delta','var')
    error('deltaIn must be a scalar')
end

  x=linspace(1-delta,1+delta,100)';   % Make a 100x1 vector over the range [1-delta,1+delta]
  %y=x.^6 - 6*x.^5 + 15*x.^4 - 20*x.^3 + 15*x.^2 -6*x + 1;
  y=(x-1).^6
  
  figure(1)
  plot(x,y,x,zeros(1,100))
  axis([1-delta 1+delta -max(abs(y)) max(abs(y))])
  
  title(sprintf('Zoomed in to delta=%e either side of 1',delta))