function [ax] = plteles(meshele,meshnod,plotcolor)
% [ax] = plteles(meshele,meshnod,plotcolor)
% or
% l = 4000:5000;
% [ax] = plteles(meshele(l,:),meshnod,'r')
% ax = axis;
%
scaler = .8;
if nargin~=3
plotcolor = 'c';
end

lf = 1:length(meshele);
x4 = [ meshnod(meshele(lf,2),2) meshnod(meshele(lf,3),2)... 
       meshnod(meshele(lf,4),2) meshnod(meshele(lf,2),2)];
xx4 = scaler*(x4-[mean(x4')' mean(x4')' mean(x4')' mean(x4')'] )+...
  [mean(x4')' mean(x4')' mean(x4')' mean(x4')'];
y4 = [ meshnod(meshele(lf,2),3) meshnod(meshele(lf,3),3)... 
       meshnod(meshele(lf,4),3) meshnod(meshele(lf,2),3)];
yy4 = scaler*(y4-[mean(y4')' mean(y4')' mean(y4')' mean(y4')'] )+...
  [mean(y4')' mean(y4')' mean(y4')' mean(y4')'];


 plot(xx4',yy4',plotcolor)

ax = axis;

