clear?all
%i點坐標(biāo)為（xi，yi）一樣的方法假設(shè)j、m點，A為面積這里為1
syms?xi?yi?xj?yj?xm?ym?u
b=[yj-ymym-yiyi-yj]
c=[-xj+xm-xm+xi-xi+xj]
%假設(shè)常數(shù)Et/（4*（1-u^2））為h或者為1，為方便計算為1
for?（i=1:3）
????for?（j=1:3）
????????k1（ij）=b（i）*b（j）+（1-u）*c（i）*c（j）/2
????????k2（ij）=u*b（