Following Chaotic Reflections

The path traveled by a photon bouncing off a grid of circular mirrors is very sensitive to the initial conditions. The paths shown are of a particle that starts at , travels east at speed 1, and reflects off mirrors of radius around integer lattice points in the plane. The paths correspond to initial conditions obtained by perturbing each of (the starting position) and (the starting direction) to five points around a circle of radius , where is controlled by the first slider. Thus there are 25 paths shown. When is not so small the paths diverge, which indicates that a interval is not small enough to guarantee accuracy. When is so small that all 25 paths coincide, it indicates that an interval of that size around the initial conditions is enough to guarantee accuracy in the final result. This illustrates how an interval-based algorithm can lead to correct trajectories. Indeed, the interval algorithm is used to compute all trajectories shown, so they are all correct.


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


The problem of determining 10 digits of the particle's position at time 10 was one of the SIAM 100-Digit Challenge problems. Using machine precision is inadequate to solve this. Interval arithmetic can be used to provide an algorithm that leads to proved-correct results, even if 10 (the time or the number of digits) is replaced by 1000 or more. The algorithm was developed by F. Bornemann and S. Wagon; see The SIAM 100-Digit Challenge: A Study in High-Accuracy Numerical Computing, by D. Laurie, F. Bornemann, S. Wagon, and J. Waldvogel, SIAM, Philadelphia, 2004.
The three snapshots show trajectories out to time 10, with initial spread of , where p=4, 11, and 14. In the last case the small spread indicates convergence to visual accuracy.
    • 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 © 2018 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+