function []=PlotTriMesh(NodesPts,triangle);
triangle=triangle([1,3,2],:);
Nx=NodesPts(1,:);
Ny=NodesPts(2,:);
Nz=NodesPts(3,:);
ntri=size(triangle,2);
X=matrix(Nx(triangle),3,ntri);
Y=matrix(Ny(triangle),3,ntri);
Z=matrix(Nz(triangle),3,ntri);
plot3d(X,Y,Z);
set(gca(),'isoview','on');
set(gca(),'box','off');
endfunction;




