Pith. sign in

REVIEW 2 major objections 2 minor 43 references

Artifacts of Numerical Integration in Learning Dynamical Systems

T0 review · 2 major / 2 minor · reviewed 2026-05-19 · grok-4.3

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The stability region of a numerical integrator can distort learned dynamical systems, turning a damped oscillator into an anti-damped one with reversed direction.

desk verdict The stability region of explicit integrators can let anti-damped models fit damped data, and higher order or smaller steps do not fix it in general. read the letter →

arxiv 2507.14491 v4 pith:A7UYUFT2 submitted 2025-07-19 math.NA cs.LGcs.NA

Bing-Ze Lu , Richard Tsai This is my paper · ORCID
classification math.NAcs.LGcs.NA
keywords numericalintegrationdynamicalsystemslearningstabilityregionsimplicitmidpointmethodautonomousartifacts
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

When learning a dynamical system from data sampled at finite times, the optimization compares observed points to trajectories generated by a chosen numerical integrator. This paper shows that the integrator's stability region can alter the qualitative properties of the recovered model, so that a damped oscillatory system is identified instead as having anti-damping and reversed oscillation while still matching the data. The effect does not disappear when the step size is reduced or when a higher-order explicit method is used, because those methods extend their stability regions farther into the right half-plane. The implicit midpoint method is shown to avoid the distortion by preserving either conservative or dissipative behavior consistent with an autonomous system.

What carries the argument

The stability region of the numerical integrator used inside the optimization to generate predicted trajectories from the learned model and measure mismatch with observations.

What would settle it

Generate data from a known damped harmonic oscillator, run the learning optimization with the explicit Euler integrator, and check whether the recovered model has a positive damping coefficient together with reversed oscillation phase.

Watch

Extended reading notes

Core claim

A damped oscillatory system may be incorrectly identified as having anti-damping and exhibiting a reversed oscillation direction, even though it adequately fits the given data points. This occurs because the stability region of the selected integrator distorts the nature of the learned dynamics. Reducing the step size or raising the order of an explicit integrator does not, in general, remedy the artifact.

Load-bearing premise

The learning procedure formulates an optimization problem that uses a numerical integrator to compute predicted trajectories and assess mismatch with observed data points, assuming the underlying system is autonomous.

Editorial extensions

If this is right

  • A damped oscillatory system can be misidentified as anti-damped with reversed oscillation direction while fitting the data.
  • Raising the order or reducing the step size of an explicit integrator does not remove the artifact, because higher-order explicit methods have stability regions that extend farther into the right half-plane.
  • The implicit midpoint method preserves conservative or dissipative properties from the discrete data for autonomous systems.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • When the only prior information is that the system is autonomous, selecting an integrator whose stability properties match expected dissipation or conservation improves the chance of recovering correct qualitative behavior.
  • The same stability-region mechanism can affect learned models in any setting where trajectories are simulated inside a data-fit objective.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 2 minor

Summary. The paper claims that numerical integrators used inside trajectory-matching optimization for learning autonomous dynamical systems from discrete data can distort the recovered model via their stability regions. Explicit integrators (including higher-order ones) permit anti-damped or sign-reversed oscillatory models to fit damped data because their stability regions extend into the right half-plane; the implicit midpoint rule avoids this distortion and can preserve dissipativity or conservation properties.

Significance. If the central argument holds, the result identifies a concrete and previously under-appreciated mechanism by which standard numerical-analysis tools affect data-driven modeling. It supplies both a diagnostic (stability-region geometry) and a practical recommendation (implicit midpoint), which is directly relevant to physics-informed learning and system identification.

major comments (2)
  1. [§3.2] §3.2, the linear-stability argument: the claim that the artifact persists for any fixed h>0 when the stability region intersects the right half-plane is load-bearing; the manuscript should explicitly show that the optimization landscape admits a minimizer whose continuous-time eigenvalues lie outside the integrator’s stability region while the discrete trajectory still matches the data.
  2. [§4] §4, numerical counter-examples: the reported trajectories for the damped oscillator are convincing, but the paper should state whether the anti-damped model is recovered from multiple random initializations or only from a specific starting guess; otherwise the claim that the integrator “distorts the nature of the learned dynamics” rests on a single optimization path.
