Analisis Numerico I (75.12-95.04-95.

13)
Profesora responsable: Miryam Sassano

Bibliografa

Burden R.L., Faires J.D. Analisis Numerico, Grupo Editorial Iberoamericano 1985.

Chapra S., Canale R.Metodos Numericos para Ingenieros, Mac Graw Hill 1985

Kincaid D., Cheney W. Analisis Numerico: las matematicas del calculo cientfico, Addison Wesley, 1994.

Zill, D. G. (2007). Ecuaciones Diferenciales con Aplicaciones, Grupo Editorial Iberoamericana

Mathews, J. H., Fink, D. K., Metodos Numericos con Matlab, Tercera Edicion, Editorial Prentice Hall, 2000.

Nakamura, S., Analisis Numerico y Visualizacion Grafica con Matlab, Prentice Hall, 1997

i
Indice general

Gua I: Teora de Errores 1

Gua II: Ecuaciones no lineales 4

Gua III: Interpolacion 7

Gua IV: Ajuste-Cuadrados Mnimos 10

Gua V: Sistemas de ecuaciones lineales 12

Gua VI: Sistemas de ecuaciones no lineales 16

Gua VII: Diferenciacion e Integracion numerica 17

Gua VIII: Ecuaciones Diferenciales 20

Gua IX: Ecuaciones diferenciales de segundo orden y problemas de valores de contorno 21

ii
 Gua I: Teora de Errores

Teora de Errores

1. Si a y b son los valores que tienen dos magnitudes y a y b son los valores medidos de dichas magnitudes, obtener
para cada una de ellas:

a) El error absoluto, el error relativo y una cota para el error relativo.

b) Obtener el er(a+b) con estos datos:a = 2, a = 2.01, b = 3 y b = 3.02.
c) Calcular el er(ab) con los datos de b.
d ) Obtener el er(a.b) cuando a = 2, a = 2.02, b = 1 y b = 1.6 . Redondear los resultados con 2 decimales.

2. Si el resultado de medir cierta longitud con una regla graduada en milmetros es x = 3cm. Dar una cota para el
error absoluto de y dar una cota para el error relativo de x.

3. Suponga que el resultado de una operacion dio: x = 0.2489, |ex | = 0.6104

a) Es posible afirmar que que tiene t = 4 decimales significativos o correctos? Justifique la respuesta.
b) Es posible afirmar que tiene t = 3 decimales significativos o correctos? Justifique la respuesta. Si es as,
diga como se debe expresar.

4. Calcular con tres decimales correctos o significativos las siguientes expresiones:

a) 1.3134
b) 0.3751 e
c) e

5. Hallar cotas para los errores inherentes propagados ( absolutos y relativos ) en los siguientes calculos, donde:
x = 2.00 , y = 3.00 y z = 4.00.

a) f = 3x + y z
xy
b) z
y
c) f = xsen( 40 )

6. Siendo x = 2.0 0.1 , y = 3.0 0.2 y z = 1.0 0.1 hallar una cota para el error absoluto inherente propagado
de la siguiente expresion:
2
xy
w= .
z

1
2

7. Calcular la expresion:

a) f = ( 2 1)6 utilizando para 2 el valor aproximado 1.4. Sustituir el valor aproximado en a) y en cada
una de las expresiones siguientes:
b) f = 1
( 21)6

c) f = (3 2 2)3
1
d) f = (32 2)3

e) f = 99 70 2
1
f) f = 9970 2

8. Se estan realizando observaciones de un satelite para determinar su velocidad. En la primera observacion la

distancia medida al satelite fue r = 30000 10km. Cinco segundos mas tarde se determina un aumento en la
distancia r = 125.0 0.5km y el cambio en la orientacion resulto de = 0.00750 0.00002radianes. Hallar
la velocidad del satelite suponiendo que el mismo se mueve en lnea recta y a velocidad constante durante ese
intervalo de tiempo; indicar la precision del resultado. Considerar exacto el lapso de 5 segundos.

9. Se dispone de un algoritmo para computar la siguiente integral:

1
ebx
Z
I(a, b) = dx
0 a + x2

Utilizando dicho algoritmo se obtuvo la siguiente tabla:

a b I
0.39 0.34 1.425032
0.40 0.32 1.408845
0.40 0.34 1.398464
0.40 0.36 1.388198
0.41 0.34 1.372950
Se pretende estimar I(z, y), siendo y, z cantidades fsicas obtenidas de un proceso de medicion: z = 0.400 0.003
e y = 0.340 0.005. Aproximar el numero buscado y dar una cota para el error inherente propagado.

10. Con cuantas cifras decimales se debe usar a 2 en la expresion z = (3 2 2)3 para que resulte |ez | 0.5103 ?

11. Con cuantas cifras decimales se debe usar a 2 en la expresion z = ( 2 1)6 para que resulte |ez | 0.5103 ?

12. Calcular, si es posible, con tres decimales correctos el resultado de la siguiente expresion: x = a con a = 1.3134
correctamente redondeado y = 3.14 redondeado a dos decimales. Cual debe ser el error de para que se
cumpla la condicion que |ex | 0.5103 ?

13. Calcular, si es posible, con 4 decimales correctos el resultado de x = a.e con a = 1.37514, correctamente
redondeado y e = 2.718. Cual debe ser el error de e para que se cumpla la condicion que la expresion del |ex |
sea menor o igual a la cota del error por redondeo del resultado en el ejercicio anterior?
 3

