Ejemplo: función jacobiana aumentada para Stiffb y Stiffr
Utilice la función
Jacob para calcular la matriz jacobiana aumentada de una ecuación diferencial ordinaria (EDO) y, a continuación, proporciónela como entrada para los solvers
Stiffb y
Stiffr.
1. Defina un sistema de cuatro incógnitas:
Las variables yi son funciones de x.
2. Defina las constantes del sistema.
3. Defina los valores iniciales.
4. Defina una función de vector D(x,y) que se corresponda con el lado derecho del sistema.
5. Utilice la función
augment para crear la función jacobiana aumentada:
6. Llame a las funciones Stiffb y Stiffr:
Las matrices devueltas contienen cinco columnas que corresponden al número de puntos y las soluciones de las cuatro incógnitas.
7. Extraiga las soluciones de las cuatro incógnitas de las matrices Stiffb y Stiffr devueltas:
8. Trace y compare las soluciones devueltas para cada incógnita de las dos funciones:
Los gráficos muestran que los dos solvers de EDO devuelven soluciones idénticas.
