cjdby1:用matlab解常微分方程 不允许用dsolve 要求写出具体算法