wave1.m

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

t=0:pi/20:4*pi;   
x=0:0.1*pi:4*pi;
m=moviein(40)
for i=1:40
   y=sin(x-t(i));
   plot(x,y);
   axis([0,12,-1,1]);
   m(:,i)=getframe;
end
movie(m,3,30)