我正在编写代码(带有R)来解决带for循环和函数的逆矩阵。 编写for-loop子句存在问题。
我正在用我的函数求解带有协因子的逆矩阵。
##This is my example matrix
A <- matrix(data = c(1, 3, -4, 5, -2, 1, -1, 3, 3, -2, 2, -1, 4,-1, -3, -1),
byrow = TRUE,
nrow = 4, n col = 4)
##This is co factor function
library(functional)
cof = function(A, i, j)
{ stop if not(n col(A) == n row(A))
det(A[-i,-j])*(-1)^(i+j)
}
## this is the problem-writing for loop
n <- nrow(A)
f <- c(n)
f[1] <- 1
for (i in 1:n) {
for (j in 1:n) {
f[i,j] <- cof(A, i, j)
}
}
f
我如何分配 我的辅因子函数(A,i,j)到a(i,j) 从11到44 带有for循环子句。
请帮助