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