matlab 二阶微分方程组解function Dy=dfun(t,y)D2y(1,1)=(0.0176*y(2)-96