ddqxn.m

四面体有限元程序,单元刚度的计算及总刚度矩征集成

%%program  ddqxn.m
m=1;  b=[0,0.1,0.1];  p=[0,0,1];
figure
axis([0  9  0  4]);
hold on
for i=1:3
   [t,y]=ode45('ddqxnfun',[0:0.001:2],[0,5,0,8],[],b(i),p(i),m);
   comet(y(:,1),y(:,3))
end