Ejercicios resueltos

Problemas trabajados con datos concretos y su solución detallada.

Fundamentos

Ejercicio: error iterativo con e^0.5

Aproximación de e0.5e^{0.5} por su serie de Taylor añadiendo términos hasta que el error iterativo porcentual baja del 0.05 %.

Interpolación

Ejercicio: Newton con datos de población

Tabla de diferencias divididas completa para el censo 1971–2011 y estimación de la población en 2005 con el polinomio de Newton de grado 4.

Diferenciación

Ejercicio: comparar fórmulas de la derivada

Cálculo de f(0.5)f'(0.5) para f(x)=x2exf(x)=x^2e^{-x} con las seis fórmulas de diferencias finitas y comparación de sus errores frente al valor exacto 0.45490.4549.

Ejercicio: Richardson sobre ln(x)

Aproximación de f(1.8)f'(1.8) para f(x)=lnxf(x)=\ln x con diferencias progresivas O(h)\mathcal{O}(h) y mejora a O(h2)\mathcal{O}(h^2) mediante extrapolación de Richardson.

Integración

Ejercicio: trapecio y punto medio

Comparación de trapecio compuesto, punto medio simple y punto medio compuesto, con errores absolutos y relativos.

EDO

Ejercicio: estabilidad de Euler explícito e implícito

Análisis completo de y=λyy'=\lambda y: factor de amplificación de cada método, condición de estabilidad h<2/λh<-2/\lambda del explícito, estabilidad incondicional del implícito y comprobación numérica con λ=10\lambda=-10.

Ejercicio: el modelo SIR con Euler, Heun y RK4

Integración del sistema epidémico SIR con los tres métodos de un paso sobre la misma malla, comparando cómo el orden del método cambia visiblemente los resultados.

Ejercicio: AB2 y estimación del orden

Resolución de un PVI logístico (Verhulst) con AB2 arrancado con Heun y estimación numérica del orden duplicando el número de subintervalos.

Sistemas lineales

Ejercicio: Jacobi y Gauss-Seidel comparados

Resolución de un sistema 4×4 con Jacobi y con Gauss-Seidel desde x⁰=0, comparando cuántas iteraciones necesita cada uno para la misma tolerancia.

Ejercicio: convergencia por radio espectral

Una matriz no diagonal dominante donde Jacobi diverge pero Gauss-Seidel converge, decidido calculando el radio espectral de cada matriz de iteración.

Ecuaciones no lineales

Ejercicio: bisección a mano

Seis iteraciones de bisección para cos²x−x=0 en [0,1], con la cadena de intervalos, la cota de error en cada paso y la predicción del número de iteraciones necesarias.

Ejercicio: Newton sobre x=cos²x

Aplicación completa del método de Newton a x=cos²x desde x0=0.3: tabla de iterados, residuos, incrementos y ACOC tendiendo al orden teórico 2.

Ejercicio: la secante a mano

El método de la secante aplicado a cos²x−x=0 desde x0=0, x1=1: cinco iteraciones con las diferencias divididas explícitas y la convergencia superlineal a la vista.

Ejercicio: comparativa numérica de métodos iterativos

Newton, Halley, Ostrowski, Traub, punto medio, Jarratt y Newton doble sobre funciones de prueba con tolerancia 10⁻¹⁰⁰: iteraciones, residuos y ACOC confirmando los órdenes teóricos.

Sistemas no lineales

Ejercicio: Newton para un sistema, a mano

Dos pasos de Newton a mano sobre el sistema x2+y2=1x^2+y^2=1, x=yx=y: montaje del jacobiano, resolución del sistema lineal 2×22\times2 de cada paso y convergencia cuadrática visible hacia (2/2,2/2)(\sqrt2/2,\sqrt2/2).

Ejercicio: Newton en un sistema 2×2 con tabla de iteraciones

Resolución completa de exey+xcosy=0e^x e^y+x\cos y=0, x+y=1x+y=1 con Newton desde [2,1][2,-1]: tabla de iterados, normas del residuo e incremento, y ACOC estabilizándose en 2.

Ejercicio: comparativa numérica en sistemas

Newton, Trapecios, Golden Ratio, NA, Jarratt y RN sobre dos sistemas de prueba con tolerancia 10⁻¹²: iteraciones, normas y ACOC, con RN (orden 6) como método más eficaz.