Estimation of Solving Time Left using ODE15s()

I'm using [t,sol] = ode15s( ... ) in order to solve a system of stiff first order ODEs. This works quite well. Used within a simulation environment it would be of great interest to estimate the time remaining.

Someone got any ideas?