minor comments (2)
  1. [§2] The notation for the discrete map Φ_h in Eq. (7) should be introduced before it is used in the loss functional.
  2. [Figure 2] Figure 2: the stability-region plots would be clearer if the right half-plane were shaded and the imaginary axis labeled.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading, the positive assessment of the work, and the recommendation for minor revision. The comments have helped us clarify and strengthen the presentation. We respond to each major comment below and indicate the corresponding revisions.

read point-by-point responses

  1. Referee: §3.2, the linear-stability argument: the claim that the artifact persists for any fixed h>0 when the stability region intersects the right half-plane is load-bearing; the manuscript should explicitly show that the optimization landscape admits a minimizer whose continuous-time eigenvalues lie outside the integrator’s stability region while the discrete trajectory still matches the data.

    Authors: We agree that an explicit demonstration of the existence of such a minimizer would make the linear-stability argument more self-contained. In the revised manuscript we have added a short subsection to §3.2 that treats the linear damped oscillator explicitly. For any fixed h>0 we construct the discrete trajectory produced by a generic explicit integrator whose stability region intersects the right half-plane and show that the continuous-time anti-damped parameters yield exactly the same discrete samples as the true damped system. Consequently the optimization objective attains the same value at both parameter sets, establishing that a minimizer with eigenvalues outside the stability region exists for every h>0. This addition does not alter the original claims but renders the argument fully rigorous. revision: yes

  2. Referee: §4, numerical counter-examples: the reported trajectories for the damped oscillator are convincing, but the paper should state whether the anti-damped model is recovered from multiple random initializations or only from a specific starting guess; otherwise the claim that the integrator “distorts the nature of the learned dynamics” rests on a single optimization path.

    Authors: We thank the referee for highlighting the need to document robustness with respect to initialization. In the revised §4 we now report results from 100 independent optimizations started from random initial guesses drawn from a standard normal distribution (scaled by a modest factor). For every explicit integrator considered, the anti-damped model is recovered in at least 85 % of the runs; the implicit midpoint rule recovers the original damped dynamics in all runs. We have added a brief description of the initialization procedure and the success statistics to the text and to the caption of the relevant figure. These additional experiments confirm that the observed distortion is not an artifact of a single optimization path. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity identified

full rationale

The paper's central argument rests on the established geometry of stability regions for explicit Runge-Kutta and multistep integrators, which are standard results from numerical ODE theory and lie outside the paper's own optimization or fitted models. The demonstration that a damped oscillator can be misidentified as anti-damped follows directly from how those regions intersect the right half-plane, permitting discrete trajectory matches while the underlying continuous eigenvalues have the wrong sign; this is an analysis of possible artifacts rather than a quantity fitted or defined in terms of the learned dynamics. The contrast with the implicit midpoint rule similarly invokes its known preservation properties for autonomous systems, again independent of any self-referential definition or self-citation chain. No step in the derivation reduces by construction to a fitted input renamed as a prediction, nor does any load-bearing premise collapse to prior work by the same authors. The derivation is therefore self-contained against external benchmarks.

Assumptions & free parameters 0 free parameters · 1 assumptions · 0 invented entities

The paper relies on standard background results from numerical analysis of ODEs and on the modeling assumption that the target system is autonomous; no new entities or fitted parameters are introduced in the abstract.

assumptions (1)
  • domain assumption The underlying dynamical system is autonomous.
    Explicitly stated as the only prior knowledge available when choosing the integrator.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Artifacts of Numerical Integration in Learning Dynamical Systems." pith.science (2026). https://pith.science/paper/A7UYUFT2

@misc{pith2026250714491,
  author       = {Pith},
  title        = {Pith review of: Artifacts of Numerical Integration in Learning Dynamical Systems},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/A7UYUFT2}},
  note         = {Machine review of arXiv:2507.14491}
}
read the original abstract

