R/riccati_root.R
riccati_root.Rd
Compute the square root of a symmetric matrix.
riccati_root(x)
squared and symmetric matrix.
cv <- matrix(c(1, 0.3, 0.3, 1), nrow = 2, byrow = TRUE) riccati_root(cv) #> [,1] [,2] #> [1,] 0.9884177 0.1517577 #> [2,] 0.1517577 0.9884177