
% Numero punti
n = 5;         

% Ampiezza intervalli
h = 2/(n-1);      

% Ascisse punti
x(1) = -1;  

for k = 2:n
    x(k) = x(k-1) + h;
end

% Ordinate punti
for k = 1:n    
    y(k) = 1 - x(k)^2;
end

% Approssimazione lunghezza
A = 0;

for k = 1:(n-1)   
    A = A + (y(k) + y(k+1))*(x(k+1) - x(k))/2;
end

% Calcola l'errore
A_ex = 4/3;
Err = abs(A_ex - A);

% Stampa i risultati a terminale
fprintf("\n   n = %d,   A = %.5f,   Err = %.6f\n", n, A, Err);
