ys.m

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

[x,y,z] = meshgrid(-2:.2:2,-2:.2:2,-2:.2:2);
 v = x .* exp(-x.^2 - y.^2 - z.^2);          
 slice(v,[5 15],15,10) 
 %mesh(x,y,z,v);
 axis([0 21 0 21 0 21]);                     
 hold on                
 colorbar('horiz')
 colorbar('vert')
 view([-25 65])