Avoiding Pitfalls In Periodic Orbits of Maps With A Piecewise Approach

Numerical computation is widely recognized to be of great importance in many fields of science. Many conclusions in nonlinear science and complex systems have been drawn upon simulation on a digital computer. Although this importance is widely recognized, the way of computers work is, in general, not known in detail.

Let us start our discussion with the following example. Let be the following equation, known as a logistic map:


The idea behind this equation is very simple. We start a simulation at stage k = 0. In this stage, we have our initial condition. Suppose that r = 4, and It is easy to calculate that in the next stage, we will get ¼. After the third iteration, the correct answer is going to be stuck — fixed — with the value of ¾.

Using a commercial software, we proceed with this calculation, and we can get the results in Figure 1. The red line and the black line are produced by the equations F(xk) = 4xk(1 – xk), while G(xk) = 4xk – 4x2k. Although the equations are mathematically equivalent, the computer produces a stable and fixed result for F and a chaotic behavior for G. Nonlinear dynamic theory and a branch of math based only on calculus and differential equations cannot explain this properly.

Figure 1: Simulation of the logistic map from two different computer realizations. Image courtesy Erivelton Geraldo Nepomuceno

In our recent study entitled “Interval computing periodic orbits of maps using a piecewise approach,” we recognize this situation and developed an approach to, at least, decrease this situation. The basic idea relies on the fact that the majority of numerical computations based on floating-point are built to guarantee a good result, or, in technical terms, a result within a precision according to the number of bits used, only for basic arithmetical operations.


The example in Figure 1 presents two features that turn this example in a very challenging one. First, the computation is a set of basic arithmetic operations. Second and most important, there is a recurrence where a very tiny error in the beginning grows exponentially. In our paper, we examine some features of the function, which allows us to set bounds of the result with more consistency than other approaches. The function is seen as a piecewise function. In general, the proposed method has produced intervals that are significantly narrower than those obtained by Intlab (a Matlab toolbox) approach. However, it is also clear that our method requires a significantly larger number of iterates.

Related Study

In an accompanying article published in the same journal “Detecting unreliable computer simulations of recursive functions with interval extensions,”, we used such ideas to show that some previous results that guarantee long simulations on a computer are not valid for current software and hardware. We also exploit the error propagation to compute the Lyapunov exponent, a usual feature to identify chaotic systems, in the article “Computation of the largest positive Lyapunov exponent using rounding mode and recursive least square algorithm”,

Future Directions

We are confident that computers are incredible tools to develop science and technology. But we are also confident that we must understand the way that they work in order to be able to analyze and avoid pitfalls.

These findings are described in the article entitled Interval computing periodic orbits of maps using a piecewise approach, recently published in the journal Applied Mathematics and ComputationThis work was conducted by Erivelton G. Nepomuceno, Heitor M. Rodrigues Junior, and Samir A.M. Martins from the  Federal University of São João del-ReiMatjaž Perc from the University of Maribor and Beihang University, and Mitja Slavinec from the University of Maribor.




Free-Living Amoebae – A Novel Waterborne Risk Pathway That Now Includes Enteric Viruses

When it comes to waterborne pathogens and host-specific pathogens like human enteric viruses, it is easy to assume that they […]

NASA Astronauts Just Found Life On The Surface Of The International Space Station

On Earth, life always seems to find a way to survive in all sorts of environments. From the coldest place […]

The End Of The Dinosaurs Marked The End Of Most Nocturnal Mammals

We wake up from a night’s sleep and go about our day until it becomes dark again and we head […]

Tropical Savanna: Animals Plants & Climate

A tropical savanna is an environment characterized by rolling hills, tall grasslands, and sparse trees. Large parts of Africa and […]

How An Echolocating Bat’s Brain May Solve The Microsecond Timing Problem

On summer evenings, before it gets completely dark, you may have spotted some flying creatures overhead. Although these creatures could […]

How Sound Spreads Is Reflected In Large Open Spaces

In a room, we typically hear that a voice is louder than in an open area. This is the effect […]

Communist Countries Around The World In 2019

Currently in the world today there are five countries classified as communist countries: China, Cuba, Laos, Korea, and Vietnam. But […]

Science Trends is a popular source of science news and education around the world. We cover everything from solar power cell technology to climate change to cancer research. We help hundreds of thousands of people every month learn about the world we live in and the latest scientific breakthroughs. Want to know more?