pendg.m

求解二阶微分代数方程组源代码 这是求解二阶微分代数方程组DAE

function pendg(t,w,wd)
% Used by pendrun.m
% Tony Roberts, 24 July 1998, [e-mail]
x=w(1); y=w(2); u=w(3); v=w(4); s=w(5);
plot([0 -x],[0 -y],-x,-y,'o',-x-[0 u]/8,-y-[0 v]/8)
axis([-1 1 -1 1])
axis('equal')
axis('off')
drawnow