function Nst=DiffShapeFun(s,t)   
    
    N1s=(s/4 - 1/4)*(t + 1) + ((t + 1)*(s - t + 1))/4; 
    N2s=-(s/4 - 1/4)*(t - 1) - ((t - 1)*(s + t + 1))/4; 
    N3s=((t - 1)*(t - s + 1))/4 - (s/4 + 1/4)*(t - 1); 
    N4s=(s/4 + 1/4)*(t + 1) + ((t + 1)*(s + t - 1))/4; 
    N5s=t^2/2 - 1/2; 
    N6s=s*(t-1); 
    N7s=1/2 - t^2/2; 
    N8s=-s*(t + 1); 
    % ¹طسع t اَµ¼ 
    N1t=(s/4 - 1/4)*(s - t + 1) - (s/4 - 1/4)*(t + 1); 
    N2t=-(s/4 - 1/4)*(t - 1) - (s/4 - 1/4)*(s + t + 1); 
    N3t=(s/4 + 1/4)*(t - s + 1) + (s/4 + 1/4)*(t - 1); 
    N4t=(s/4 + 1/4)*(t + 1) + (s/4 + 1/4)*(s + t - 1); 
    N5t=t*(s - 1); N6t=s^2/2 - 1/2; 
    N7t=-t*(s + 1); N8t=1/2 - s^2/2; 
    Nst=[N1s N1t;N2s N2t;N3s N3t;N4s N4t;N5s N5t;N6s N6t;N7s N7t;N8s N8t;];
end