vcov()/SE calculation straight from Hessian

when possible (add use.hessian argument to vcov and summary)
closes #47