Covariance() now produces the appropriately shaped matrix. Gamma
function has been roughly introduced, though needs to be double-checked
with fisher.nb