15 lines
570 B
Mathematica
15 lines
570 B
Mathematica
|
function Nst=DiffShapeFun(s,t)
|
|||
|
|
% <EFBFBD><EFBFBD><EFBFBD><EFBFBD> s <EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
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;
|
|||
|
|
|
|||
|
|
% <EFBFBD><EFBFBD><EFBFBD><EFBFBD> t <EFBFBD><EFBFBD><EFBFBD><EFBFBD>
|
|||
|
|
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);
|
|||
|
|
|
|||
|
|
Nst=[N1s N1t;N2s N2t;N3s N3t;N4s N4t];
|
|||
|
|
end
|