close all;
clear all;
clc;

% Vettore contenente i valori di delta richiesti
deltas = [1e-1,1e-2,5e-3,2.5e-3];

% Ciclo sui delta
for delta = deltas

    % Calcolo dei 100 punti equidistanti
    x = linspace(1-delta,1+delta,100);

    % Calcolo delle due espressioni per i 100 punti nell'intervallo
    % descritto
    y1 = (1-x).^6;
    y2 = x.^6 - 6*x.^5 + 15*x.^4 - 20*x.^3 + 15*x.^2 - 6*x + 1;

    % Plot delle due espressioni sovrapposte
    figure
    plot(x,y1,'r')
    hold on
    plot(x,y2,'b')
    legend('y1','y2')
    
    % Calcolo degli errori assoluti e relativi
    err_ass = abs(y1-y2);
    err_rel = abs(err_ass./y1);

    fprintf('\n \t Delta: %1.2e : ErrMaxAss: %1.3e  ErrMaxRel: %1.3e \n',delta,max(err_ass),max(err_rel))

end