ddqxnfun.m

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

function ydot=ddqxnfun(t,y,flag,b,p,m)
ydot=[y(2)
   -b/m*y(2)*(y(2).^2+y(4).^2).^(p/2);
   y(4);
   -9.8-b/m*y(4)*(y(2).^2+y(4).^2).^(p/2)];