多个ode方程如何写入

newton_grc

小弟正在通过模型建立在重力场情况下，理想气体的运动情况 ；把粒子看成经典粒子，符合力学定律，模拟粒子的位置、速度随时间的变化。
现在有一个问题，如何能够一次写入数量巨大的ode方程阿？

newton_grc

就是这样的方程
ydot=[y(2); -1;y(4);0;
          y(6);-1;y(8);0
          ……                 
          ];