14. Dibujar las graficas de proceso, determinar la condicion del problema y el termino de estabilidad para los
siguientes casos:

a) u = a + a
b) v = 2a
c) u = a + a + a
d ) v = 3a

Suponer que a posee cierto error inherente y que los errores inherentes de 2 y 3 son nulos y pueden ser repre-
sentados exactamente. Comparar los resultados y extraer conclusiones. Efectuar los calculos con a = 0.6992 en
el sistema [G(10, 4, 1), red].

15. Demostrar con los siguientes ejemplos que en [G(10, 5, 1), red] no valen la ley asociativa ni la distributiva:

a) 0.98765 + 0.012424 0.0065432

b) (4.2832 4.2821)5.7632

16. Calcular u2 v 2 en el sistema de numeracion [G(10, 4, 1), red] con u = 43.221 y v = 43.11, utilizando los siguientes
algoritmos:

a) uu vv
b) (u + v)(u v)

Indicar cual algoritmo es mas conveniente y justificar.

17. Indicar cual de los siguientes algoritmos es mas estable para calcular la menor raz de la ecuacion:

x2 2x + a = 0 0<a1

a) 1 = 1 a; 2 = 1 ; x = 3 = 1 2 ;
a
b) 1 = 1 a; 2 = 1 ; 3 = 1 2 ; x = 4 = 3

18. Preguntas Teoricas.

a) Identificar y describir las principales fuentes de errores.

b) Indicar que significa que un algoritmo es menos estable que otro. Ejemplificar.
c) Indicar como hara para determinar la influencia de los errores de redondeo sobre los resultados de un
algoritmo utilizando la computadora.
d ) Pretendemos evaluar:

X xn
S(x) =
n=0
n!
Hallar estimaciones de los errores de truncamiento, de los inherentes propagados y de los de redondeo
propagados al efectuar los calculos.
 Gua II: Ecuaciones no lineales

Ecuaciones no lineales

1. Dada la funcion f (x) = ex (sen(x) + cos(x) 2x 2), hallar una de las races en el intervalo [2.5; 0.5] por el
metodo de la biseccion, trabajando con 5 (cinco) decimales, con una tolerancia de 105 . Informar la cantidad de
iteraciones que son necesarias para obtener la tolerancia indicada.

2. Las siguientes ecuaciones tienen una raz en el intervalo (0; 1, 6). Determinarlas con un error menor que 0.02 por
el metodo de la biseccion.

a) x.cos(x) = ln(x).
b) 2.x ex = 0.
c) e2x = 1 x.
x2
3. Sea f (x) = 4 sen(x), se desea encontrar la primer raz positiva de f (x).

a) Hallar un intervalo de partida para utilizar el metodo de la biseccion.

b) Estimar el numero de aproximaciones necesarias para hallar la raz con una tolerancia para el error absoluto
de 0.02. Calcular la raz.
c) Si la tolerancia de 0.02 es sobre el error relativo, cuantas aproximaciones se requieren?
d ) Sabiendo que la raz buscada a 5 decimales es = 1.93375 obtener conclusiones sobre la performance del
metodo.
x2
4. Aplicar el metodo del punto fijo para encontrar la mejor aproximacion de la raz de la funcion 4 sen(x) en el
6
intervalo I = [1.5; 2] con 6 decimales y tolerancia de 10 . Cuantas iteraciones se necesitan? Hacer graficos de
g(x), la funcion identidad y los sucesivos terminos de la sucesion generada.
x
5. La ecuacion e 4 = x tiene dos races reales.
xk
a) Verificar que solo una de ellas puede obtenerse con el metodo de punto fijo en la forma xk+1 = e 4 . Explicar
porque no puede obtenerse la otra raz de esta manera.
b) Obtener una aproximacion a la raz que s puede obtenerse con ese esquema iterativo con un error menor
que 103 . Representar graficamente la marcha hacia el punto atractor verificando que es en escalera.

x41
6. Tiene races la funcion f (x) = x2 25 ? Fundamentar.

4
 5

7. Se desea hallar la primera raz positiva de la ecuacion x = cos(x) con el metodo de Newton Raphson.

a) Plantee el metodo iterativo correspondiente para el problema, usando el metodo de punto fijo.
b) Estudie las propiedades de convergencia del metodo propuesto. Encuentre explcitamente un intervalo de
convergencia.
c) Encuentre el cero buscado con una tolerancia para el error relativo de 1010 .

8. Estudiar la convergencia del metodo de Newton Raphson aplicado a la ecuacion x2 1 = 0. Elegir como valor
inicial x0 = 2 y calcular aproximaciones de la raz con precision sucesivamente creciente.

9. Determinar la raz no nula de la ecuacion x = 1 e2x , usando el metodo de Newton Raphson con 4 decimales
significativos. Verificar las condiciones de convergencia del metodo en el intervalo elegido.

10. Determinar la raz de la ecuacion xln(x) 1 = 0, usando el metodo de Newton Raphson con 5 decimales
significativos.

11. Aplicar el metodo de Newton Raphson para determinar una raz compleja de la ecuacion x2 + 1 = 0; comenzar
las iteraciones con x0 = 1 + i.

12. Sea la ecuacion f (x) = x3 9x2 + 24x 20 = 0.

a) Aplicar el metodo de Newton Raphson con x0 = 1, 5. Detener el proceso cuando se obtengan 2 decimales
significativos correctos.
b) Aplicar el metodo de Newton Raphson para el caso de races multiples y las mismas condiciones del punto
a.

