Linguagem r, criar um programa para resolver um sistema

Estou tentando criar um programa, usando a lógica da programação no R, para resolver um sistema de equações qualquer (vindo do excel). todavia, creio que a lógica esteja parcialmente errada: PRECISO MUITO DE AJUDA PARA FAZER RODAR sist_equacao = as.matrix(sist_equacao) class(sist_equacao) n = nrow(sist_equacao) c = ncol(sist_equacao) sist_equacaoD1 = matrix(0, nrow = n, ncol = c) i = 1 while (i <= n) { for (j in 1:c) { sist_equacaoD1[i,j] = sist_equacao[i,j]/sist_equacao[i,i] } i == i + 1 for (j in 1:c) { sist_equacaoD1[i,j] = sist_equacaoD1[i,j] - (sist_equacaoD1[i,j] * sist_equacaoD1[i-1,j]) } }