    El subprograma "ECUACION":
-Resuelve sistemas de N ecuaciones con N incgnitas
-Calcula determinantes
-Calcula inversas

    Se recomienda poner el programa en el rea P2 y llamarlo con
    GOSUB #2
    siempre que algn programa necesite resolver un sistema de ecuaciones
    (p.e., el programa ARTICULA lo hace, por lo que necesita tener el
    subprograma en el rea de programas P2)

    Los datos se pasan en una matriz de MxN elementos, [DIM A(M,N)]
donde N es el nmero de filas y M el de columnas; tanto M como N
empiezan en 0 [la matriz va de A(0,0) hasta A(M,N)]

    Si M=N, entonces se calcula el determinante, que se devuelve en la
variable K; el contenido de A() tras ejecutarse el programa no tienen
significado.

    Si M>N, entonces se resuelve el (o los) sistemas de ecuaciones
asociados, por el mtodo de Gauss-Jordan; el determinante de la
matriz de coeficientes se devuelve en K y los resultados se devuelven
en A(N+1,i) hasta A(M,i). Se recuerda que para calcular una inversa
basta con resolver N sistemas de ecuaciones en los que la matriz de
coeficientes es A y los trminos independientes son todos cero menos
uno que vale la unidad; para entenderlo, ver los ejemplos:

Ejemplo 1:
Hallar el determinante de
         
 1  2  3 
 4  5  6 
 7  8  9 
         
entonces M=2,N=2 [as, A vara de (0,0) a (2,2)]
A(0,0)=1,A(0,1)=2,A(0,2)=3
A(1,0)=4,A(1,1)=5,A(1,2)=6
A(2,0)=7,A(2,1)=8,A(2,2)=9
As, tras ejecutar el programa, K=0 (determinante nulo)

Ejemplo 2:
Resolver el sistema
3x+5y-6z=2
4x   + z=0
    y-4z=5
entonces M=3,N=2
A(0,0)= 3,A(1,0)= 5,A(2,0)=-6,A(3,0)= 2
A(0,1)= 4,A(1,1)= 0,A(2,1)= 1,A(3,1)= 0
A(0,2)= 0,A(1,2)= 1,A(2,2)=-4,A(3,2)= 5
tras ejecutar el programa, K=53 (determinante matriz de coeficientes)
y los resultados:
x=A(3,0)= 0.433962...
y=A(3,1)=-1.943396...
z=A(3,2)=-1.735849...

Ejemplo 3:
Calcular la inversa de
                                                      
   1  0  2                          1  0  2  1  0  0 
A= 0  1  0     orlando la  matriz:  0  1  0  0  1  0 
   2  0  1                          2  0  1  0  0  1 
                                                      
entonces M=5 (3+3-1) y N=2
A(0,0)=1,A(1,0)=0,A(2,0)=2,A(3,0)=1,A(4,0)=0,A(5,0)=0
A(0,1)=0,A(1,1)=1,A(2,1)=0,A(3,1)=0,A(4,1)=1,A(5,1)=0
A(0,2)=2,A(1,2)=0,A(2,2)=1,A(3,2)=0,A(4,2)=0,A(5,2)=1
tras ejecutar el programa, K=-3 (determinante de la matriz)
y los resultados:
A(3,0)=-0.333...,A(4,0)= 0       ,A(5,0)= 0.666...
A(3,1)= 0       ,A(4,1)= 1       ,A(5,1)= 0
A(3,2)= 0.666...,A(4,20= 0       ,A(5,2)=-0.333...
de lo que deducimos que
                 
        -1  0  2 
1/A=1/3* 0  3  0 
         2  0 -1 
                 

Variables que utiliza el programa:
E,F,I,J,K,Z