% Calcola (1-x)^6 con le due formule:
% y1=(1-x)^6
% y2=1-6x+15x^2-20x^3+15x^4-6x^5+x^6
% e confronta y1 e y2 in un intorno di uno
%
k=0;
for delta=[0.1, 0.01, 0.005, 0.0025]
    h=delta/100;
    x=1-delta:h:1+delta;
    y1=(1-x).^6;
    y2=1-6*x+15*x.^2-20*x.^3+15*x.^4-6*x.^5+x.^6;
    k=k+1;
    subplot(2,2,k)
    plot(x,y1,'r',x,y2,'g')
    axis([1-delta 1+delta -max(abs(y2)) max(abs(y2))])
end