经典R-K方法（四阶） y+1=y,+h(飞+2k2+2k3+k4)/6, k1=f(X,y,), k2=f(x,+h/2,y,+hk/2), k3=f(x,+h/2,y+hk2/2), k4=f(x,+h,y:+hk3): 局部截断误 差为O(h)经典R-K方法（四阶） 1 1 2 3 4 1 2 1 3 2 4 3 ( 2 2 ) / 6, ( , ), ( / 2, / 2), ( / 2, / 2), ( , ). i i i i i i i i i i y y h k k k k k f x y k f x h y hk k f x h y hk k f x h y hk  + = + + + +  =    = + +  = + +   = + + 局部截断误 差为O(h5 )