13. Hallar la raz negativa de la funcion f (x) = x2 x 2, utilizando el metodo de Newton Raphson y aritmetica
de 4 digitos. Estimar el error que se comete entre dos iteraciones consecutivas.

14. Considerar la ecuacion ex = sen(x).

a) Verificar que esta ecuacion tiene infinitas soluciones reales negativas. Que puede concluirse de la distancia
entre soluciones consecutivas muy alejadas de x = 0?
b) Usando el metodo de Newton obtenga aproximaciones con error menor que 106 para las tres races mas
cercanas a x = 0.

15. Se desea hallar la raz negativa de la funcion f (x) = x3 1, 9x2 1, 05x+2, 745 con 6 dgitos de precision. Utilizar
el metodo de Newton Raphson, partir de x0 = 1 y no superar las 10 iteraciones, verificar que la convergencia
es lineal.
x2
16. Aplicar el metodo de la secante para hallar la raz no nula de f (x) = 4 sen(x) con una tolerancia del 0.1 %.

17. Usando el metodo de Newton Raphson, hallar la raz negativa de f (x) = 3x2 3. Luego comparar con los
otros metodos: biseccion, regula falsi , secante y punto fijo. Cual metodo es mejor? Trabajar con 4 decimales y
redondeo, con tolerancia de 104 .
1
18. La funcion f (x) = sen(x) 2 x tiene 2 ceros en el intervalo [0; 2]. Uno es x = 0, se desea hallar el otro. Para
ello se utiliza un metodo de punto fijo basado en la funcion de iteracion g(x) = x f (x). Las figuras muestran
g y g 0 en el intervalo [0; 2].

a) Hallar, un intervalo que contenga al cero buscado como unico cero de f (x).
6

b) Mostrar que en dicho intervalo el metodo propuesto converge.

c) Hallar el cero con una tolerancia del 1 % para el error relativo entre dos pasos consecutivos.
d ) Hallar el orden de convergencia del metodo y la constante asintotica del error.

19. Dada la profundidad h y el perodo T de una ola, su longitud de onda l surge de la relacion de dispersion
2. 2.
2 = g.k.tangh(k.h), donde = T es la pulsacion, g es la aceleracion de la gravedad y k = l es el numero
de onda. Conociendo g = 9.8 sm2 y h = 4m, se desea calcular cual es la longitud de onda correspondiente a un ola
con T = 5seg.

a) Utilizar el metodo de punto fijo para calcular la solucion.

b) Utilizar el metodo de Newton Raphson para calcular la solucion con 4 dgitos de precision.
c) Partir del resultado obtenido en a).

20. Para un tiro oblicuo considerando el amortiguamiento viscoso del aire (fuerza resistente de modulo proporcional
a la primera potencia de la velocidad) el alcance L satisface:

AL + Bln(1 CL) = 0

Determinar el alcance L (con error relativo menor al 1 %) para los siguientes casos:

a) A = 2; B = 10; C = 0.1.
b) A = 4; B = 10; C = 0.1.
 Gua III: Interpolacion

Interpolacion

1. Encontrar el polinomio de grado 3 que pasa por los siguientes puntos utilizando la formula de Lagrange:
x 0 1 2 4
y 1 1 2 5

2. Hallar los valores de 1, 01 y 1, 28 a partir de la siguiente tabla, por interpolacion de Lagrange y de Newton
con tres dgitos significativos.
x 1, 00 1, 05 1, 10 1, 15 1, 20 1, 25 1, 30

x 1, 00000 1, 02470 1, 04881 1, 07238 1, 09544 1, 11803 1, 14017

3. Calcular f (3) utilizando la formula de Newton, dada la siguiente tabla:

x 1 2 4 5
f (x) 0 2 12 21

a) Tomar los puntos 1,2 y 4 luego los puntos 2,4 y 5.

b) Repetir a pero usando el polinomio de Lagrange.
c) Aproximar por un polinomio de grado 3.

4. Calcular f (0) utilizando la formula de Newton, dada la siguiente tabla:

x 0, 1 0, 2 0, 4 0, 8
f (x) 64987 62055 56074 43609
Notar que la formula de interpolacion se utiliza para extrapolar.

5. Encontrar el polinomio de grado 3 que pasa por los siguientes puntos utilizando la formula de Newton:
x 1 1, 25 1, 50 1, 75 2, 00
y 5, 10 5.79 6.53 7.45 8.46

6. Hallar un polinomio Q de grado 3 tal que: Q(0) = 0, Q0 (0) = 1, Q(1) = 3, Q0 (1) = 6.

7. Se conocen los siguientes datos acerca de la funcion f (x): f (0) = 1, f (1) = 2, f (2) = 5, f 0 (0) = 0 y f 0 (2) = 4.

a) Hallar el polinomio interpolante que verifica esos datos mediante el metodo de Hermite.
b) Hallar la funcion Spline de orden 2 que verifica esas condiciones.

7
8

8. Aproximar f (x) = sen(ex 2) en x = 0.9 con los siguientes datos:

x f (x) f 0 (x)
0.8 0.22363362 2.1691753
1.0 0.65809197 2.0466965

1
9. Hallar el polinomio interpolante de grado 2 para f (x) = x por medio de la formula de Lagrange, utilizando los
nodos: x0 = 2, x1 = 2.5 y x2 = 4. Graficar la curva y su aproximacion. Analizar los errores para x = 0.5 y x = 31 .

