%Plot the original STL mesh:
function [xco,yco,zco]=Show_STL(filename,transp,colorfacet,colorline)



[stlcoords] = READ_stlascii(filename);
zco = squeeze( stlcoords(:,1,:) )';
xco = squeeze( stlcoords(:,2,:) )';
yco = squeeze( stlcoords(:,3,:) )';
% xco=xco+10;
hpat = patch(zco,xco,yco,colorfacet,'EdgeColor',colorline,'EdgeAlpha',0.1);
xlabel('x');
ylabel('y');
zlabel('z');
alpha(hpat,transp);
axis equal

end