请教各位大侠解达芬方程的问题

yhp1980

  用matlab解单自由度的达芬方程，如果想在每等间隔（0.02s）加人工加速度，并想在每等间隔得出一个解，那么，在程序上应该如何体现呢？<BR>    另外，有没有高手编过多自由度达芬体系的求解程序，望指教一二。