Euler中点法、4阶龙格库塔法和对角隐式辛积分R-K方法，Matlab单步求解。