10. Un coche que viaja en una carretera recta es cronometrado en algunos puntos. Los datos obtenidos se dan en
la siguiente tabla. Utilice un polinomio de Hermite para predecir la poscion del coche y su velocidad cuando
t = 10 segundos.
T iempo (seg) 0 3 5 8 13
Distancia (m) 0 67, 5 114, 9 186, 9 297, 9
V elocidad (m/s) 22, 5 23, 1 24 22, 2 21, 6

11. Se desea hallar una funcion polinomica para aproximar a la funcion f (x) = ex cos(x), en el intervalo [0; 2]

a) Tabular f (x) en los nodos x = 0, 0.5, 1 y 2 y hallar el polinomio interpolante por el metodo de Newton,
trabajar con 5 dgitos.
b) Agregar el nodo x = 1.5 para hallar una expresion aproximada para el error. Utilizarla para estimar el error
en x = 0.1, 0.3, 1.2, y 1.7.
c) Comparar los errores estimados en el punto anterior con los valores correctos, calculados como diferencia
entre el valor correcto de f (x) y el obtenido por medio del polinomio interpolante.

12. Se tiene la funcion f (x) = ex , de la cual se proveen los siguientes valores:

x 0 0.5 1 2
f (x) 1 1, 64872 2, 71828 7, 38906

a) Estimar f (0.25) utilizando interpolacion de Lagrange con los nodos x0 = 0.0 y x1 = 0.5.
b) Estimar f (0.75) utilizando interpolacion de Lagrange con los nodos x0 = 0.5 y x1 = 1.0.
c) Estimar f (0.25) y f (0.75) utilizando interpolacion de Lagrange con los nodos x0 = 0.0, x1 = 0.5 y x2 = 2.0.
d ) Estimar los errores de truncamiento de los calculos realizados en los punto a), b) y c) en base a la formula:

f (n+1) ()
f (x) = f (x) + (x x0 ) . . . (x xn )
(n + 1)!
e) Compararlos con los valores exactos calculados a partir de los valores reales de la funcion f (0.25) = 1.28403
y f (0.75) = 2.11700.
f ) Indicar que aproximaciones resultaron mas precisas y por que.

13. Se desea hallar una funcion de interpolacion polinomica para aproximar la funcion sen2 (x) en el intervalo [0; ]

a) Construir un polinomio de Hermite en los nodos: 0, 2 y trabajar con cuatro decimales. Estimar el error
cometido en la construccion del polinomio.
b) Estimar el error cometido en 0.2, 0.5 y 1 utilizando el punto extra: 5 4 en la tabla de interpolacion de
Hermite.
 9

14. Se desea interpolar una Spline cubica para una funcion tabulada en 4 nodos. Explicar cuantas son las incognitas
y cuales las ecuaciones que completan el planteo del problema.

15. Construya el trazador cubico libre con los siguientes datos:

x f (x)
a) 8.3 17.56492
8.6 18.50515

x f (x)
0.1 0.62049958
b) 0.2 0.28398668
0.3 0.00660095
0.4 0.24842440

16. Los datos del ejercicio anterior se generaron usando las siguientes funciones. Utilice los trazadores cubicos con-
0
struidos en el ejercicio anterior a fin de aproximar f (x) y f (x). Calcule el error.
0
a) f (x) = x.ln(x), aproxime f (8.4) y f (8.4).
0
b) f (x) = x.cos(x) 2x2 + 3x 1, aproxime f (0.25) y f (0.25).

17. Un trazador cubico sujeto S de la funcion f esta definido en el intervalo [1; 3] por:
(
S0 (x) = 3(x 1) + 2(x 1)2 (x 1)3 si 1 x 2
S= 2 3
S1 (x) = a + b(x 2) + c(x 2) + d(x 2) si 2 x 3
0 0
Dadas f (1) = f (3), encontrar a, b, c y d.

18. Un trazador cubico natural S esta definido por:

(
S0 (x) = 1 + B(x 1) D(x 1)3 si 1 x 2
S= 3 2 3
S1 (x) = 1 + b(x 2) 4 (x 2) + d(x 2) si 2 x 3

Si S interpola los datos (1; 1), (2; 1) y (3; 0) obtener: B, D, b y d.

 Gua IV: Ajuste-Cuadrados Mnimos

Ajuste-Cuadrados Mnimos

1. Aproximar los datos con un polinomio de grado 2, por cuadrados mnimos y graficar la solucion:
x 0.00 0.25 0.50 0.75 1.00
y 1.0000 1.2840 1.6487 2.1170 2, 7183

a) Calcular los errores para cada dato de la tabla.

b) Calcular el polinomio interpolante de Lagrange de grado 2, en los nodos 0, 0.5 y 1.
c) Graficar la solucion y calcular los errores para cada dato de la tabla.
d ) Comparar los resultados obtenidos y determinar cual solucion aproxima mejor a la curva en el intervalo
[0; 1].
e) Comparar los resultados obtenidos con la funcion f (x) = ex en los puntos: 0.1, 0.2, 0.3, 0.6, 0.7, 0.8 y 0.9,
calcular los errores y obtener conclusiones.

2. Determinar las lneas rectas que aproximan la curva y = ex , segun los siguientes metodos, comparar los resultados
y calcular los errores en x = 1. Utilizar 3 decimales y redondeo.

xi 1 0.5 0 0.5 1
a) Cuadrados mnimos sobre la siguiente tabla:
yi 0, 368 0, 607 1 1, 649 2, 718

