qtzt.m

来自「四面体有限元程序,单元刚度的计算及总刚度矩征集成」· M 代码 · 共 8 行

subplot(1,2,1)
sphere(40)
axis equal
subplot(1,2,2)
t=0:0.1*pi:2*pi;
[X,Y,Z]=cylinder(4*cos(t));
mesh(X,Y,Z)
axis square