Trabajo 11, abril 17

La semana que viene estaré fuera de la ciudad.
A mi regreso deberán entregar (no aceptaré su entrega un día después) el siguiente trabajo:

Implemente en Fortran uno de los tres siguientes algoritmos iterativos para resolver sistemas lineaes:

0. Gauss-Seidel
1. Jacobi
2. SOR

La numeración la debe seguir módulo 3 con respecto a la tabla de calificaciones última que aparece en este blog.

Pueden consultar el libro:
Principles of Petroleum Reservoir Engineering
Gian Luigi Chierici

O cualquier otro, incluso en google obtendrán información.

Debe presentar un bosquejo del algoritmo, en que se basa, dónde se usa y cómo detener el algoritmo. Deben presentar una rutina cuyas entradas, al menos para Jacobi son:
subroutine jacobi(lda,A,n,b,x,tol,maxiter,iband)

La rutina debe estar suficientemente documentada siguiendo las ideas de la documentación de cualquier rutina de LINPACK (seré estricto en esto).

Deben acompañar a la rutina un programa principal de prueba, donde la matriz A es de Hilbert, de orden 5 o 10 (le piden esa info al usuario), el lado derecho es el vector b=(1,2,3,...n) y pruebe el punto inicial x0=(1,1,...1).

Este trabajo es muy simpático y requiere mucha paciencia, creatividad y esmero de su parte.

NOTA: El trabajo lo pueden hacer entre 2 personas con la misma numeración modulo 3, solamente 2 personas.

Cualquier cosa... en el blog.

Les dejo el archivo donde vienen los métodos señalados de manera escueta: iterative_methods.pdf

Suerte.