Nota: Este artigo é sobre Sistema de equações lineares. Para outros significados, veja Sistema (desambiguação).

Esta página cita fontes, mas que não cobrem todo o conteúdo. Ajude a inserir referências (Encontre fontes: ABW  • CAPES  • Google (notícias • livros • acadêmico)). (Junho de 2021)

Em Matemática, um sistema de equações lineares (abreviadamente, sistema linear) é um conjunto finito de equações lineares aplicadas num mesmo conjunto, igualmente finito, de variáveis.^[1] Por exemplo, $${\displaystyle {\begin{alignedat}{7}3x&&\;+\;&&2y&&\;-\;&&z&&\;=\;&&1&\\2x&&\;-\;&&2y&&\;+\;&&4z&&\;=\;&&-2&\\-x&&\;+\;&&{\tfrac {1}{2}}y&&\;-\;&&z&&\;=\;&&0&\end{alignedat}}}$$

é um sistema de três equações com três variáveis (x, y e z). Uma solução para um sistema linear é uma atribuição de números às incógnitas que satisfazem simultaneamente todas as equações do sistema. Uma solução para o sistema acima é dada por $${\displaystyle {\begin{alignedat}{2}x&\,=\,&1\\y&\,=\,&-2\\z&\,=\,&-2\end{alignedat}}}$$

já que esses valores tornam válidas as três equações do sistema em questão. A palavra "sistema" indica que as equações devem ser consideradas em conjunto, e não de forma individual, isto é, cada uma das equações precisa ser satisfeita.

Em matemática, a teoria de sistemas lineares é a base e uma parte fundamental da álgebra linear, um tema que é usado na maior parte da matemática moderna. Deve-se observar que, em primeiro lugar, a equação linear é, necessariamente, uma equação polinomial. Em diversos ramos da matemática aplicada e ciências naturais, podemos encontrar vários usos de sistemas lineares. Exemplos são a física, a economia, a engenharia, a biologia, a geografia, a navegação, a aviação, a cartografia, a demografia e a astronomia.^[2]

Algoritmos computacionais são para encontrar soluções constituem uma parte importante da álgebra linear numérica, e desempenham um papel proeminente nas áreas de aplicação da álgebra linear. Tais métodos têm uma grande importância para obter soluções rápidas e acuradas.^[3] Pode-se muitas vezes aproximar um sistema de equações não-lineares por um sistema linear, uma técnica chamada de linearização e útil ao elaborar modelos matemáticos ou realizar simulações computacionais de um sistema mais complexo.

O sistema linear também pode ser conceituado como um sistema de equações do primeiro grau, ou seja, um sistema no qual as equações possuem apenas polinômios em que cada parcela tem apenas uma incógnita. Em outras palavras, num sistema linear, não há potência diferente de um ou zero e tampouco pode haver multiplicação entre incógnitas.

Muitas vezes, os coeficientes das equações são números reais ou complexos e as soluções são procuradas no mesmo conjunto de números, mas a teoria e os algoritmos aplicam os coeficientes e soluções em qualquer campo.

O tipo mais simples de sistema linear envolve duas equações e duas variáveis: $${\displaystyle {\begin{alignedat}{5}x&&\;+\;&&3y&&\;=\;&&1&\\2x&&\;+\;&&-y&&\;=\;&&-2&.\end{alignedat}}}$$ Um método para resolver tal sistema é do seguinte modo: em primeiro lugar, resolva a equação superior para x em termos de y: $${\displaystyle x=1-3y.}$$ Agora substitua essa expressão para x na equação inferior: $${\displaystyle 2\left(1-3y\right)-y=-2.}$$ Isto resulta numa única equação envolvendo apenas a variável ${\displaystyle y.}$ Resolvendo, obtemos ${\displaystyle y=4/7,}$ e voltando para a equação e substituindo y por seu valor, vem que ${\displaystyle x=-5/7.}$ Este método se generaliza para sistemas com variáveis adicionais (veja "eliminação de variáveis" abaixo, ou o artigo sobre álgebra elementar).

Um sistema geral de m equações lineares com n incógnitas pode ser escrito como:

${\displaystyle {\begin{alignedat}{7}a_{11}x_{1}&&\;+\;&&a_{12}x_{2}&&\;+\cdots +\;&&a_{1n}x_{n}&&\;=\;&&&b_{1}\\a_{21}x_{1}&&\;+\;&&a_{22}x_{2}&&\;+\cdots +\;&&a_{2n}x_{n}&&\;=\;&&&b_{2}\\\vdots \;\;\;&&&&\vdots \;\;\;&&&&\vdots \;\;\;&&&&&\;\vdots \\a_{m1}x_{1}&&\;+\;&&a_{m2}x_{2}&&\;+\cdots +\;&&a_{mn}x_{n}&&\;=\;&&&b_{m}.\\\end{alignedat}}}$

Aqui, ${\displaystyle x_{1},x_{2},\ldots ,x_{n}}$ são as incógnitas, ${\displaystyle a_{11},a_{12},\ldots ,a_{mn}}$ são os coeficientes do sistema e ${\displaystyle b_{1},b_{2},\ldots ,b_{m}}$ são os termos constantes.^[4]

Muitas vezes, os coeficientes e as incógnitas são números reais ou complexos, mas pode-se encontrar também números inteiros e racionais, já que são polinômios e elementos de uma estrutura algébrica abstrata.

Um ponto de vista extremamente útil é que cada incógnita é um peso para um vetor coluna em uma combinação linear.${\displaystyle x_{1}{\begin{bmatrix}a_{11}\\a_{21}\\\vdots \\a_{m1}\end{bmatrix}}+x_{2}{\begin{bmatrix}a_{12}\\a_{22}\\\vdots \\a_{m2}\end{bmatrix}}+\cdots +x_{n}{\begin{bmatrix}a_{1n}\\a_{2n}\\\vdots \\a_{mn}\end{bmatrix}}={\begin{bmatrix}b_{1}\\b_{2}\\\vdots \\b_{m}\end{bmatrix}}}$

Isso permite que seja exercida toda a linguagem e teoria dos espaços vetoriais (ou, mais geralmente, módulos). Por exemplo: o espaço vetorial gerado é o conjunto de todas as possíveis combinações lineares dos vetores sobre o lado esquerdo, e as equações têm uma solução apenas quando o vetor da mão direita se encontra nesse espaço vetorial gerado. Se cada vetor desse espaço vetorial gerado tem exatamente uma expressão como uma combinação linear dos vetores dados à esquerda, então qualquer solução é única. De qualquer maneira, o espaço vetorial gerado tem uma base de vetores linearmente independentes que garantem exatamente uma expressão; e o número de vetores nessa base (a sua dimensão) não pode ser maior do que m ou n, mas pode ser menor. Isto é importante porque, se tivermos m vetores independentes, a solução é garantida, independentemente do lado direito, o que não ocorre de outra forma.

A equação vetorial é equivalente a uma equação matricial da forma $${\displaystyle A\mathbf {x} =\mathbf {b} }$$ onde A é uma matriz m×n, x é um vetor coluna com n elementos e b é um vetor coluna com m elementos.

$${\displaystyle A={\begin{bmatrix}a_{11}&a_{12}&\cdots &a_{1n}\\a_{21}&a_{22}&\cdots &a_{2n}\\\vdots &\vdots &\ddots &\vdots \\a_{m1}&a_{m2}&\cdots &a_{mn}\end{bmatrix}},\quad \mathbf {x} ={\begin{bmatrix}x_{1}\\x_{2}\\\vdots \\x_{n}\end{bmatrix}},\quad \mathbf {b} ={\begin{bmatrix}b_{1}\\b_{2}\\\vdots \\b_{m}\end{bmatrix}}}$$

Agora, o número de vetores em uma base para o espaço vetorial gerado é expresso como o posto da matriz.

O método matricial pode ser usado para resolver sistemas de equações lineares em que o número de incógnitas é igual ao número de equações, ou seja, sistemas de equações linear com uma matriz quadrada dos coeficientes das incógnitas.

Outra condição para a aplicabilidade do método matricial é a não degeneração da matriz de coeficientes para incógnitas, ou seja, o determinante da matriz deve ser diferente de zero.

O sistema de equações lineares, quando as condições acima são satisfeitas, pode ser representado em forma de matriz, e então resolvido encontrando a matriz inversa da matriz do sistema.

A solução dos sistemas de equações lineares pelo método da matriz baseia-se na seguinte propriedade da matriz inversa: o produto da matriz inversa e da matriz original é igual à matriz identidade. ^[5]

Que seja necessário resolver um sistema de equações lineares:

a₁₁x₁ + a₁₂x₂ + a₁₃x₃ = b₁

a₂₁x₁ + a₂₂x₂ + a₂₃x₃ = b₂

a₃₁x₁ + a₃₂x₂ + a₃₃x₃ = b₃

Escrevemos este sistema de equações em forma de matriz:

| a11  a12  a13 |   | x1 |   | b1 | 
| a21  a22  a23 | x | x2 | = | b2 | 
| a31  a32  a33 |   | x3 |   | b3 |

Exemplo 1. Resolva o sistema de equações lineares pelo método da matriz:

2x₁ - x₂ = 1

3x₁ + 5x₂ = -1

A solução consiste nas seguintes etapas:

Etapa 1. Compomos as seguintes matrizes.

Matriz de coeficientes das incógnitas:

A = | 2  -1 |  
    | 3   5 |  

Matriz de incógnitas:

X = | x1 |   
    | x2 |  

Matriz de membros:

B = |  1 |  
    | -1 |  

Verificar se a matriz de coeficientes para incógnitas é diferente de zero, ou seja, aplicar o método da matriz:

| A | = 2 . 5 - (-1) . 3 = 13

O determinante desta matriz não é igual a zero, portanto, podemos aplicar o método da matriz.

Etapa 2. Encontre a matriz inversa de coeficientes para incógnitas:

A-1 = 1/13 . |  5  1 |  
             | -3  2 |  

Etapa 3. Encontre a matriz de valores desconhecidos:

X = 1/13 . | 5   1 |  . |  1 | =
           | -1  2 |    | -1 |

  = 1/13 . |  4 | = |  4/13 |
           | -5 |   | -5/13 |

Então, temos a solução:

x₁ = 4/13, x₂ = 5/13

Vamos checar:

2 . 4/13 - (-5/13) = 1

3 . 4/13 + 5 . (-5/13) = -1

Portanto, a resposta está correta.

O sistema linear é correlacionada à álgebra linear e o entendimento mais profundo dos sistemas é dependente do domínio desta matéria.^[6]

Sendo assim, é importante o entendimento dos espaços vetoriais, dos isomorfismos, das transformações lineares, da interpolação de Lagrange, da decomposição de um polinômio em fatores primos, de anéis comutativos, do teorema da decomposição primária, da forma de Jordan e das formas bilineares.

Um sistema linear, partindo da premissa de que tem resultado existente e determinado e não há dependência entre as equações, deve ter o mesmo número de equações e de incógnitas. O número de variáveis (incógnitas) também é chamado de quantidade de dimensões do problema. O número de dimensões está relacionado ao espaço vetorial. Por outro lado, os números que são subsumidos às incógnitas das equações podem ser de vários universos. Em geral, se resolvem sistemas para números reais, mas também existem sistemas para números complexos e ainda para outros tipos de números. Assim, para n dimensões no conjunto dos números reais, diz-se que se trabalha no conjunto ℝⁿ.

Para que o resultado de um sistema seja existente e determinado, não pode haver redundância, o que é chamado também dependência entre as matrizes que representam as equações.

A história dos sistemas de equações lineares começa no oriente. Em 1683, num trabalho do japonês Seki Kowa, surge a ideia de determinante^[7] (como polinômio que se associa a um quadrado de números).

O uso de determinantes no Ocidente começou dez anos depois num trabalho de Leibniz, ligado também a sistemas lineares.

A conhecida regra de Cramer é na verdade uma descoberta do escocês Colin Maclaurin (1698-1746), datando provavelmente de 1729, embora só publicada postumamente em 1748 no seu Treatise of algebra.

O suíço Gabriel Cramer (1704-1752) não aparece nesse episódio de maneira totalmente gratuita. Cramer também chegou à regra independentemente.

O francês Étienne Bézout (1730-1783), autor de textos matemáticos de sucesso em seu tempo, tratou do assunto, sendo complementado posteriormente por Laplace, em Pesquisas sobre o cálculo integral e o sistema do mundo.

O termo determinante, com o sentido atual, surgiu em 1812 num trabalho de Cauchy sobre o assunto. Neste artigo, apresentado à Academia de Ciências, sugeriu a notação que hoje é aceita como convenção.

Já o alemão Jacobi fez a leitura dessa teoria da forma como atualmente se estuda.

