求解微分方程组，用matlab或者mathmatica都可以。