Extract Jacobian from ODE-Solver for each timestep

I have a system of ODE'S:

dy(t)/dt=f(y(t)) - y and f are vectors (14,1)

When I solve this (stiff) system via one of the solvers, for instance ode15s I get a solution vector Y(nx14) corresponding to a time vector T(nx1).

Now I know that the solver calculates the Jacobian of f(y(t)) during each step to see whether the tolerances have been met.

My question is: Is it possible to extract this Jacobian calculated by numjac as additional data for each time step while the solver is running?

If not can I somehow get the Jacobian later on?