matlab开发-ROCK4。刚性方程的显式求解