% plottest.m     Script, not a function
% plots lots of versions of mesh to see if there is an error
figure(1)
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

figure(2)
 [ax] = ploteles(meshele,meshnod,.8);

hold on
for i = 1:length(bcs)
lab = ['-',int2str(i),':',int2str(bcs(i,1))];
text(meshnod(bcs(i,1),2),meshnod(bcs(i,1),3),lab)
end


