compute_laurent function

Laurent Expansion of Inverse of Linear Matrix Function

This function computes the first two terms of the Laurent expansion of the inverse of a linear matrix function.

compute_laurent(V, zapsmall = TRUE)

Arguments

• V: for some integer m >= 1, an array of dimension (m, m, 2), where V[, , 1] is the intercept and V[, , 2] is the slope of the linear matrix function.
• zapsmall: logical: should zapsmall be called on the result? Default TRUE.

Returns

array of dimension (m, m, 2), where W[, , 1] corresponds to the exponent -1, and W[, , 2] corresponds to the exponent 0.