set.seed(200100)
L <- matrix(rnorm(100), nrow = 10)
d <- rnorm(10)
equilibrium_output(L, d)
#> [,1]
#> [1,] -0.55849992
#> [2,] -0.28172718
#> [3,] 0.09158624
#> [4,] 0.29651982
#> [5,] -0.98302870
#> [6,] 0.18353881
#> [7,] 0.38595864
#> [8,] -1.01658063
#> [9,] 0.71849251
#> [10,] -1.01093063