b) Tomando la lnea tangente a y = ex en el punto medio del intervalo [0; 1], es decir, la aproximacion de
Taylor de primer orden.
c) Tomando la lnea tangente a y = ex en el punto medio del intervalo [1; 1], es decir, la aproximacion de
Taylor de primer orden.

3. Se tiene la siguiente tabla de datos:

x 6 8 10 12 14 16 18 20 22 24
y 3.8 3.7 4.0 3.9 4.3 4.2 4.2 4.4 4.5 4.5

a) Encontrar una funcion lineal que aproxime estos datos por cuadrados mnimos. Utilizar esta curva para
suavizar los datos.
b) Repetir el punto anterior con una funcion cuadratica.
c) Comparar los resultados.

10
 11

x 1 2 3 4
4. Obtener una formula del tipo: f (x) = a.emx a partir de los siguientes datos:
y 7 11 17 27

5. Dada la siguiente coleccion de datos, elegir una curva de aproximacion y analizar los errores respecto de los
valores dados:
x 1.00 1.25 1.50 1.75 2.0
y 5.10 5.79 6.53 7.45 8.46

6. Construir las aproximaciones indicadas, calcular los errores y obtener conclusiones:

a) Aproximacion polinomica de grado 1.

b) Aproximacion polinomica de grado 2.
c) Aproximacion polinomica de grado 3.
d ) Aproximacion de la forma: b.eax .
e) Aproximacion de la forma: b.xa .
Para las siguientes tablas:

x 4.0 4.2 4.5 4.7 5.1 5.5 5.9 6.3 6.8 7.1
i)
y 102.56 113.18 130.11 142.05 167.53 195.14 224.87 256.73 299.50 326.72

x 0.2 0.3 0.6 0.9 1.1 1.3 1.4 1.6

ii)
y 0.050446 0.098426 0.332770 0.726600 1.097200 1.569700 1.848700 2.501500

7. Para 5 instantes de tiempo se observaron los siguientes valores de un parametro fsico:

t 2 1 0 1 2
u u2 u1 u0 u1 u2
Mostrar que, si los datos se ajustan por una parabola (t), la aproximacion en t=0 es:
1
(0) = (3u2 + 12u1 + 17u0 + 12u1 3u2 )
35

8. Hallar el polinomio aproximante de segundo grado para la funcion: f (x) = sen(.x) en el intervalo [0; 1]. Graficar
la funcion y su aproximacion.

9. El nivel de agua en el Mar del Norte esta determinado principalmente por la marea llamada
M 2, cuyo perodo es de aproximadamente 12 horas. Se han realizado las siguientes mediciones:
t (horas) 0 2 4 6 8 10
H(t) (m) 1, 0 1, 6 1, 4 0, 6 0, 2 0, 8

a) Ajustar la serie de mediciones usando el metodo de los cuadrados mnimos y la funcion:

2..t
H (t) = h0 + a1 sen( )
12
b) Calcular errores que permitan estimar la precision de la aproximacion realizada en a).
c) Utilizar ahora la funcion:
2..t 2..t
H (t) = h0 + a1 sen(
) + a2 cos( )
12 12
Repetir b) para la nueva funcion aproximante. Comparar y obtener conclusiones.
 Gua V: Sistemas de ecuaciones lineales

Sistemas de ecuaciones lineales

1. Resolver el sistema A.x = b utilizando eliminacion de Gauss sin pivoteo, donde:

1 2 3 4 2

1 4 9 16 , b = 10

A= 1 8
27 64 44

1 16 81 256 190

2. Calcular la inversa de la matriz A resolviendo el sistema A.X = I, utilizando eliminacion de Gauss, siendo I la
matriz identidad y X la matriz inversa de A.

2 1 2

A=
1 2 3

4 1 2

3. Dada la siguiente descomposicion LU de la matriz A efectuada utilizando pivoteo parcial:

1 0 0 4 1 0 2

L = 0.5 1 0 , U = 0 2.5 1 p = 3


0.5 0.2 1 0 0 8 1

a) Resolver el sistema de ecuaciones A.x = b, siendo:

b = (1 2 7)t

b) Obtener la matriz A y verificar la solucion obtenida en a).

4. Considerar la matriz A definida segun:

1
ai,j = 1 i; j 4
i+j1
Considerar el sistema Ax = b donde:

b = (0.58333 0.21667 0.11666 0.07381)t

Resolver el sistema utilizando eliminacion de Gauss con pivoteo parcial operando con 5 decimales. Investigar las
caractersticas de la matriz y obtener conclusiones.

12
 13

5. Dada la matriz A del problema anterior y

b = (2.66666 1.50000 1.06666 0.83334)t

Resolver Ax = b aplicando la descomposicion LU de A, trabajando con 5 decimales y redondeo. Obtener

conclusiones.

6. Sea el sistema de ecuaciones lineales:

" # " # " #
0.003153 15.28 x 14.98
. =
0.009413 45.6 y 44.75

a) Obtener la solucion numerica utilizando dos algoritmos: eliminacion de Gauss con pivoteo parcial y elimi-
nacion de Gauss con pivoteo total.
b) Estimar el numero de condicion de la matriz de los coeficientes.
c) En base a los resultados obtenidos en los puntos a) y b), indicar cuales de las siguientes afirmaciones son
correctas y por que:
1) El primer algoritmo esta mal condicionado.
2) El segundo algoritmo esta mal condicionado.
3) El problema esta mal condicionado.

