Sommerfeld–Runge method