Rootfinding and interpolation with Runge-Kutta-Sarafyan methods