7. Describir como se simplifica el algoritmo del metodo de eliminacion de Gauss para el caso particular en que la
matriz de coeficientes es simetrica definida positiva.

8. Dado el siguientes sistema de ecuaciones lineales donde la matriz A es no singular:

" # " # " #
a11 a12 x b1
. =
a21 a22 y b2

a) Establecer cuando el metodo de Jacobi diverge.

b) Demostrar que si el metodo de Jacobi converge, entonces el metodo de Gauss Seidel lo hace mas rapido.

9. Sea el sistema de ecuaciones lineales:

" # " # " #
0.01235 2.387 x 1.37
. =
5.462 0.008406 y 10.85

a) Resolverlo por el metodo de Jacobi. Efectuar las modificaciones necesarias para garantizar la convergencia.
Trabajar con 5 decimales y redondeo.
b) Explicar la convergencia o no de los algoritmos del punto a) en terminos de la norma de la matriz de
iteracion.

10. Resolver el siguiente sistema utilizando el metodo de Gauss-Seidel, iterando hasta que la maxima diferencia entre
dos valores sucesivos de x, y o z sea menor que 0.02. Indicar si esto ultimo significa que la solucion obtenida
esta en un intervalo de radio 0.02 alrededor de la solucion exacta.

10x + 2y + 6z = 28


x + 10y + 4z = 7

2x 7y 10z = 17

14

11. Resolver el siguiente sistema utilizando el metodo de Gauss-Seidel:

a+d=2


a + 4b d = 4


a+c=2

c+d=2

12. Considerar el sistema poco denso de ecuaciones:

2a b = 1


a + 2b c = 1


b + 2c d = 1

c + 2d = 1

Mostrar que el sistema permanece poco denso cuando se lleva a la forma triangular utilizando el metodo de
eliminacion de Gauss. Hallar la solucion por Gauss y luego por Gauss-Seidel.

13. Dado el siguiente sistema de ecuaciones:

3.21x + 0.943y + 1.02z = 2.3
0.745x 1.29z = 0.74

0.875x 2.54y + 0.247z = 3.39

a) Efectuar las modificaciones necesarias para poder garantizar la convergencia utilizando el metodo de Gauss-
Seidel.
b) Resolverlo iterando hasta obtener una precision de 3 dgitos significativos, sin exceder un maximo de 5
iteraciones.
c) Determinar como influye un error absoluto de 0.01 en el primer coeficiente de la primera ecuacion sobre los
valores calculados de x, y y z.

14. Dado el sistema A.x = b, construir un algoritmo que halle el vector solucion mediante el metodo de Gauss-Seidel.

15. Sea el siguiente sistema:

5x + 4z = 1
3x + 6y = 1.20


9y + 10z = 0.8

Que metodo usara para resolverlo y por que?

16. El siguiente codigo en Octave implementa la aproximacion a la solucion de un sistema de 5 ecuaciones lineales
con 5 incognitas con un metodo iterativo:
e = 1; x = zeros(1,5); y = zeros(1,5)
k = 0; s =103
while e > s
y(1) = (5 + x(2))/10
y(2) = (4 + x(1) + x(3))/10
y(3) = (4 + x(2) + x(4))/10
y(4) = (4 + x(3) + x(5))/10
y(5) = (5 + x(4))/10
 15

f=abs(y(1)-x(1))
for i =2:5
if abs(y(i)-x(i))>f
f=abs(y(i)-x(i))
end
end
e=f
x=y
k = k +1
end

a) Cual es el sistema de ecuaciones y cual el metodo iterativo? Justifique su respuesta.

b) Cual es el criterio de parada y como esta implementado? Cuales son los significados de las variables e y
k? Nota: zeros(n,m) genera una matriz nula de n filas por m columnas.

17. El siguiente codigo implementa una aproximacion a la solucion de un sistema de ecuaciones lineales:
xn1=1;xn2=1;xn3=0;xn4=0
i=0;e=1
tol=0.01
while (e > tol)
i=i+1;x1=xn1; x2=xn2; x3=xn3;x4=xn4
xn1=(13 - x2 + x3 - x4)/4
xn2=(-8 - x1 + x3 + 2 x4)/(-5)
xn3=(-2 - 2* x1 + x2 - 2 x4)/(-6)
xn4=(-5 - x1 + x2 - x3)/(8)
e=max([abs(xn1-x1) abs(xn2-x2) abs(xn3-x3) abs(xn4-x4)])
end

a) Escriba el sistema de ecuaciones lineales. Que metodo esta implementado? Sera convergente el pro-
cedimiento? Explique.
b) Modifique algunas lneas del codigo para que quede implementado el otro tipo de metodo iterativo presentado
en este curso. Haga dos pasos de este otro metodo con su calculadora y presente su respuesta en una tabla
redondeando a dos decimales pero operando con toda la precision de su calculadora.
 Gua VI: Sistemas de ecuaciones no lineales

Sistemas de ecuaciones no lineales

1. Sea el siguiente sistema de ecuaciones no lineal:

