Altillo.com > Exámenes > UBA - Ingeniería > Análisis Numérico

Análisis Numérico I

Examen Final

Cat: Menéndez - Cavaliere - Perez Berro - Tarela

22 / 07 / 2002

Altillo.com

Problema 1)
Dado el siguiente problema de valores de contorno (PVC) definido en el intervalo [0 , 2]

con  y(0) = 1  e  y'(2) = -y(2)

a) Aplicar el método de las diferencias y plantear la ecuación en diferencias correspondiente a un error de discretización de orden 2.
b) Escribir un programa (en pseudo-código bien estructurado o en el lenguaje de programación empleado para realizar los TP's) que calcule e imprima la solución aproximada del PVC utilizando el planteo descripto en el punto anterior. Suponer que el único dato que se ingresa al programa es la cantidad N de intervalos en los que se divide el dominio del problema la cual se pide que no sea inferior a 2. El programa deberá llamar al procedimiento TRIDIAG (M , T) que resuelve sistemas de ecuaciones lineales (SEL) tridiagonales. La variable M, que es de tipo entero, almacena a la dimensión del SEL mientras que T es un arreglo de M filas y 5 columnas que almacena en variables de punto flotante los siguientes datos:

Columna 1: (i , 1) = ai , i-1 con i = 2, 3, ..., M donde a es la matriz de coeficientes del SEL.
Columna 2: (i , 2) = ai , i con i = 1, 2, ..., M
Columna 3: (i , 3) = ai , i+1 con i = 1, 2, ..., M-1
Columna 4: (i , 4) = bi con i = 1, 2, ..., M donde b es el vector de términos independientes del SEL.
Columna 5: (i , 5) = xi con i = 1, 2, ..., M donde x es la solución del SEL.

c) Calcular la solución aproximada para el caso de N = 2.

Problema 2)
a) Analizar si el siguiente método numérico para integrar problemas de valores iniciales (PVI) de primer orden es consistente e indicar cual es el orden de su error de discretización.

un+1 = un + k/3 [f(tn , un) + 2 f(tn + 3/4 k , un + 3/4 k f(tn, un))] donde y' = f(t,y)

b) Aplicar el método anterior para calcular la solución del siguiente PVI en t = 3 usando un paso k = 0.5

y' = 1 + (t - y)2 con  y(2) = 1 cuya solución es  y = t+1/(1 - t)

c) Obtener un resultado mas preciso sin reducir h usando extrapolación de Richardson.

Pregunta 1)
Indicar cual es el polinomio de mayor grado que puede ser integrado por los métodos de Rectángulo, Trapecio, Simpson y Gauss.

Pregunta 2)
¿Cuál método de extrapolación conviene ser usado si luego se planea aumentar el grado del polinomio interpolante a traves de la incorporación de nuevos puntos de interpolación. Justificar la respuesta.