vcov.bellreg function

Variance-covariance matrix for a bellreg model