本程序使用四阶龙格库塔法解常微分方程,使用者可根据自身需要修改要求解的方程