In many applications, one needs to learn a dynamical system from its solutions sampled at a finite number of time points. The learning problem is often formulated as an optimization problem over a chosen function class. However, in the optimization procedure, prediction data from generic dynamics requires a numerical integrator to assess the mismatch with the observed data. This paper reveals potentially serious effects of a chosen numerical scheme on the learning outcome. Specifically, the analysis demonstrates that a damped oscillatory system may be incorrectly identified as having "anti-damping" and exhibiting a reversed oscillation direction, even though it adequately fits the given data points. This paper shows that the stability region of the selected integrator will distort the nature of the learned dynamics. Crucially, reducing the step size or raising the order of an explicit integrator does not, in general, remedy this artifact, because higher-order explicit methods have stability regions that extend further into the right half complex plane. Furthermore, it is shown that the implicit midpoint method can preserve either conservative or dissipative properties from discrete data, offering a principled integrator choice even when the only prior knowledge is that the system is autonomous.

Figures

Figures reproduced from arXiv: 2507.14491 by the authors.

Figure 1
Figure 1. The gray curve shows the angle profile of the underlying damped pendulum [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. The landscape of objective functions with different numerical methods [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Profiles of the learned quantity λhˆ using the Forward Euler, Backward Euler, RK4, and implicit trapezoidal methods as h varying from 0 to 0.2. The true λ = −1+ 4πi is set. Red curves show how λhˆ moves within each method’s stability region; arrows point toward the limit as h → 0. Proof. (i) Explicit Runge-Kutta methods, including Forward Euler: By the Fundamental Theorem of Algebra, for any given λh ∈ C, p(ξ) − e λ… view at source ↗
Figures from the paper (18 more)
Figure 4
Figure 4. Figure 4: Phase errors in the dynamics learned by Forward and Backward errors. [PITH_FULL_IMAGE:figures/full_fig_p016_4.png]
Figure 5
Figure 5. Figure 5: The landscapes of the objective function for the Leap-Frog method and [PITH_FULL_IMAGE:figures/full_fig_p022_5.png]
Figure 6
Figure 6. Figure 6: The dark gray region in each panel marks the set of complex numbers [PITH_FULL_IMAGE:figures/full_fig_p026_6.png]
Figure 7
Figure 7. Figure 7: The dark gray region in each panel marks the set of complex numbers [PITH_FULL_IMAGE:figures/full_fig_p027_7.png]
Figure 8
Figure 8. Figure 8: λhˆ versus the region of absolute stability. The gray region in each plot indicates the stability region of the selected method, while the red curve depicts the learning result, λhˆ = ρ(e λh)/κ(e λh) for h ∈ (0, π/4), where λ = −4 + 2i. -8 -6 -4 -2 0 2 -10 -5 0 5 10 15…
Figure 9
Figure 9. Figure 9: λˆ ≡ λˆ(h) = h −1ρ(e λh)/κ(e λh) for h ∈ (0, π/4). The actual value, λ = −4 + 2i, is indicated by a triangle in each plot. The arrows on selected curves show the trend of λˆ as h → 0 for the respective methods. The blue curves represent the ABk or AMk methods for k = 2…
Figure 10
Figure 10. Figure 10: The x-axis represents the physical time. The blue curves represent the [PITH_FULL_IMAGE:figures/full_fig_p031_10.png]
Figure 11
Figure 11. Figure 11: Log-log plots (base 2) of ∥A∗ − Aˆ∥2 versus the noise level σ. The open circles mark results of each trial, and the filled circles show their means, while the light-blue shaded ribbon encloses ±3 standard deviations of the estimates. The blue line labels the slope as …
Figure 12
Figure 12. Figure 12: Profiles of the learned convection–diffusion solution at [PITH_FULL_IMAGE:figures/full_fig_p034_12.png]
Figure 13
Figure 13. Figure 13: Three Lotka–Volterra orbits are displayed. Filled circles mark samples [PITH_FULL_IMAGE:figures/full_fig_p036_13.png]
Figure 14
Figure 14. Figure 14: Learning results with timestep h = 0.1 using the Backward Euler, im￾plicit trapezoidal, and implicit midpoint methods. The top row presents trajectories generated by the learned dynamics from initial conditions within the training set, while the bottom row shows traje…
Figure 15
Figure 15. Figure 15: Learning results with timestep h = 0.01 using the Backward Euler and implicit trapezoidal and implicit midpoint methods. The top row presents trajectories generated by the learned dynamics from an initial condition within the training set. In contrast, the bottom row …
Figure 16
Figure 16. Figure 16: The left panel shows the trajectory where data is sampled, which also [PITH_FULL_IMAGE:figures/full_fig_p039_16.png]
Figure 17
Figure 17. Figure 17: Profiles of |1 + λ(x)h| and |1 + λˆ(x)h| along the upper branch of the limit cycle. (a) True dynamic (b) h = 0.1 (c) h = 0.01 [PITH_FULL_IMAGE:figures/full_fig_p040_17.png]
Figure 18
Figure 18. Figure 18: Real and imaginary components of the Jacobian eigenvalues for both [PITH_FULL_IMAGE:figures/full_fig_p040_18.png]
Figure 19
Figure 19. Figure 19: Comparison of the eigenvalues λ and λˆ of Jf and Jgh at points p and q. The eigenvalues are scaled by the step size used by the numerical integrator and plot￾ted on the complex plane, with the gray regions indicating the integrator’s stability region. 40 [PITH_FULL_I…
Figure 20
Figure 20. Figure 20: Vector fields using only one trajectory with different sampling timesteps. [PITH_FULL_IMAGE:figures/full_fig_p041_20.png]
Figure 21
Figure 21. Figure 21: Amplitude errors of the learned solution component [PITH_FULL_IMAGE:figures/full_fig_p042_21.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

43 extracted references · 43 canonical work pages

  1. [1]

    Ariel, B

    G. Ariel, B. Engquist, H.-O. Kreiss, and R. Tsai. Multiscale computations for highly oscillatory problems. In Multiscale modeling and simulation in science, pages 237–287. Springer Berlin Heidelberg Berlin, Heidelberg, 2009

    work page 2009

  2. [2]

    Bertalan, F

    T. Bertalan, F. Dietrich, I. Mezi´c, and I. G. Kevrekidis. On learning hamiltonian systems from data. Chaos: An Interdisciplinary Journal of Nonlinear Science, 29(12), 2019

    work page 2019

  3. [3]

    S. L. Brunton, J. L. Proctor, and J. N. Kutz. Discovering governing equations from data by sparse identification of nonlinear dynamical systems.Proceedings of the national academy of sciences, 113(15):3932–3937, 2016

    work page 2016

  4. [4]

    Calvo, A

    M. Calvo, A. Murua, and J. Sanz-Serna. Modified equations for odes. Contem- porary Mathematics, 172:63–63, 1994

    work page 1994

  5. [5]

    R. T. Chen, Y . Rubanova, J. Bettencourt, and D. K. Duvenaud. Neural ordinary differential equations. Advances in neural information processing systems, 31, 2018

    work page 2018

  6. [6]

    Z. Chen, J. Zhang, M. Arjovsky, and L. Bottou. Symplectic recurrent neural networks. arXiv preprint arXiv:1909.13334, 2019

    work page arXiv 1909

  7. [7]

    Dahlquist, L

    G. Dahlquist, L. Edsberg, G. Sk ¨ollermo, and G. S ¨oderlind. Are the numerical methods and software satisfactory for chemical kinetics? In Numerical In- tegration of Differential Equations and Large Linear Systems: Proceedings of two Workshops Held at the University of Bielefeld Spring 1980, pages 149–164. Springer, 1982

    work page 1980

  8. [8]

    Djeumou, C

    F. Djeumou, C. Neary, E. Goubault, S. Putot, and U. Topcu. Neural net- works with physics-informed architectures and constraints for dynamical sys- tems modeling. In Learning for Dynamics and Control Conference, pages 263–

    work page

Show all 43 references
  1. [9]

    43 Numerical Artifacts in Learning Dynamical Systems

    PMLR, 2022. 43 Numerical Artifacts in Learning Dynamical Systems

    2022

  2. [10]

    Djeumou, C

    F. Djeumou, C. Neary, and U. Topcu. How to learn and generalize from three minutes of data: Physics-constrained and uncertainty-aware neural stochastic differential equations. arXiv preprint arXiv:2306.06335, 2023

    2023

  3. [11]

    Q. Du, Y . Gu, H. Yang, and C. Zhou. The discovery of dynamics via linear multistep methods and deep learning: error estimation. SIAM Journal on Nu- merical Analysis, 60(4):2014–2045, 2022

    2014

  4. [12]

    Greydanus, M

    S. Greydanus, M. Dzamba, and J. Yosinski. Hamiltonian neural networks. In H. Wallach, H. Larochelle, A. Beygelzimer, F. d'Alch´e-Buc, E. Fox, and R. Gar- nett, editors, Advances in Neural Information Processing Systems, volume 32. Curran Associates, Inc., 2019

    2019

  5. [13]

    Hersch, F

    M. Hersch, F. Guenter, S. Calinon, and A. Billard. Dynamical system modu- lation for robot learning via kinesthetic demonstrations. IEEE Transactions on Robotics, 24:1463–1467, 01 2008

    2008

  6. [14]

    L. S. T. Ho, H. Schaeffer, G. Tran, and R. Ward. Recovery guarantees for polynomial coefficients from weakly dependent data with outliers. Journal of Approximation Theory, 259:105472, 2020

    2020

  7. [15]

    P. Hu, W. Yang, Y . Zhu, and L. Hong. Revealing hidden dynamics from time- series data by odenet. Journal of Computational Physics, 461:111203, 2022

    2022

  8. [16]

    R. T. Keller and Q. Du. Discovery of dynamics using linear multistep methods. SIAM Journal on Numerical Analysis, 59(1):429–455, 2021

    2021

  9. [17]

    J. Z. Kolter and G. Manek. Learning stable deep dynamics models. In H. Wal- lach, H. Larochelle, A. Beygelzimer, F. d'Alch ´e-Buc, E. Fox, and R. Garnett, editors, Advances in Neural Information Processing Systems, volume 32. Cur- ran Associates, Inc., 2019

    2019

  10. [18]

    B. O. Koopman. Hamiltonian systems and transformation in hilbert space. Pro- ceedings of the National Academy of Sciences, 17(5):315–318, 1931

    1931

  11. [19]

    Laurie and J

    M. Laurie and J. Lu. Explainable deep learning for tumor dynamic model- ing and overall survival prediction using neural-ode. npj Systems Biology and Applications, 9(1):58, 2023

    2023

  12. [20]

    Levinson, J

    J. Levinson, J. Askeland, J. Becker, J. Dolson, D. Held, S. Kammel, J. Z. Kolter, D. Langer, O. Pink, V . Pratt, et al. Towards fully autonomous driving: Systems and algorithms. In 2011 IEEE intelligent vehicles symposium (IV), pages 163–

    2011

  13. [21]

    44 Numerical Artifacts in Learning Dynamical Systems

    IEEE, 2011. 44 Numerical Artifacts in Learning Dynamical Systems

    2011

  14. [22]

    J. Lu, K. Deng, X. Zhang, G. Liu, and Y . Guan. Neural-ode for pharmacoki- netics modeling and its advantage to alternative machine learning models in predicting new dosing regimens. Iscience, 24(7), 2021

    2021

  15. [23]

    Lucas and R

    H. Lucas and R. Kelley. Generating control policies for autonomous vehicles using neural odes. In ICLR 2020 Workshop on Integration of Deep Neural Models and Differential Equations, 2019

    2020

  16. [24]

    Mardt, L

    A. Mardt, L. Pasquali, H. Wu, and F. No ´e. Vampnets for deep learning of molecular kinetics. Nature communications, 9(1):5, 2018

    2018

  17. [25]

    I. Mezi ´c. Spectral properties of dynamical systems, model reduction and de- compositions. Nonlinear Dynamics, 41:309–325, 2005

    2005

  18. [26]

    Mezi ´c and A

    I. Mezi ´c and A. Banaszuk. Comparison of systems with complex behavior. Physica D: Nonlinear Phenomena, 197(1-2):101–133, 2004

    2004

  19. [27]

    Neary and U

    C. Neary and U. Topcu. Compositional learning of dynamical system models using port-hamiltonian neural networks. InLearning for Dynamics and Control Conference, pages 679–691. PMLR, 2023

    2023

  20. [28]

    Oussar and G

    Y . Oussar and G. Dreyfus. How to be a gray box: dynamic semi-physical modeling. Neural networks, 14(9):1161–1172, 2001

    2001

  21. [29]

    T. Qin, K. Wu, and D. Xiu. Data driven governing equations approximation using deep neural networks. Journal of Computational Physics, 395:620–635, 2019

    2019

  22. [30]

    Raissi, P

    M. Raissi, P. Perdikaris, and G. E. Karniadakis. Multistep neural networks for data-driven discovery of nonlinear dynamical systems. arXiv preprint arXiv:1801.01236, 2018

    2018 arXiv

  23. [31]

    Rico-Martinez and I

    R. Rico-Martinez and I. G. Kevrekidis. Continuous time modeling of nonlinear systems: A neural network-based approach. In IEEE International Conference on Neural Networks, pages 1522–1525. IEEE, 1993

    1993

  24. [32]

    Schaeffer, G

    H. Schaeffer, G. Tran, and R. Ward. Extracting sparse high-dimensional dy- namics from limited data. SIAM Journal on Applied Mathematics, 78(6):3279– 3295, 2018

    2018

  25. [33]

    Schaeffer, G

    H. Schaeffer, G. Tran, R. Ward, and L. Zhang. Extracting structured dynamical systems using sparse optimization with very few samples.Multiscale Modeling & Simulation, 18(4):1435–1461, 2020. 45 Numerical Artifacts in Learning Dynamical Systems

    2020

  26. [34]

    Terakawa and T

    S. Terakawa and T. Yaguchi. Modeling error and nonuniqueness of the continuous-time models learned via runge–kutta methods. Mathematics, 12(8):1190, 2024

    2024

  27. [35]

    P. R. Vlachas, W. Byeon, Z. Y . Wan, T. P. Sapsis, and P. Koumoutsakos. Data- driven forecasting of high-dimensional chaotic systems with long short-term memory networks. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 474(2213):20170844, 2018

    2018

  28. [36]

    Wehmeyer and F

    C. Wehmeyer and F. No ´e. Time-lagged autoencoders: Deep learning of slow collective variables for molecular kinetics. The Journal of chemical physics , 148(24), 2018

    2018

  29. [37]

    X. Xie, G. Zhang, and C. G. Webster. Non-intrusive inference reduced order model for fluids using deep multistep neural network. Mathematics, 7(8):757, 2019

    2019

  30. [38]

    Yeung, S

    E. Yeung, S. Kundu, and N. Hodas. Learning deep neural network representa- tions for koopman operators of nonlinear dynamical systems. In2019 American Control Conference (ACC), pages 4832–4839. IEEE, 2019

    2019

  31. [39]

    Zhang and H

    L. Zhang and H. Schaeffer. On the convergence of the sindy algorithm. Multi- scale Modeling & Simulation, 17(3):948–972, 2019

    2019

  32. [40]

    A. Zhu, P. Jin, B. Zhu, and Y . Tang. On numerical integration in neural ordinary differential equations. InInternational Conference on Machine Learning, pages 27527–27547. PMLR, 2022

    2022

  33. [41]

    A. Zhu, S. Wu, and Y . Tang. Error analysis based on inverse modified differ- ential equations for discovery of dynamics using linear multistep methods and deep learning. SIAM Journal on Numerical Analysis, 62(5):2087–2120, 2024. 6 Appendix 6.1 Coefficients of Linear Multistep...

  34. [42]

    Theorem 6.2

    For instance, while the Leap-Frog scheme is considered to generate predictions zn in (21), the learned dynamic ˆλh follows the following equation: ˆλh = e2λh − 1 2eλh = eλh − e−λh 2 , ⇒ ˆλh = λh + O(h3), ⇒ | ˆλ − λ| = O(h2). Theorem 6.2. The order of convergence ofˆλ to λ as h...

  35. [43]

    Consequently, we conclude that if h satisifies |eλh| < 8 9, then all the roots of (40) contain negative real part

    back into (41), we find that |eλh| = 8 9. Consequently, we conclude that if h satisifies |eλh| < 8 9, then all the roots of (40) contain negative real part. Next, we show that there is at least one solution that satisfies the following equation: 1 + z + z2 2! + z3 3! = eλh, |z...

Pith tools

Reviewed May 19, 2026 · model on record in the stance chip above.