wn+1′=wn′+16⁢(k1+2⁢k2+2⁢k3+k4)+O⁡(h5),