用龙格库塔法求解耦合微分方程
 

各位大侠好：
小弟最近用龙格库塔法求解微分方程时，遇到下列情况：方程 dYi(t)/dt=5*Yi(t)+4*Xi(t)+[X1(t)+X2(t)+X3(t)+.....+X10(t)], 式中，Yi和Xi中的i为下标，i=1,2,3....10。Yi为待求量,其初值为Y0, Xi有一范围。这个公式求解的难点在于，dYi(t)/dt=5*Yi(t)+4*Xi(t)直接用龙格库塔法求初值问题， 好解，但加上后面那一坨后，就不知道怎么玩了。请各位大虾不吝赐教，万分感谢！