Integrating "Beyond Infinity" and Back

Consider the simple ordinary differential equation with initial condition . The solution is obviously given by . It tends to infinity as tends to 1. We now consider a time-discrete approximation to this solution as provided by the asynchronous leapfrog method. Since the evolution step formula does not involve potentially undefined operations, any such discrete approximation is well defined for all its values. Of course, the approximated values grow dramatically and will soon transcend what can be represented even with Mathematica's arbitrary-precision numbers. Since the asynchronous leapfrog method is a reversible integration method, we should be able to go along each finite discrete trajectory back to its initial point. If the final point was "close to infinity", the computation needs to be done with a large number of digits in order to come back to its initial point. This is what the present Demonstration studies. Values are input by means of setter bars instead of sliders since updating the curve for higher precision and may take a few seconds. (This does not affect Autorun, which uses machine precision throughout.) The reversed trajectory is marked by red dots and one easily sees (if the box 'show also the reversed trajectory' is activated) whether the reversed trajectory reaches the initial point. In all cases where it fails to do so, increasing the value of precision will finally solve the problem. In some cases, precision up to 2000 is needed.
As the integration method is set up here, it creates a trajectory that always grows, even if we have crossed the singularity of the exact solution. Referring to this situation as going beyond infinity, as I did in the title, is, of course, not to be taken literally. A slight modification of the asynchronous leapfrog method can allow the method to "integrate across the singularity" and to approximate the exact solution well also for . This will be analyzed in a future Demonstration.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


Snapshot 1: Maximum stress, reaching the largest numbers. The integration is from to with ; precision 1000 turns out to be insufficient for reversibility.
Snapshot 2: Reversible trajectory that gets lost in machine precision. Here the integration is from to with ; due to the larger value of , there are fewer steps in the critical domain and the values are less gigantic; so a precision of 250 is sufficient for getting a reversible trajectory.
Snapshot 3: Using machine precision, together with the remaining settings of Snapshot 2, spoils reversibility.
[1] Ulrich Mutze An Asynchronous Leap-Frog Method, 2008., http://www.ma.utexas.edu/mp_arc/c/08/08-197.pdf.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2017 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+