Internal itp functions
itp_cpp(f, pars, a, b, ya, yb, epsilon, k1, k2, for_rk, inc) itp_r(f, ..., a, b, ya, yb, epsilon, k1, k2, for_rk, inc)
These functions are not intended to be called by the user.