for i=1:Ne
    PT(i,1)=det([1 x(T(i,1)) y(T(i,1)) z(T(i,1));
                 1 x(T(i,2)) y(T(i,2)) z(T(i,2));
                 1 x(T(i,3)) y(T(i,3)) z(T(i,3));
                 1 x(T(i,4)) y(T(i,4)) z(T(i,4))]);
    if (PT(i,1)<0)
        PT(i,1)=-PT(i,1);
        app=T(i,1);
        T(i,1)=T(i,2);
        T(i,2)=app;
    end
    if (PT(i,1)<TOLLERANZA)
       fprintf('Errore.... tetraedro con volume nullo\n'); 
    end
end

% VettBuoni=find(PT(:,1)>TOLLERANZA);
% T=T(VettBuoni,:);
% PT=PT(VettBuoni,:);
% Ne=length(VettBuoni);