A solução de um sistema linear é a atribuição de valores às variáveis x₁, x₂, ..., x_n de modo a satisfazer ambas equações. O grupo de todas as soluções possíveis é chamado de conjunto solução.^[8]

Um sistema linear pode comportar-se em qualquer uma das três formas possíveis:^[9][10]

1. Sistema Possível Determinado – S.P.D.: É o sistema que possui apenas uma única solução possível, quando o determinante é diferente de 0 (D ≠ 0). Esse sistema pode ser chamado de sistema possível e suas equações de equações compatíveis.
2. Sistema Possível Indeterminado – S.P.I.: É o sistema que possui múltiplas soluções e o seu determinante é igual a zero (D = 0).
3. Sistema Impossível – SI: É o sistema que não admite uma solução, sendo designado por sistema impossível e suas equações como equações incompatíveis. Nesse caso, o determinante principal é igual a zero (D = 0) e um ou mais determinantes secundários são diferentes de zero (D ≠ 0). ^[11]

Existem inúmeros métodos de resolução de sistemas, como distintas propriedades matemáticas. Para aplicações computacionais complexas, veja algebra liner computacional. Alguns métodos papel e lapis existem.

O método de eliminação de Gauss, que também é conhecido como escalonamento é um método que consiste de manipular o sistema através de algumas operações elementares, que transforma a matriz estendida em uma matriz triangular, que após isso, a solução é obtida via substituição regressiva.

Seja Ax = b um sistema linear. O método de eliminação de Gauss para se encontrar a solução do sistema consiste nas seguintes etapas

• Etapa 1: Obter a matriz aumentada na forma ${\displaystyle {\begin{bmatrix}A|b\end{bmatrix}}}$ que representa o sistema de equações.
• Etapa 2:Transformar a matriz aumentada ${\displaystyle {\begin{bmatrix}A|b\end{bmatrix}}}$ em uma matriz aumentada na forma ${\displaystyle {\begin{bmatrix}{\bar {A}}|{\bar {b}}\end{bmatrix}}}$ onde ${\displaystyle {\bar {A}}}$ é uma matriz triangular superior.
• Etapa 3: Resolver o sistema linear ${\displaystyle {\begin{bmatrix}{\bar {A}}|{\bar {b}}\end{bmatrix}}}$ da Etapa 2 por substituição regressiva.

Considere o sistema linear de 3 equações abaixo: $${\displaystyle {\begin{alignedat}{7}a_{11}x_{1}&&\;+\;&&a_{12}x_{2}&&\;+\;&&a_{13}x_{3}&&\;=\;&&b_{1}&\qquad (L_{1})\\a_{21}x_{1}&&\;+\;&&a_{22}x_{2}&&\;+\;&&a_{23}x_{3}&&\;=\;&&b_{2}&\qquad (L_{2})\\a_{31}x_{1}&&\;+\;&&a_{32}x_{2}&&\;+\;&&a_{33}x_{3}&&\;=\;&&b_{3}&\qquad (L_{3})\end{alignedat}}}$$

A matriz aumentada A do sistema é:

${\displaystyle {\begin{bmatrix}A|b\end{bmatrix}}^{(0)}}$ = ${\displaystyle \left[{\begin{array}{ccc|c}a_{11}&a_{12}&a_{13}&b_{1}\\a_{21}&a_{22}&a_{23}&b_{2}\\a_{31}&a_{32}&a_{33}&b_{3}\end{array}}\right]}$

Deseja-se zerar todos os elementos da primeira coluna abaixo da diagonal principal. Assim, sendo ${\displaystyle a_{11}\neq 0,}$ define-se as constantes ${\displaystyle k=a_{21}/a_{11}}$ e ${\displaystyle w=a_{31}/a_{11}}$ e faz-se as seguintes operações lineares:

$${\displaystyle L_{2}^{(1)}\leftarrow L_{2}-k.L_{1}}$$

$${\displaystyle L_{3}^{(1)}\leftarrow L_{3}-w.L_{1}}$$

Obtendo-se:

$${\displaystyle {\begin{bmatrix}A|b\end{bmatrix}}^{(1)}}$$ = ${\displaystyle \left[{\begin{array}{ccc|c}a_{11}^{(1)}&a_{12}^{(1)}&a_{13}^{(1)}&b_{1}^{(1)}\\0&a_{22}^{(1)}&a_{23}^{(1)}&b_{2}^{(1)}\\0&a_{32}^{(1)}&a_{33}^{(1)}&b_{3}^{(1)}\end{array}}\right]}$

