function example_wiggle(n);

f=inline('(1+x.^2).^(-1)','x');

x=linspace(-5,5,n);
y=f(x);

X=linspace(-5,5,200);
Y=f(X);
plot(X,Y,'b-','LineWidth',3);

DD=DDtable(x,y);
c=diag(DD);
hold on;
plotnewtoninterp(x,y,c,'r-','LineWidth',2);