Now let's suppose I would integrate my problem in auxiliary variable, let' say E. It is only theoretically. (e.g. E = t^2 or E = sin(t),...). I am able to expres t0 e in tf in terms of E. How can I solve my ODE now? It is sufficient to multiply y' with dt/dE and use my Runge Kutta?