function d=distanza(P1,P2)

d=sqrt(sum((P1-P2).^2));

% --- codice equivalente ---
% d=0;
% for ii=1:3
%     d=d+(P1(ii)-P2(ii))^2; % operazione tra scalari.
% end
% d=sqrt(d);