Added a None return if the Hessian cannot be calculated due to numerical instabilities. Note that as it stands now, there is a problem with the hessian calculations of the time series
Solved a bug in calculating the Hessian, need to advance the timestep at the end