Output Jacobian from of ODE solver's (ODE15s) first iteration to feed into subsequent iterations to increase efficiency

My system requires that I iterate through calls to ODE15s. Ideally, I would supply a Jacobian so that the ODE solver doesn't have to recalculate it each time it starts, but my system is too complex so I'd rather avoid this approach. However, if I could somehow output the solvers approximation of the Jacobian after the first call to ODE15s, then feed it back into the subsequent calls to ODE15s, I might be able to greatly increase the computational efficiency.

Does anybody know how I can output a Jacobian from ODE15s?