function A = area(y, dx, N, r, R)
    % TODO:
    % - commentare i codici passaggio per passaggio
    % - rendere N, r e R parametri globali
    A = (y(1) + r) .* sqrt(dx^2 + (y(1) - r)^2);
    for k=1:N-2
        A = A + (y(k+1) + y(k)) .* sqrt(dx^2 + (y(k+1) - y(k))^2);
    end
    A = pi * (A + (R + y(end)) .* sqrt(dx^2 + (R - y(end))^2));
    return