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