
% 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
L = 0;

for k = 1:(n-1)   
    L = L + sqrt( (x(k+1)-x(k))^2 + (y(k+1)-y(k))^2 );
end

% Calcola l'errore
L_ex = 2.9578857;
Err = abs(L_ex - L);

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