急问：龙格库塔法解耦合微分方程

caicai43

各位大侠好：
    小弟最近用龙格库塔法求解微分方程时，遇到下列情况：方程 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)直接用龙格库塔法求初值问题， 好解，但加上后面那一坨 Xi 的各项之和后，就不知道怎么玩了。请各位大虾不吝赐教，万分感谢！