stvmat function

Variance matrix for space-time clustering