Este es un método implícito: el valor aparece en ambos lados de la ecuación, y para calcularlo realmente, hay que resolver una ecuación que será generalmente no lineal. Un método posible para resolver esta ecuación es el método de Newton. Se puede utilizar el método de Euler para obtener una estimación bastante buena para la solución, que puede utilizarse a su vez como valor inicial del método de Newton.[2]
Justificación
Integrando la ecuación diferencial de a , se tiene que
La regla del trapecio indica que la integral en el lado derecho puede ser aproximada como
Ahora, combinando ambas fórmulas y utilizando y se obtiene la regla trapezoidal para resolver ecuaciones diferenciales ordinarias.[1]
Análisis de errores
Se deduce del análisis de error de la regla trapezoidal para cuadratura que el error local de truncado de la regla trapezoidal para resolver ecuaciones diferenciales puede ser delimitado como:
Por lo tanto, la regla trapezoidal es un método de segundo orden. Este resultado se puede utilizar para demostrar que el error global es cuando el tamaño del paso tiende a cero (véase cota superior asintótica para el significado de esta afirmación).[3]
Estabilidad
En color rosa, la región de estabilidad del método trapezoidal.
La región de absoluta estabilidad para la regla trapezoidal viene dada por:
Esto incluye el semiplano izquierdo, por lo que la regla trapezoidal es A-estable. La segunda barrera de Dahlquist establece que la regla trapezoidal es la más precisa entre los métodos lineales multipaso estable. Más precisamente, un método lineal de varios pasos que es A-estable tiene como máximo el orden dos, y la constante de error de un método lineal multipaso estable de segundo orden no puede ser mejor que la constante de error de la regla trapezoidal.[2]
De hecho, la región de estabilidad absoluta para la regla trapezoidal es precisamente el semiplano izquierdo. Esto significa que si la regla trapezoidal se aplica a la ecuación de prueba lineal y = λyy, la solución numérica tiende a cero si y solo si la solución exacta también lo hace.