What doyoumean by "third order". I would be inclined to say that this algorithm uses the value of f at three difference points, (x, y), (x+ h/2, x+ h/2*k1), and (x+ h, y+ h(-k1+ 2k2)) and so is "third order" but you appear to have a different definition of "order".