En este capítulo vamos a recordar un conjunto de preliminares matemáticos con el cálculo matemático que vamos a utilizar durante el curso.
Un conocimiento sólido de cálculo es esencial para comprender el análisis de las técnicas numéricas.
Límite de una función en un punto
El conjunto formado por las funciones \(f(x)\) que son continuas en un determinado dominio \(D\) se denota por \({\cal C}(D)\).
Ejemplo
Consideremos la sucesión \(x_n=1+\frac{1}{n}\). Tenemos que \(\displaystyle\lim_{n\to\infty}1+\frac{1}{n}=1\). Veámoslo:
sea \(\epsilon >0\), hemos que hallar un valor \(n_0\) tal que si \(n\geq n_0\), \(|x_n-1|=\left|1+\frac{1}{n}-1\right|=\frac{1}{n}<\epsilon\).
Consideremos como valor \(n_0\) el primer entero que supera \(\frac{1}{\epsilon}\), es decir \(n_0 >\frac{1}{\epsilon}\).
Así, si \(n\geq n_0\), se cumplirá que \(\frac{1}{n}\leq \frac{1}{n_0}<\epsilon\), tal como queríamos ver.
Ejemplo
Veamos usando el teorema anterior que la función \(f(x)=\sin\left(\frac{1}{x}\right)\), si \(x\neq 0\) y \(f(0)=0\), no es continua en \(x_0=0\).
Para verlo, construiremos dos sucesiones \(x_n\) e \(y_n\), ambas con \(\displaystyle\lim_{n\to \infty}x_n =\lim_{n\to\infty}y_n=0\), pero \(\displaystyle\lim_{n\to\infty}f(x_n)\neq \lim_{n\to\infty}f(y_n)\). Como debe verificarse que para toda sucesión \(z_n\) con \(\displaystyle\lim_{n\to \infty}z_n =0\), \(\displaystyle\lim_{n\to\infty}f(z_n)=f(0)\), llegamos a una contradicción.
Sea \(x_n =\frac{1}{\pi n}\). La sucesión \(x_n\) cumple \(\displaystyle\lim_{n\to \infty}x_n=0\) pero \(\displaystyle\lim_{n\to \infty} f(x_n)=\lim_{n\to \infty} \sin\left(\frac{1}{x_n}\right)=\lim_{n\to \infty} \sin\left(\frac{1}{\frac{1}{\pi n}}\right)=\lim_{n\to \infty}\sin (\pi n)=0.\)
Sea ahora \(y_n=\frac{1}{\frac{\pi}{2}+2\pi n}\). La sucesión \(y_n\) cumple \(\displaystyle\lim_{n\to \infty}y_n=0\) pero \(\displaystyle\lim_{n\to \infty} f(y_n)=\lim_{n\to \infty} \sin\left(\frac{1}{y_n}\right)=\lim_{n\to \infty} \sin\left(\frac{1}{\frac{1}{\frac{\pi}{2}+2\pi n}}\right)=\lim_{n\to \infty}\sin \left(\frac{\pi}{2}+2\pi n\right)=1.\)
Para definir funciones suaves necesitamos el concepto de diferenciabilidad:
Derivada de una función en un punto
El conjunto formado por las funciones que son \(n\) veces derivables en un dominio \(D\subseteq\mathbb{R}\) se denota por \({\cal C}^n (D)\).
El conjunto formado por las funciones que admiten cualquier derivada en un dominio \(D\subseteq\mathbb{R}\) se denota por \({\cal C}^\infty (D)\).
Los teoremas que vienen a continuación son claves para estimar el error cometido al evaluar funciones reales de variable real.
Teorema de Rolle
Ejemplo
Consideremos la función \(f(x)=\sin x\) definida en el intervalo \([0,\pi]\).
Tenemos que \(f(0)=f(\pi)=0\).
Por tanto existe un punto \(c\in (0,\pi)\) tal que \(f'(c)=\cos c=0\).
Este punto \(c\) vale \(c=\frac{\pi}{2}\) ya que \(\cos\left(\frac{\pi}{2}\right)=0\).
Teorema del valor medio
Ejemplo
Consideremos la función \(f(x)=x^2 \sin x\) definida en el intervalo \(\left[-\frac{\pi}{2},\frac{\pi}{2}\right]\).
Por tanto existe un punto \(c\in \left(-\frac{\pi}{2},\frac{\pi}{2}\right)\) tal que \(f'(c)=\frac{f\left(\frac{\pi}{2}\right)-f\left(-\frac{\pi}{2}\right)}{\frac{\pi}{2}-\left(-\frac{\pi}{2}\right)}=\frac{\frac{\pi^2}{4}+\frac{\pi^2}{4}}{\pi}=\frac{\frac{\pi^2}{2}}{\pi}=\frac{\pi}{2}\).
Usando que \(f'(c)=2c\sin c+c^2\cos c=\frac{\pi}{2}\), el valor de \(c\) será un cero de la ecuación \(2c\sin c+c^2\cos c-\frac{\pi}{2}=0\).
Durante el curso veremos técnicas para hallar soluciones aproximadas de ecuaciones de este tipo. Hay dos valores de \(c\) en el intervalo \(\left(-\frac{\pi}{2},\frac{\pi}{2}\right)\) que cumplen la condición anterior y valen aproximadamente \(c\approx\pm 0.792863\).
Teorema de los valores extremos
El siguiente teorema se demuestra aplicando sucesivamente el Teorema de Rolle a las funciones \(f,\ f',\ldots,f^{(n-1)}\):
Demostración.
Para demostrar el Teorema de Rolle generalizado hay que aplicar el Teorema de Rolle a las funciones \(f,\ f',\ldots, f^{(n-1)}\) en el sentido siguiente:
Ejemplo
Consideremos la función \(f(x)=(x-x_0) (x-x_1) (x-x_2) (x-x_3)\), con \(x_0<x_1<x_2<x_3\) definida en el intervalo \(\left[x_0,x_3\right]\).
Como \(f(x_0)=f(x_1)=f(x_2)=f(x_3)=0\), existirá un punto \(c\in (x_0,x_3)\) tal que \(f'''(c)=0\). Veamos cuál es este punto \(c\).
El valor de \(f'''(x)\) vale: \[ f'''(x)=24 x-6(x_0+x_1+x_2+x_3). \] Haciendo \(f'''(c)=0\), tenemos que el valor \(c\) vale: \(\displaystyle c=\frac{x_0+x_1+x_2+x_3}{4}\).
El Teorema del valor intermedio nos dice que si una función es continua, dicha función alcanza todos los valores entre \(f(a)\) y \(f(b)\):
Ejemplo
Consideremos la función anterior \(f(x)=x^2 \sin x\) definida en el intervalo \(\left[-\frac{\pi}{2},\frac{\pi}{2}\right]\).
Tenemos que \(f\left(-\frac{\pi}{2}\right)=-\frac{\pi^2}{4}\), \(f\left(\frac{\pi}{2}\right)=\frac{\pi^2}{4}\).
El Teorema del valor intermedio nos dice que para cualquier \(k\in\left(-\frac{\pi^2}{4},\frac{\pi^2}{4}\right)\), existe un valor \(c\in \left(-\frac{\pi}{2},\frac{\pi}{2}\right)\) tal que \(f(c)=c^2\sin c=k\).
Por ejemplo, si \(k=0\), el valor de \(c\) sería \(c=0\); si \(k=\frac{\pi}{2}\in \left(-\frac{\pi^2}{4},\frac{\pi^2}{4}\right)\), existe un valor \(c\in\left(-\frac{\pi}{2},\frac{\pi}{2}\right)\) tal que \(c^2\sin c=\frac{\pi}{2}\). Dicho valor se tiene que hallar usando métodos numéricos, concretamente usando técnicas de cálculo aproximado de ceros de funciones como veremos más adelante.
La integral de Riemann se define básicamente para formalizar el cálculo de áreas de funciones: