function N=ShapeFun(s,t)             % 4级子程序，矩形单元形函数矩阵
%% 四边形八结点等参单元形函数矩阵 
% 角点
N1=1/4*(1-s)*(1+t)*(-s+t-1); 
N2=1/4*(1-s)*(1-t)*(-s-t-1); 
N3=1/4*(1+s)*(1-t)*(s-t-1); 
N4=1/4*(1+s)*(1+t)*(s+t-1); 

% 边中点 
N5=1/2*(1-t^2)*(1-s); 
N7=1/2*(1-t^2)*(1+s); 
N6=1/2*(1-s^2)*(1-t); 
N8=1/2*(1-s^2)*(1+t); 
N=[N1 0 N2 0 N3 0 N4 0 N5 0 N6 0 N7 0 N8 0;
0 N1 0 N2 0 N3 0 N4 0 N5 0 N6 0 N7 0 N8];