function ReadSTLTest()
%读取和显示STL文件
VerMat=ReadSTLFile();
[pt,pgon]=ConstrucTopology(VerMat);
showpgon(pt,pgon);
function showpgon(pt,pgon)
%显示三角网格
view(3);
pat=zeros(1,length(pgon));
hold on
for i=1:length(pgon)
ptpgoni=pt(pgon{i},:);
pat(i)=patch(ptpgoni(:,1) , ptpgoni(:,2) , ptpgoni(:,3) ,[1.000 0.7000 0.3000], 'facelighting','flat','edgecolor',...
[0 1.00 1.00], 'edgelighting','flat');
end
hold off
axis equal