Zmatrix function

Computes the Z-matrix