Agora, deve-se zerar todos os elementos da segunda coluna abaixo da diagonal principal. Sendo o pivô o elemento ${\displaystyle a_{22}^{(1)}}$ e a linha pivô a linha 2 de ${\displaystyle {\begin{bmatrix}A|b\end{bmatrix}}^{(1)},}$ supõe-se ${\displaystyle a_{22}^{(1)}\neq 0,}$ e define-se uma nova constante ${\displaystyle v=a_{32}^{(1)}/a_{22}^{(1)}.}$ Realizando a operação

$${\displaystyle L_{3}^{(2)}\leftarrow L_{3}^{(1)}-v.L_{2}^{(1)}}$$

obtém-se:

$${\displaystyle {\begin{bmatrix}A|b\end{bmatrix}}^{(2)}}$$ = ${\displaystyle \left[{\begin{array}{ccc|c}a_{11}^{(2)}&a_{12}^{(2)}&a_{13}^{(2)}&b_{1}^{(2)}\\0&a_{22}^{(2)}&a_{23}^{(2)}&b_{2}^{(2)}\\0&0&a_{33}^{(2)}&b_{3}^{(2)}\end{array}}\right]}$

• Note que ${\displaystyle {\begin{bmatrix}A|b\end{bmatrix}}^{(2)}}$ é uma matriz aumentada cuja matriz ${\displaystyle A}$ é uma matriz triangular superior.

Resolve-se o sistema ${\displaystyle {\begin{bmatrix}A|b\end{bmatrix}}^{(2)}.}$ Assim:

$${\displaystyle x_{3}=b_{3}^{(2)}/a_{33}^{(2)}}$$ ${\displaystyle ,\ a_{33}^{(2)}\neq 0}$ $${\displaystyle x_{2}=(b_{2}^{(2)}-(a_{23}^{(2)}x_{3}))/a_{22}^{(2)}}$$ $${\displaystyle x_{1}=(b_{1}^{(2)}-(a_{12}^{(2)}x_{2})-a_{13}^{(2)}x_{3})/a_{11}^{(2)}}$$

Assim, encontra-se a solução ${\displaystyle {\begin{Bmatrix}x_{1},\ x_{2},\ x_{3}\end{Bmatrix}}}$ do sistema ${\displaystyle {\begin{bmatrix}A|b\end{bmatrix}}^{(2)},}$ que é a mesma solução de ${\displaystyle {\begin{bmatrix}A|b\end{bmatrix}}.}$

• Observação: o método da Eliminação de Gauss só poderá ser usado para resolver sistemas lineares associados a matrizes escalonadas reduzidas com elementos das suas diagonais principais não-nulos, ou seja, ${\displaystyle a_{11}^{(1)},a_{22}^{(2)},a_{33}^{(3)},...,a_{nn}^{(n)}\neq 0.}$

O método da substituição consiste em isolar uma incógnita em qualquer uma das equações, obtendo igualdade com um polinômio. Então deve-se substituir essa mesma incógnita em outra das equações pelo polinômio ao qual ela foi igualada.

Consiste em compararmos as duas equações do sistema, após termos isolado a mesma variável (x ou y) nas duas equações. e as equações ficam mais detalhadas.

O método de adição, é uma solução de sistemas lineares. Porem é mais usado quando se tem duas equações e duas incógnitas, Basicamente deve-se somar as equações a fim de cancelarmos uma variável no sistema.

Exemplo:

      21x + y = 13
      29x - y = 12

Podemos organiza-la para facilitar, deixando-a assim:

     21x + 29x + y - y = 13 + 12Texto em subscrito

Resultando em:

     21x + 29x = 13 + 12
     50x = 25
     x = 50 / 25 => 2

Assim encontramos o resultado da variável X, porem não acabou, temos que substituir o x na segunda equação

     29x - y = 12
     29(2) - y = 12 
     58 - y = 12
     y = - 12 + 58
     y = 46

Sendo assim, o conjunto solução é:

     S ( x , y ) = (2 , 46)

Os métodos mais utilizados computacionalmente para resolver sistemas lineares envolvem fatorizações de matrizes. O mais conhecido, a eliminação de Gauss, origina a fatoração LU. Resolver o sistema Ax=b é equivalente a resolver os sistemas mais simples Ly=b e Ux=6.

