<- matrix(c(-1, 0, 0, 0, 2, -1, 0, 0, 0, 0, -3, 1, 0, -1, 1, 0, 0, 1, -1, 0, 0,
A 1, 2, 0, 0, 0, 0, -1/2), nrow = 7, byrow = TRUE) # matriz A de estequiometria de reações
rownames(A) <- list("glicose", "piruvato", "H2O", "CoA", "acetil CoA", "CO2", "O2") # etiquetas dos reagentes
# matriz A de reações A
[,1] [,2] [,3] [,4]
glicose -1 0 0 0.0
piruvato 2 -1 0 0.0
H2O 0 0 -3 1.0
CoA 0 -1 1 0.0
acetil CoA 0 1 -1 0.0
CO2 0 1 2 0.0
O2 0 0 0 -0.5
<- c(-1, 0, 6, 0, 0, 6, -6) # vetor y de rotas
y
solve(t(A) %*% A) %*% t(A) %*% y
[,1]
[1,] 1
[2,] 2
[3,] 2
[4,] 12