(
f (x, y) = x2 + y 2 4 =0
g(x, y) = x.y 1 =0

Resolverlo por el metodo de Newton con:(x0 y 0 )t = (2 0)t

2. Resolver el siguiente sistema usando el metodo de Newton, trabajar con 4 decimales y redondeo partiendo de:
(x0 y 0 )t = (0.3 0)t
2
(
1.021 xy = 4.953
5.040x.y = 0.05440

3. Resolver el siguiente sistema usando el metodo de Newton, trabajar con 4 decimales y redondeo partiendo de:
(x0 y 0 )t = (1 2)t (
3.11x.(y 1) = 8.73
0.749x + 121y = 2.08

4. Resolver el siguiente sistema de ecuaciones no lineales por el metodo de Newton trabajando con tres dgitos
significativos, partiendo de la aproximacion: (x0 y 0 )t = (1 3)t .
(
x2 y 2 = 11.20
x + y = 1.83

5. Repetir el ejercicio anterior usando el metodo del descenso mas rapido.

~ (0) = (1 1 3). Realizar

6. Resolver el siguiente sistema no lineal usando el metodo de Newton, partiendo de: X
tres iteraciones
x1 x2 x3 = 4.188


x1 + x2 + x3 = 3.677


x + 1.258x = 0
1 2

7. Resolver el siguiente sistema usando el metodo de Newton, hasta obtener una precision tal que ||ek+1 || < 106
(tomar la norma infinito), ek+1 es el error absoluto entre dos iteraciones consecutivas.

1
3x cos(y.z) 2 = 0


x2 81(y + 0.1)2 + sen(z) + 1.06 = 0

exy + 20z + 103 = 0

3

16
 Gua VII: Diferenciacion e Integracion numerica

Diferenciacion e Integracion numerica

1. Usar las formulas de diferencia progresiva y de diferencia regresiva para determinar las aproximaciones con que
se completan las siguientes tablas:
0
x f (x) f (x)
0.5 0.4794
a)
0.6 0.5646
0.7 0.6442
0
x f (x) f (x)
0.0 0.00000
b)
0.2 0.74140
0.4 1.3718

2. Los datos del ejercicio 1) se tomaron de las siguientes funciones. Calcular los errores reales y obtener una cota
para el error:

a) f (x) = sen(x)
b) f (x) = ex 2x2 + 3x 1

3. Se tiene la siguiente tabla de valores para la funcion seno:

x sen(x)
0.920 0.79560
0.950 0.81342
1.00 0.84147

a) Estimar el valor de la derivada de la funcion en x = 1 utilizando dos aproximaciones en diferencias en atraso

y luego obtener un valor mas preciso por extrapolacion de Richardson.
b) Construir una formula de aproximacion en diferencias de segundo orden para la derivada en x = 1, y
utilizarla para hallar un nuevo valor.
c) Discutir sobre la precision de los valores hallados en los puntos anteriores. Estimar el error de truncamiento
y compararlo con el valor exacto de ese error, obtenido por diferencia con el valor de cos(1).

17
18

4. Calcular la siguiente integral utilizando las formulas del trapecio y de Simpson con pasos h = 4 y h = 2,
respectivamente: Z 2
sen(x)
I= e(2 2 )
dx
2 0
Obtener conclusiones sobre la precision obtenida.

5. Integrar la funcion f (x) = x entre los argumentos 1.00 y 1.30, segun las formulas del trapecio y de Simpson.
Obtener conclusiones sobre la precision obtenida.

6. Integrar la funcion f (x) = sen(x) entre 0 y 2 a partir de la siguiente tabla:
5
x 0 12 6 4 3 2 2
sen(x) 0.0000 0.2587 0.5000 0.7071 0.8660 0.9659 1.0000

Usar el metodo del trapecio y el metodo de Simpson. Comparar los resultados con el valor exacto.

7. Integrar la funcion f (x) = sen(x)

x entre 0 y 0.8 por el metodo de Romberg con 5 dgitos significativos.
R2
8. Evaluar la integral I = 1 dx x mediante el metodo de Romberg hasta obtener 4 dgitos de precision.
R3
9. Hallar una formula de cuadratura para la integral I = 0 f (x)dx utilizando 4 nodos equiespaciados e interpolacion
polinomial sobre el intervalo. Utilizar el metodo de los coeficientes indeterminados. Si M es una cota superior
para |f (0)IV |, hallar una estimacion del error por truncamiento en terminos de M .
R5
10. Evaluar la integral: I = 1 ln(x)dx utilizando la formula de Simpson con un error no mayor a 0.01. Sabiendo
que la regla de Simpson aplicada a un intervalo generico (xi1 ; xi+1 ) es:
Z xi+1
h h5
f (x)dx = [f (xi1) + 4f (xi ) + f (xi+1 )] f IV ()
xi1 3 90

Donde xi1 <  < xi+1 se pide:

a) Obtener una expresion para acotar el error de truncamiento global sobre todo el intervalo de integracion.
b) Utilizando la expresion hallada en el punto anterior, determinar un paso h que garantice la cota de error
establecida y efectuar el calculo utilizando este valor de h.
c) Comparar el resultado de la integracion numerica con el valor exacto de la integral y verificar que la
diferencia esta acotada por el error estimado.
R
11. Evaluar la integral: I = 0 ln(2 cos(x))dx

a) Utilizar el metodo de Romberg comenzando con el paso h = 2. Trabajar con 5 decimales de precision.
Afinar el paso de calculo no mas de 2 veces.
b) Calcular la influencia de los errores en los valores del integrando sobre los valores obtenidos por la Regla
del Trapecio.
c) El valor exacto de la integral es:
3
ln(1 + ) = 1.959759164
2
Dar una explicacion de la alta precision obtenida con la Regla del Trapecio.
1
12. Integrar usando la regla del trapecio la funcion f (x) = x entre 2 y 5 con un incremento h de 0.2. Repetir usando
el metodo de Romberg para las filas que usted crea necesarias y que acerquen con un error relativo menor a 106
el resultado con respecto al valor analtico de la integral.
 19

