西交大计算方法上机可用。RK-4 该算法用Runge-Kutta法解初值问题。a,b—求解区间的左右端点,y0—初 值,步长选用h。依次由公式计算出K1,K2,K3,K4的值并计算出出每步的y(n)。 同时每算依次就更新依次t(n)的取值,最终输出ti和y