Hi I’m still in the newbie category so I dare ask this question as a “First Steps” question. I read this blog: https://www.johndcook.com/blog/2020/02/08/arenstorf-orbit/ … about an equation that was used to construct orbits for the Apollo program. It looked interesting, and not too complex, …

[image] rmz: x′′ = x + 2y′ + μ′*(x + μ)/D1 - μ*(x - μ′)/D2 y′′ = y + 2x′ + μ′*y/D1 - μ*y/D2 Looking at the link, it seems like some of the signs are off here; it should be x′′ = x + 2y′ - μ′*(x + μ)/D1 - μ*(x - μ′)/D2 y′′ = y - 2x′ - μ′*y/D1 - μ*y/D2 With that fix, I don’t get insta…

Second order differential equation blows up numerically, and I have no clue why. Any ideas?

New to Julia

ChrisRackauckas February 10, 2020, 11:43pm 5

Check with a smaller dt.

Topic		Replies	Views
Vern9 from differentialequations.jl and ode89 from matlab Modelling & Simulations	15	496	May 25, 2024
KenCarp47 is giving wrong result no matter how tight the tolerance is. Is there a way to resolve this? New to Julia	29	2223	January 25, 2022
Numerical instability of partial differential equations General Usage diffeq	6	1006	May 17, 2020
ODE system solvable with forward Euler but not general solvers (early stiffness) New to Julia diffeq	9	1365	April 2, 2022
Handling Instability When Solving ODE Problems Modelling & Simulations	14	32858	August 23, 2024

Second order differential equation blows up numerically, and I have no clue why. Any ideas?

Related topics