I need to implement a ode function inside other because the solver methods are different and I want transfer values between them.

The problem consist in the implementation of a fractional controller to a integer order system. To solve the system I can use the ode45 but to implement the fractional differential equation of the controller I'm using the ode of Garrappa based in Predictor-Corrector of Diethelm.