function [x,y]=BilinearQuadElementcoordinate(NEI,NEJ,LENGTH,BREADTH)
%   給每一個結點配置Ｘ、Ｙ方向的坐標，按橫向分布
for J=1:NEJ+1
    for I=1:NEI+1
        X((J-1)*(NEI+1)+I)=(I-1)*LENGTH/NEI;
        Y((J-1)*(NEI+1)+I)=(J-1)*BREADTH/NEJ;
    end
end
x=X;
y=Y;