model.matrix.systemfit function

Construct Design Matrices for Systems of Equations