martes, 29 de abril de 2008

The Coffee cooling problem

El trabajo para el 6 de mayo, es investigar sobre el problema del enfriamiento de la taza de café.

Hay varias formulacioes del problema y algunas muy ingeniosas en Internet.

Debe indicar en qué consiste, y una posible solución al problema, ya sea en su forma continua o discreta. De ser posible trace algunas gráficas del comportamiento de la temperatura del café a lo largo del tiempo.

Revise la ley de enfriamiento de Newton... ahí está la clave.

domingo, 27 de abril de 2008

jueves, 17 de abril de 2008

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.

martes, 15 de abril de 2008

Trabajo 10, abril 15

En la página 100 del libro Scientific Computation, viene un ejercicio, el 2.3 sobre una estructura con fuerzas en equilibrio en los nodos.

En clase lo resolvimos e interpretamos su solución. Ahora bajo 3 diferentes numeraciones de las fuerzas que se ejercen sobre los nodos, identifique el sistema Ax=b asociado a casa uno de ellos, resuelva el sistema y comente sobre la sparcidad de la matriz y su número de condición cond(A).

Se entrega el jueves 17 de abril a primera hora.

miércoles, 9 de abril de 2008

Trabajo 9, abril 8

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.

jueves, 3 de abril de 2008

True or False

El trabajo para el próximo martes es por grupos. Los grupos con numeración módulo 6 de acuerdo a la lista de calificaciones, deben resolver los problemas con numeración módulo 6 de la lista de problemas:treuorfalse.pdf

No olviden que el martes pasaré a un representante a la pizarra a resolver uno de los problemas que en grupo presentan.

Les dejo las notas del capítulo2, conserven el pdf que seguiremos resolviendo ejercicios de ahí.
Chau.

miércoles, 2 de abril de 2008

Calificaciones de Trabajos, Abril 3

A continuación listo sus calificaciones, así como una nueva numeración para sus hacer trabajos módulo el número que se les indique.

Con el mouse pulse sobre la palabra: calificaciones, para acceder al archivo.