Y +1/6 (K +2k +2k +K) H Z +1/6 (L +2l +2l +L) H F (X, Y, Z) F (X +1/2h, y +1/2k, Z +1/2l) F (X +1/2h, y +1/2k, Z +1/2l) F (X +H, y +K, Z +L)
Y +1/6 (K +2k +2k +K) H Z +1/6 (L +2l +2l +L) H F (X, Y, Z) F (X +1/2h, y +1/2k, Z +1/2l) F (X +1/2h, y +1/2k, Z +1/2l) F (X +H, y +K, Z +L)
Y +1/6 (K +2k +2k +K) H Z +1/6 (L +2l +2l +L) H F (X, Y, Z) F (X +1/2h, y +1/2k, Z +1/2l) F (X +1/2h, y +1/2k, Z +1/2l) F (X +H, y +K, Z +L)
5
g(x,y,z) = dz/dx = -8y-0.6z
y(0) = 4 y'(0)= 0 z(0)
yi+1=yi+1/6(k0+2k1+2k2+k3)h
zi+1=zi+1/6(l0+2l1+2l2+l3)h
k0=f(xi,yi,zi)
k1=f(xi+1/2h,yi+1/2k0h,zi+1/2l0h)
k2=f(xi+1/2h,yi+1/2k1h,zi+1/2l1h)
k3=f(xi+h,yi+k2h,zi+l2h)
l0=g(xi,yi,zi)
l1=g(xi+1/2h,yi+1/2k0h,zi+1/2l0h)
l2=g(xi+1/2h,yi+1/2k1h,zi+1/2l1h)
l3=g(xi+h,yi+k2h,zi+l2h)