bug.m

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

x=0.6;
u=2.6:0.001:4;

for j=1:150
   x=u.*(x-x.^2);
end

for i=1:100
   x=u.*(x-x.^2);
   plot(u,x,'r.','marker','.','markersize',1)
   hold on;
end