Deducción: Richardson sube el orden

De la forma del error a la fórmula de extrapolación: por qué combinar N1(h)N_1(h) y N1(h/2)N_1(h/2) elimina el término O(h)\mathcal{O}(h).

Eliminar el término O(h)

  1. Partimos del error con todas las potencias y evaluamos también con paso h/2h/2:

    M=N1(h)+k1h+k2h2+M=N1 ⁣(h2)+k1h2+k2h24+\begin{aligned} M&=N_1(h)+k_1 h+k_2 h^2+\cdots\\ M&=N_1\!\left(\tfrac{h}{2}\right)+k_1\tfrac{h}{2}+k_2\tfrac{h^2}{4}+\cdots \end{aligned}
  2. Hacemos 2(segunda)(primera)2\cdot(\text{segunda})-(\text{primera}). El término k1hk_1h se cancela y queda un error O(h2)\mathcal{O}(h^2):

    M=N1 ⁣(h2)+[N1 ⁣(h2)N1(h)]k2h22M=N_1\!\left(\tfrac{h}{2}\right)+\left[N_1\!\left(\tfrac{h}{2}\right)-N_1(h)\right]-k_2\tfrac{h^2}{2}-\cdots
  3. Llamamos N2(h)N_2(h) a la parte conocida: es una aproximación de orden 2. Repetir con los pesos 13\tfrac13, 115\tfrac{1}{15}, … sube más el orden.

    N2(h)=N1 ⁣(h2)+[N1 ⁣(h2)N1(h)]N_2(h)=N_1\!\left(\tfrac{h}{2}\right)+\left[N_1\!\left(\tfrac{h}{2}\right)-N_1(h)\right]