Calculo numérico

´ ´ ELEMENTOS DE CALCULO NUMERICO Ricardo G. Dur´n, Silvia B. Lassalle y Julio D. Rossi a

´ Indice General
Cap´ ıtulo 1. Punto flotante y redondeo 1. Punto flotante 2. Redondeo 3. Ejercicios Cap´ ıtulo 2. Normas y condicionamiento de una matriz. 1. Ejercicios Cap´ ıtulo 3. Resoluci´n de sistemas lineales. o 1. M´todos directos e 2. M´todos iterativos e 3. Ejercicios Cap´ ıtulo 4. Resoluci´nde ecuaciones no lineales. o 1. M´todo de bisecci´n. e o 2. M´todo regula falsi e 3. M´todo de Newton-Raphson. e 4. M´todo de punto fijo e 5. M´todo de la secante e 6. Ejercicios Cap´ ıtulo 5. Interpolaci´n o 1. Interpolaci´n de Lagrange o 2. Error de interpolaci´n o 3. Forma de Newton 4. Polinomios de Tchebychev - Minimizaci´n del Error o 5. Interpolaci´n de Hermite o 6. Interpolaci´n porpolinomios a trozos o 7. Ejercicios Cap´ ıtulo 6. Polinomios ortogonales y aproximaci´n por cuadrados m´ o ınimos. 1. Preliminares 2. Soluci´n de los Problemas de Aproximaci´n o o 3. Ejercicios Cap´ ıtulo 7. Integraci´n num´rica o e 1. F´rmulas de Newton-Cˆtes o o 2. Estimaci´n del error o 3. F´rmulas de cuadratura compuestas o
3

1 1 3 11 15 29 35 35 40 60 65 65 68 69 76 79 83 87 87 90 92 95 102 104107 111 112 118 127 131 132 139 143

4

´ INDICE GENERAL

4. 5. 6.

Convergencia de los m´todos de cuadratura e Cuadratura Gaussiana Ejercicios

148 150 154 159 162 164 166 170 176 177

Cap´ ıtulo 8. Resoluci´n de ecuaciones diferenciales ordinarias. o 1. M´todos de Euler y Taylor de orden k e 2. M´todos de Runge-Kutta e 3. An´lisis de los Errores a 4. M´todos multipaso lineales e 5.M´todos de paso variable e 6. Ejercicios

CAP´ ıTULO 1

Punto flotante y redondeo
El objeto de este cap´ ıtulo es analizar la representaci´n de los n´meros en una computadora y la o u propagaci´n de los errores de redondeo al realizar c´lculos. o a Como la cantidad de informaci´n que puede guardarse en una computadora es finita, la m´quina o a trabajar´ s´lo con un conjunto finito de n´meros. A´stos los llamaremos n´meros de m´quina. a o u e u a En consecuencia, toda vez que de nuestros datos o c´lculos surja un n´mero que no pertenece a a u este conjunto finito, ´ste deber´ ser reemplazado por una aproximaci´n (el n´mero de m´quina e a o u a m´s cercano). Este reemplazo da lugar a lo que llamamos errores de redondeo. a Al realizar c´lculos estos errores de redondeo se propagan y esto puedellevar a resultados a totalmente incorrectos como veremos en algunos ejemplos simples. En las aplicaciones del c´lculo num´rico es pr´cticamente imposible determinar exactamente la a e a magnitud de los errores de redondeo. Lo que si puede hacerse, y es de fundamental importancia, es identificar las posibles causas de que los errores se propaguen m´s de lo admisible. Esto a permite mejorar losalgoritmos o determinar que m´todo es m´s conveniente para resolver un e a problema. Un claro ejemplo de esto, que veremos m´s adelante, aparece cuando se utiliza el a m´todo de eliminaci´n de Gauss para resolver un sistema de ecuaciones lineales. En este caso, e o el an´lisis de la propagaci´n de errores permite determinar la forma m´s eficiente de aplicar el a o a m´todo. e Por otra parte, esfundamental distinguir cuando la propagaci´n excesiva de errores se debe a o que el algoritmo utilizado es “malo” o inestable o a que el problema en s´ mismo est´ “mal ı a condicionado”. En el primer caso se puede (se debe!) tratar de mejorar el m´todo de resoluci´n e o mientras que en el segundo caso el problema es m´s esencial. Los ejemplos que presentaremos a ilustrar´n estos dos casos. a

1. Puntoflotante En lo que sigue supondremos que los n´meros de m´quina son los que aparecen en la pantalla. u a Esto no es exacto pues en realidad la computadora opera internamente con los n´meros desarrou llados en base 2 y no en base 10. Este abuso de lenguaje es s´lo para mayor claridad (el lector o podr´ observar que todo nuestro an´lisis puede repetirse trabajando en base 2). a a Observemos...