function displayNorms(srf,rot1,rot2) dims=size(srf) xDim=dims(1);yDim=dims(2);zDim=dims(3); % Display isosurface of torus and customize it hiso=patch(isosurface(srf,0.0),'FaceColor',[1,.75,.65],'EdgeColor','none'); view(rot1,rot2); lightangle(45,30); set(gcf,'Renderer','zbuffer'); lighting phong %Extract the triangles and vertices [f vert]=reducepatch(isosurface(srf,0.0)); [X,Y,Z]=meshgrid(1:xDim,1:yDim,1:zDim); % Display the normal vectors [Dxc,Dyc,Dzc]=calculateNormals(srf); coneplot(X,Y,Z,Dxc,Dyc,Dzc,vert(:,1),vert(:,2),vert(:,3),'quiver');