% plotmesh.m     Script, not a function
% just does the basics for meshnod
if length(meshnod)>1000 
plot(meshnod(:,2),meshnod(:,3),'.y',...
     meshnod(bndele(:,2),2),meshnod(bndele(:,2),3),'.r')
else
plot(meshnod(:,2),meshnod(:,3),'oy',...
    meshnod(bndele(:,2),2),meshnod(bndele(:,2),3),'*r')
end
axis('equal')
if exist('filename')
title([filename,'.nod'])
end

