[有限元] 使用 Matlab 证明十结点三角形单元的完备性和协调性

>> 
syms xi xj xm yi yj ym x y real;
ai=xj*ym-xm*yj;
aj=xm*yi-xi*ym;
am=xi*yj-xj*yi;
bi=yj-ym;
bj=ym-yi;
bm=yi-yj;
ci=-(xj-xm);
cj=-(xm-xi);
cm=-(xi-xj);
delta=1/2*det([1 xi yi;1 xj yj;1 xm ym]);
Ni=1/(2*delta<