四阶变步长龙格库塔matlab程序(齿轮系统动力学),适用于机械动力学