##############################################################################
#                          FÍSICA COMPUTACIONAL II                           #
#                                   por                                      #
#                         Francisco  Carlos  Lavarda                         #
##############################################################################


Tutorial de MAXIMA - Sistemas de Equações Lineares
==================================================

Objetivo: Resolução de sistemas de equações lineares


LINSOLVE
--------
** Digite a primeira equação chamada eq1:
eq1: 94*x1+87*x2+37*x3=17108$

** Digite a segunda equação:
eq2: 19*x1+0*x2+13*x3=2945$

** Digite a terceira equação:
eq3: -96*x1+3*x2-42*x3=-12453$

** Use a função LINSOLVE para resolver o sistema linear [eq1,eq2,eq3] nas
variáveis [X1,X2,X3]: 
linsolve([eq1,eq2,eq3],[x1,x2,x3]);


ALGSYS
------
** Também é possível resolver o mesmo sistema de equações com ALGSYS:
algsys([eq1,eq2,eq3],[x1,x2,x3]);
A diferença é que ALGSYS também pode resolver sistemas de equações
não-lineares.


LINSOLVE_BY_LU (DVJ)
--------------
Esta função ainda não é documentada e não está disponível em versões mais
antigas do Maxima. Alguma documentação pode ser encontrada na WWW via
buscas.

** Vamos usar o mesmo exemplo 3x3 acima. Primeiramente defina a matriz A dos
coeficientes:
a:matrix([94,87,37],[19,0,13],[-96,3,-42]);

** Defina o vetor B das constantes:
b:[17108,2945,-12453];

** use o comando:
linsolve_by_lu(a,b);

Na resposta vem junto um "false": trata-se de uma variável que indica que
não foi utilizado algum método de cálculo numérico para resolver o problema.
(RMR)