13. Repetir el ejercicio anterior para la funcion f (x) = ex x2 entre 3 y 6.

14. Integrar por la regla de los trapecios compuesta la funcion f (x) = sen(x) + e2x con x en el intervalo [0.80, 2].
Hacerlo con n = 6. Repetir el ejercicio usando la regla de Simpson 1/3. Comparar errores.

15. Integrar mediante la regla de Simpson 1/3 compuesta la funcion del ejercicio 12, comparar con el resultado
analtico. Usar h = 0.5.

16. Integrar mediante la regla de Simpson 1/3 compuesta la funcion f (x) = 2x2 + 5x entre 1 y 2.5 con incremento
h de 0.05. Comparar con el valor exacto.
 Gua VIII: Ecuaciones Diferenciales

Ecuaciones Diferenciales

1. Resolver usando el metodo de Euler, el siguiente problema de valores iniciales :

(
dy 2
dt = y t + 1 0 t 2
y(0) = 0.5
Hacer una tabla comparando los valores en las iteraciones. Usar h = 0.2 y h = 0.1. Obtener la solucion para
t = 0.6.

2. Idem con el metodo de Runge Kutta del punto medio.

3. Idem con el metodo de Runge Kutta de orden 4.

4. Resolver usando el metodo de Runge Kutta de orden 4 el siguiente problema de valores iniciales:
(
dy 2 t
dt = (t + 1) 0.5e 0t2
y(0) = 0.5
Con los pasos del ejercicio 1.

5. Resolver usando el metodo de Euler y Runge Kutta del punto medio los siguientes problemas de valores iniciales:

a) (
dy
dt = (t + 1)3 0.4et 1 t 2
y(1) = 0.2

b) (
dy
dt = t2 0.5ln(t) 0 t 2
y(1) = 0.2

c)
10t2
(
dy
dt = t+4 0t2
y(1) = 0.2

6. Responder: Que hecho justifica el uso de metodos multipasos respecto a los de un solo paso como los hasta ahora
empleados?

7. Resolver por el metodo de Adams Bashforth el problema del punto 4.

8. Resolver por el metodo de Adams Moulton el problema del punto 4.

9. Comparar en una tabla los resultados obtenidos al resolver el 4, con los metodos usados en dicho ejercicio y en
los ejercicios 7 y 8. Sacar conclusiones.

20
Gua IX: Ecuaciones diferenciales de segundo orden y problemas de
valores de contorno

Ecuaciones diferenciales de segundo orden y problemas de valores de contorno

1. Considerar el siguiente problema de valores iniciales:

00
y (t) + sen(y(t)) = 0


y(0) = 2
y 0 (0) = 0

Realizar dos pasos del metodo de Euler con h = 0.1.

2. Considerar el siguiente PVI: 00 0

y (t) = 0.2y (t) + t

y(0) = 0
y 0 (0) = 0

Realizar dos pasos del metodo de Euler con h = 0.1.

3. Considerar el siguiente problema de valores iniciales:

00
y (t) = f (t) t [0; T ]


y(0) = 0
y 0 (0) = 0

La funcion f (t) viene dada en la siguiente tabla de valores:

tk 0.0 0.1 0.2 0.3 0.4 0.5
f (tk ) 1.0 1.0 2.0 2.0 0.0 0.0

a) Obtener el sistema de ecuaciones en diferencias que resulta cuando se utiliza el metodo de Euler para
aproximar la solucion del sistema de ecuaciones diferenciales de primer orden equivalente al problema de
valor inicial considerado.
b) Obtener la aproximacion a y(tk ) con ese metodo para k tomando los valores de 0 a 6 si se elige el paso
h = 0.1. Presentar una tabla redondeando los resultados a dos decimales pero operando en su calculadora
con la maxima precision.

4. Resolver el siguiente problema de valores de contorno utilizando diferencias finitas:

00
x (t) = x + 4et t [0; 0.5]


x(0) = 1

x(0.5) = 2e0.5

21
22

5. Se tiene el siguiente problema con valores de contorno:

0

00 (x )2
x (t) = x

t [1; 6.0]
x(1) = 4


x(6.0) = 6.0

0
a) Hallar la funcion que expresa a x(6) en funcion de x (1).

(Sugerencia: La solucion general de la ecuacion de segundo orden es x(t) = A. B + t)
0
b) Con el valor obtenido para x (1) en el tem a) hallar una aproximacion de la solucion exacta de x(t) usando
como paso h = 1.

6. Resolver el siguiente problema por el metodo del disparo, usando el paso h que crea conveniente:
00
2 0 2
x (t) + 2.(1 t )x (t) + x = 1 t [1; 1]


0
x (1) = 0


x(2) = 1

7. Sea el siguiente problema: 00

x (t) + x(t) = 0 t [0; 1]


0
x (1) = 0


x(0) = 1

1
Resolver el problema para un paso de calculo h = 2 y h = 31 .

8. Obtener una aproximacion de la curva que debe tener una cuerda para que una bolilla se desplace desde el punto
de coordenadas (1; 7) hasta el punto (0; 1) en el menor tiempo posible. Se propone plantear dicho problema fsico
en la forma del siguiente problema de valores en la frontera:
00
x (t) = x(t) t [0; 1]


x( 1) = 7


x(0) = 1