A Regra de Cramer é uma fórmula explícita para a solução de um sistema de equações lineares, com cada variável dada por um quociente de dois determinantes. Por exemplo, a solução para o sistema

$${\displaystyle {\begin{alignedat}{7}x&\;+&\;3y&\;-&\;2z&\;=&\;5\\3x&\;+&\;5y&\;+&\;6z&\;=&\;7\\2x&\;+&\;4y&\;+&\;3z&\;=&\;8\end{alignedat}}}$$

é dada pela

$${\displaystyle x={\frac {\,\left|{\begin{matrix}5&3&-2\\7&5&6\\8&4&3\end{matrix}}\right|\,}{\,\left|{\begin{matrix}1&3&-2\\3&5&6\\2&4&3\end{matrix}}\right|\,}},\;\;\;\;y={\frac {\,\left|{\begin{matrix}1&5&-2\\3&7&6\\2&8&3\end{matrix}}\right|\,}{\,\left|{\begin{matrix}1&3&-2\\3&5&6\\2&4&3\end{matrix}}\right|\,}},\;\;\;\;z={\frac {\,\left|{\begin{matrix}1&3&5\\3&5&7\\2&4&8\end{matrix}}\right|\,}{\,\left|{\begin{matrix}1&3&-2\\3&5&6\\2&4&3\end{matrix}}\right|\,}}.}$$

Para cada variável, o denominador é a determinante da matriz de coeficientes, enquanto o numerador é o determinante de uma matriz na qual cada coluna foi substituída pelo vetor de termos constantes.

Embora a regra de Cramer seja importante teoricamente, tem pouco valor prático para grandes matrizes, uma vez que o cálculo de grandes determinantes é um processo computacionalmente custoso.

Além disso, a regra de Cramer tem pobres propriedades numéricas, tornando-a inadequada para resolver até mesmo pequenos sistemas de forma confiável, a menos que as operações forem executadas em aritmética racional com precisão ilimitada.

O método de Jacobi é obtido a partir do sistema linear:

a₁₁x₁ + a₁₂x₂ + ... + a_1nx_n = y₁

a₂₁x₁ + a₂₂x₂ + ... + a_2nx_n = y₂

...

a_n1x₁ + a_n2x₂ + ... + a_nnx_n = y_n

Agora, isolando x₁ da primeira linha, obtemos:

x₁^(k+1) = y₁ - (a₁₂x₂^(k) + ... + a_1nx_n^(k)) / a₁₁

Fazendo isso para os próximos x_i de cada equação linear, o resultado será algo como:

x₁^(k+1) = y₁ - (a₁₂x₂^(k) + ... + a_1nx_n^(k)) / a₁₁

x₂^(k+1) = y₂ - (a₂₁x₂^(k) + a₂₃x₃^(k) + ... + a_2nx_n^(k)) / a₂₂

...

x_n^(k+1) = y₂ - (a_n1x₁^(k) + ... + a_n,n-2x_n-2^(k) + a_n,n-1x_n-1^(k)) / a_nn

Uma forma mais abreviada de escrever a equação seria:

x₍₁₎ = aproximação inicial

x_i^(k+1) = (y_i - Σ_{j=1j ≠ i}ⁿ a_ijx_j^(k)) / a_ii

Logo, resolvendo o sistema linear abaixo via Método de Jacobi:

10x + y = 23

x + 8y = 26

Iniciando com x⁽¹⁾ = y⁽¹⁾ = 0.

x^(k+1) = 23 - y^k / 10

y^(k+1) = 26 - x^k / 8

x⁽²⁾ = 23 - y⁽¹⁾ / 10 = 2,3

y⁽²⁾ = 26 - x⁽¹⁾ / 8 = 3,25

x⁽³⁾ = 23 - y⁽²⁾ / 10 = 1,975

y⁽³⁾ = 26 - x ⁽²⁾ / 8 = 2,9625^[12]

Referências

• Álgebra Linear - Um Livro Colaborativo, mantido pelo Instituto de Matemática e Estatística da Universidade Federal do Rio Grande do Sul.
• Solução de sistemas linear em Cálculo Numérico - Um Livro Colaborativo, mantido pelo Instituto de Matemática e Estatística da Universidade Federal do Rio Grande do Sul