This function converts the (q x 1) vector of constants and (q x qp) matrix of autoregressive coefficients into (qp x qp) matrix belonging to the companion form
companionMat(phi, p, q)
Arguments
phi: matrix of dimension (q x qp) containing autoregressive coefficients.