似乎可能大概也许可以把x,y,z合成一个三维向量，然后写成左边导数右边矩[…]阵乘向量，然后通过矩阵特征值求解？貌似ode有标准解法，可惜多年不用，都还给老师了。。。

re写成关于 x 的ode， 然后laplace transform, 因式分解成 sum of C/（A*S+B），再inverse laplace transform

。。。。。。。按照牛排的提示就可以了，写成矩阵形式。 如果觉得那样不熟悉，可以把三个变量转化成两个变量。 从方程组很容易看出三个方程相加为0，所以x+y+z=常数c，利用已知初始条件，x=-y-z+X，或者其他等式，消掉一个变量，就可以得到两个线性方程的式子了。另外，弱问，你的f，g那些是常系数还是函数？

。。。。。。。按照牛排的提示就可以了，写成矩阵形式。 如果觉得那样不熟悉，可以把三个变量转化成两个变量。 从方程组很容易看出三个方程相加为0，所以x+y+z=常数c，利用已知初始条件，x=-y-z+X，或者其他等式，消掉一个变量，就可以得到两个线性方程的式子了。另外，弱问，你的f，g那些是常系数还是函数？

【学术讨论】求解线性ODE 　 有没有大神能不能给个symbolic solution? 有没有人会用mathsmatica? 其实笔算应该可以但是我好久不用ODE已经不知道怎么做了。 感觉解出来的形式应该是 x(t) = A　 -　　 B.e^(-Ct)　 -　　 D.e^(-Et) 　 (more...)

Draft solution (not verified yet) 　 　 PS: if having time, will verify with python later because no access to Mathematica (more...)

Forgot the mention, x(0)=X is acting[…] as a scaling factor. In the draft solution it is taken as 1.

按照二次方程解法的惯例http://zh.wikipedia.org/wiki/%E4%BA%8C%E6%AC%A1%E6%96%B9%E7%A8%8B