function [Nnode,Nele]=MeshPlot(Node,Element)
%plot mesh
figure
hold on;
patch('Faces',Element,'Vertices',Node,'FaceColor','red')

axis equal
title('Undeformed Geometry');

num=0;
for i=1:numel(Node(:,1))
    num=num+1;
    plot(Node(i,1),Node(i,2),'bo')
    text(Node(i,1),Node(i,2),0,num2str(num));
end
Nnode=size(Node,1);
Nele=size(Element,1);