El siguiente trabajo lo pueden leer de la página 60 del libro de Kahanner.
La inversa de una matriz cuadrada no singular, se puede calcular resolviendo una coleccion de sistemas lineales de la forma:
A x_i = e_i
donde e_i es el vector canónico, de entradas cero excepto en la i-ésima posición. x_i representa la i-ésima columna de la matriz inversa. Esto es, estamos calculando cada una de las columnas resolviendo un sistema de ecuaciones.
Debe hacer una rutina invert que calcule la matriz inversa de A, y de igual manera, calcule el número de condición de la matriz.
Debe usar las rutinas dgefa.for y degesl.for de la biblioteca linpack, así como las rutinas de la biblioteca blas se requieran por esas dos primeras. Estas biliotecas las pueden adquirir en Netlib.
Su ejercicio se puede titular: use bibliotecas finas escritas en fortran.
Suscribirse a:
Enviar comentarios (Atom)
2 comentarios:
que tal prof.
Pudiera subir la bibliografia completa del libro para el procediiento que usted suguiere
Es un libro clásico, basta escribir kahanner en google para que aparezca la referencia.
Computer Methods for Mathematical Computations.
Geroge E. Forsythe, Michael A. Malcolm and Cleve B. Moler.
Prentice Hall, 1977.
QA297.F568
ISBN 0-13-165332-6
Saludos.
Publicar un comentario