In [1]:
# rm(list=ls()) 
options(OutDec = ",") 
#==============================================================================
# Funcao para a solucao de sistemas lineares dado um sistema triangular 
# inferior
#==============================================================================
sistema.triang <- function(L,b){
    n    <- length(b)
    x    <- rep(NA,n)
    x[1] <- b[1]/L[1,1]
    for(i in 2:n)
        x[i] <- (b[i]-sum(L[i,1:(i-1)]*x[1:(i-1)]))/L[i,i]
    return(x)
}
In [2]:
#==============================================================================
# Exemplo 
#==============================================================================
L <- matrix(c(1,0,0,2,1,0,2,4,1),3,3,byrow=T)
b <- c(3,4,5)
print("===============================")
print(L)
print("===============================")
print(b)

[1] "==============================="
     [,1] [,2] [,3]
[1,]    1    0    0
[2,]    2    1    0
[3,]    2    4    1
[1] "==============================="
[1] 3 4 5
In [3]:
#==============================================================================
sistema.triang(L,b) ## 3 -2 7
solve(L,b)          ## 3 -2 7
  1. 3
  2. -2
  3. 7
  1. 3
  2. -2
  3. 7
In [4]:
# rm(list=ls()) 
#==============================================================================
# Fim
#==============================================================================