function CheckFaceNormals(TR)
%CheckFaceNormals: This function checks normals of triangles.
%   
%   Input: TR: Triangulation matrix



figure, trisurf(TR, 'FaceColor', 'b')

IC = incenter(TR)
FC = faceNormal(TR)

hold on  
quiver3(IC(:,1),IC(:,2),IC(:,3), ...
     FC(:,1),FC(:,2),FC(:,3),0.5,'color','r');
end