Метод Якоби в линейных системах

Question

Здравствуйте. Я бы хотел использовать метод jacobi в своей собственной функции, чтобы получить лучшие предположения для b1, b2 и b3 на отображаемом изображении.

Пока что я попробовал эту функцию, которая давала мне ответы [2, 4, 1.625].Правильные ответы, однако: [-1, 3, 2].Что я делаю не так?

jacobi <- function(x, y){
  b <- matrix(0,3)
  n <- nrow(b)
  som <- 0
    for (i in 1:n) {
      for (j in 1:n) {
        if (i != j){
          som = x[i,j]*b[j]+som
        }
      }
     u[i] <- (y[i]-som)/x[i,i]
     b<-u
    }
  print(b)
}

D <- matrix(c(4,3,-2,2,-5,3,3,2,8),3,3)
E <- matrix(c(8,-14,27))
jacobi(D,E)