function [ax] = ploteles(meshele,meshnod,scaler)
% [ax] = ploteles(meshele,meshnod)
% or
% [ax] = ploteles(meshele,meshnod,scaler)
%   scaler defaults to .8,,This is the shrinkage to create 
%    inner triangles
% ax = axis;
%
if nargin~=3
scaler = .8;
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','c',meshnod(:,2),meshnod(:,3),'ro')

ax = axis;

