alpha <- 3
Y <- matrix(c(43,41,45,38,40,42,39,43), nrow = 8); Y
X <- matrix(c(1,0,1,0,1,1,0,1,0,1,0,1,0,0,1,0), nrow = 8); X
Z <- diag(8); Z
A <- diag(8); A

Xt <- t(X); Xt  #transpuesta de X
Zt <- t(Z); Zt

TI <- matrix(rbind(Xt%*%Y,Zt%*%Y)); TI
#combinación vertical, las matrices deben compartir el mismo ancho
#en una comb horizontal las matrices deben tener misma altura

MC <- rbind(cbind(Xt%*%X,Xt%*%Z),cbind(Zt%*%X,Zt%*%Z+solve(A)*alpha)); MC #combinación vertical de dos combinaciones horizontales

sol <- solve(MC,TI); sol #inversa de MC mult. por TI