Condition singularity warning on unit root check

When singularity is detected, try to compute eigenvalues to see if unit root is present and issue warning according to this check
Allow for check of Hessian at order>2 in model_diagnostics.m
Merge pull request #768 from JohannesPfeifer/model_diagnostics

Improvements to model diagnostics