求線性方程式解

GNU-R:

P <- array(0, dim=c(3,3)) P[1,1] <- 1 P[1,2] <- -1 P[1,3] <- 1 P[2,1] <- 3 P[2,2] <- 2 P[2,3] <- 1 P[3,1] <- 4 P[3,2] <- 2 P[3,3] <- 2 print (P) V <- array(0, dim=c(3,1)) V[1,1] <- 4 V[2,1] <- 2 V[3,1] <- 8 print (V) print(solve(P, V))