#--------------------------------------------
# Soluzione dell'esempio guida 
# (par. 2.9 circuiti_in_corrente_continua.pdf)
# con metodi di algebra lineare in R
#
# GdA 11/11/2021
#----------------------------------------------

# matrice dei coefficienti
A <- rbind(c(-45, 10, 0),
           c(45, 0, 10),
           c(1, 1, -1))
print(A)
cat("\n")

# vettore dei 'termini noti' 
f <- c(-2.1, 1.9, 0)
print(f)
cat("\n")

# matrice inversa
iA = solve(A)
print(iA)

# soluzione
cat("\n I1, I2, I3:\n")
I = as.vector( iA %*% f )
print(I)