# Questions tagged [ode]

Ordinary Differential Equations (ODEs) contain functions of only one independent variable, and one or more of their derivatives with respect to that variable. This tag is intended for questions on modeling phenomena with ODEs, solving ODEs, and other related aspects.

336 questions
Filter by
Sorted by
Tagged with
3k views

### What's the state of the art in parallel ODE methods?

I'm currently looking into parallel methods for ODE integration. There is a lot of new and old literature out there describing a wide range of approaches, but I haven't found any recent surveys or ...
4k views

### How can the gravitational n-body problem be solved in parallel?

How can the gravitational n-body problem be solved numerically in parallel? Is precision-complexity tradeoff possible? How does precision influence the quality of the model?
2k views

### How does one test a numerical ODE solver implementation?

I'm about to start working on a software library of numerical ODE solvers, and I'm struggling with how to formulate tests for the solver implementations. My ambition is that the library, eventually, ...
2k views

### What does “symplectic” mean in reference to numerical integrators, and does SciPy's odeint use them?

In this comment I wrote: ...default SciPy integrator, which I'm assuming only uses symplectic methods. in which I am refering to SciPy's odeint, which uses ...
6k views

### Why are higher-order Runge–Kutta methods not used more often?

I was just curious as to why high-order (i.e. greater than 4) Runge–Kutta methods are almost never discussed/employed (at least to my knowledge). I understand it requires greater computational time ...
4k views

### The definition of stiff ODE system

Consider an IVP for ODE system $y'=f(x,y)$, $y(x_0)=y_0$. Most commonly this problem is considered stiff when Jacobi matrix $\frac{\partial f}{\partial y}(x_0,y_0)$ has both eigenvalues with very ...
1k views

### Why does the numerical solution of an ODE move away from an unstable equilibrium?

I wish to simulate the behaviour of a double-pendulum-like system. The system is a 2-degrees-of-freedom robot manipulator that is not actuated and will, therefore, behave mostly like a double-pendulum ...
1k views

### Numerical methods for discontinuous r.s. ODEs

what are state of art methods for numerical solution of ODEs with discontinuous right side? I'm mostly interested piecewise-smooth right side functions, e.g. sign. I'm trying to solve the equation of ...
1k views

### Optimal ODE method for fixed number of RHS evaluations

In practice, the runtime of numerically solving an IVP $$\dot{x}(t) = f(t, x(t)) \quad \text{ for } t \in [t_0, t_1]$$ $$x(t_0) = x_0$$ is often dominated by the duration of evaluating the right-...
227 views

### Best practices for describing agent-based models

I work fairly heavily in mathematical biology/epidemiology, where most of the modeling/computational science work is still dominated by sets of ODEs, admittedly sometimes fairly elaborate sets of them....
213 views

### Algorithms for linear system of ODEs

I wonder: what is the best algorithm to solve $$\frac{du}{dt} = Au$$ Where $A$ is a real $n\times n$ matrix. A is not explicitly time dependent, usually sparse but not ...
2k views

### Is there an open source set of ODE solvers for C that use the native C99 complex type?

I've been using GSL as the foundation of many of my simulations, but it's a little bit overkill for my purposes and it defines its own complex type for legacy reasons. Rather than code my own Runge-...
3k views

### Choice of step size using ODEs in matlab

Hey there and thanks for giving time to look at my question. This is a updated version of my question which I posted earlier in physics.stackexchange.com I'm currently studying a 2D exciton spinor ...
6k views

### solving coupled ODEs with initial-value and final-value constraints

The essence of my question is the following: I have a system of two ODEs. One has an initial-value constraint and the other has a final-value constraint. This can be thought of as a single system with ...
480 views

### Runge-Kutta and Reusing Datapoints

I am trying to implement the fourth order Runge-Kutta method for solving a first order ODE in Python i.e. $\frac{dy}{dx} = f(x,y)$. I understand how the method works, but am trying to write an ...
1k views

### Numerically stable explicit solution of small linear system

I have an inhomogeneous linear system $$Ax=b$$ where $A$ is a real $n\times n$ matrix with $n\leq 4$. The nullspace of $A$ is guaranteed to be of zero dimension so the equation has a unique ...
842 views

### How do you improve the accuracy of a finite difference method for finding the eigensystem of a singular linear ODE

I am attempting to solve an equation of the type: $\left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x)$ Where $f(x)$ has a simple pole at $0$, for the ...
835 views

### Can I use an explicit time stepping scheme to determine numerically whether an ODE is stiff?

I have an ODE: $u'=-1000u+sin(t)$ $u(0)=-\frac{1}{1000001}$ I know that this particular ODE is stiff, analytically. I also know that if we use an explicit (forward) time stepping method (...
136 views

4k views

### ODEs vs DAE vs ADE?

I am totally confused between ODEs which I am familiar with, and differential algebraic equations (DAE) and Algebraic Differential Equations (ADE). Are they the same but just different names or what ...
2k views

### Is the shooting method the only general numerical method for solving nonlinear boundary value ODEs?

During my wandering in Mathematica.se, I gradually noticed that a certain kind of differential equation solving problem is "troubling" us all the time, that is, the boundary value problem (BVP) of ...
127 views

### How can I numerically solve an ODE to $N$ provably correct digits?

Suppose we have an initial value problem of the form $$\frac{\mathrm{d} \mathbf{x}}{\mathrm{d} t} = f(\mathbf{x}) \qquad \mathbf{x}(0) = \mathbf{x}_0$$ where $\mathbf{x}_0 \in \mathbb{R}^n$ is known ...
5k views

427 views

### ODE: How to measure stiffness if the Jacobian has zero eigenvalues?

Say you have a system of ODE's where the Jacobian has one zero eigenvalue; what does that tell you about the stiffness of the system? This case doesn't seem to be discussed in the cases I have been ...
328 views

### Numerical investigation of stability of motion (confinement)

I am trying to find the required specifications of a RF trap, in which a proton can be confined.(trap dimensions, voltage frequency and amplitude used, etc). I have to solve the equations of motion ...