function BilinearQuadElementShowmesh(X,Y,ME,NE,ND)
%    LinearTriangleElement的作用是畫出有限元網格、顯示各個結點和單元的編號。
%    在結點上標出的是結點號，為了防止單元和結點的編號混淆，我把單元的編號放在單元中心位置，且用方括號表示，如[2]，表示第二個單元，
plot(0,0);
hold on;
for I=1:NE
       for J=1:ND
        MEE=ME(J,I);
        x(J)=X(MEE);
        y(J)=Y(MEE);
        text(x(J),y(J),int2str(MEE));
    end
    x(ND+1)=x(1);
    y(ND+1)=y(1);
    plot(x,y,'r');
    s=0;
    t=0;
    for m=1:ND
       s=s+x(m);
       t=t+y(m);
   end
   dx=x(1)-x(2);
   dy=y(1)-y(3);
   xx=s/ND;
   yy=t/ND;
   text(xx+dx/6,yy-dy/8,['[',int2str(I),']']);
end