<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AM</journal-id><journal-title-group><journal-title>Applied Mathematics</journal-title></journal-title-group><issn pub-type="epub">2152-7385</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/am.2023.149038</article-id><article-id pub-id-type="publisher-id">AM-128072</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Analysis of Dynamical Behavior of One-Dimensional Real Maps: An Executable Dynamical Programming Software Approach
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Mohammad</surname><given-names>Sharif Ullah</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Masuda</surname><given-names>Akter</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>K.</surname><given-names>M. Ariful Kabir</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Department of Mathematics, Bangladesh University of Engineering and Technology, Dhaka, Bangladesh</addr-line></aff><aff id="aff1"><addr-line>Department of Mathematics, Feni University, Feni, Bangladesh</addr-line></aff><pub-date pub-type="epub"><day>14</day><month>09</month><year>2023</year></pub-date><volume>14</volume><issue>09</issue><fpage>652</fpage><lpage>672</lpage><history><date date-type="received"><day>23,</day>	<month>August</month>	<year>2023</year></date><date date-type="rev-recd"><day>25,</day>	<month>September</month>	<year>2023</year>	</date><date date-type="accepted"><day>28,</day>	<month>September</month>	<year>2023</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The dynamical behavior of real-world phenomena is implausible graphically due to the complexity of mathematical coding. The present article has mainly focused on some one-dimensional real maps’ dynamical behavior irrespective of using coding. In continuation, linear, quadratic, cubic, higher-order, exponential, logarithmic, and absolute value maps have been used to scrutinize their dynamical behavior, including the characteristics of the orbit of points. Dynamical programming software (
  DPS.exe) will be proposed as a new technique to ascertain the dynamical behavior of said maps. Thus, a mathematician can automatically determine one-dimensional real maps’ dynamical behavior apart from complicated programming code and analytical solutions.
 
</p></abstract><kwd-group><kwd>One-Dimensional Map</kwd><kwd> Cobweb</kwd><kwd> Orbit Diagram</kwd><kwd> Fixed Point</kwd><kwd> the Fate of the Orbit</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Mathematical equations analytically reveal the idea of numerous expectations for modeling all-natural phenomena. In this regard, dynamical systems demonstrate a significant role and have an extensive and substantial aspect expressed by prominent mathematicians [<xref ref-type="bibr" rid="scirp.128072-ref1">1</xref>] - [<xref ref-type="bibr" rid="scirp.128072-ref12">12</xref>] .</p><p>One-dimensional maps play a crucial role [<xref ref-type="bibr" rid="scirp.128072-ref13">13</xref>] [<xref ref-type="bibr" rid="scirp.128072-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.128072-ref15">15</xref>] in predicting the natural behavior and physical object’s fate more precisely. The readers advised reading [<xref ref-type="bibr" rid="scirp.128072-ref16">16</xref>] [<xref ref-type="bibr" rid="scirp.128072-ref17">17</xref>] to know the periodic and aperiodic actions in discrete one-dimensional dynamical systems and the history of one-dimensional dynamics. Many studies [<xref ref-type="bibr" rid="scirp.128072-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.128072-ref19">19</xref>] [<xref ref-type="bibr" rid="scirp.128072-ref20">20</xref>] address real maps and dynamical behaviors with different approaches. Clark et al. prove the complex box bounds for real maps [<xref ref-type="bibr" rid="scirp.128072-ref21">21</xref>] . Iwanaga and Namatame signified evacuation decision-making contagion on a real map [<xref ref-type="bibr" rid="scirp.128072-ref22">22</xref>] . Jia et al. proposed a mobility model based on a real map for VANETs to overwhelm the existing model’s disadvantages [<xref ref-type="bibr" rid="scirp.128072-ref23">23</xref>] . Joshi and Blackmore effectively modeled the discrete evolution of space, biological, and ecological sciences by exponentially decaying discrete dynamical systems [<xref ref-type="bibr" rid="scirp.128072-ref24">24</xref>] .</p><p>Furthermore, many studies [<xref ref-type="bibr" rid="scirp.128072-ref25">25</xref>] [<xref ref-type="bibr" rid="scirp.128072-ref26">26</xref>] [<xref ref-type="bibr" rid="scirp.128072-ref27">27</xref>] [<xref ref-type="bibr" rid="scirp.128072-ref28">28</xref>] investigated one-dimensional map characteristics under different conditions. Sushko et al. discussed some basic concepts and definitions of non-smooth one-dimensional maps [<xref ref-type="bibr" rid="scirp.128072-ref29">29</xref>] . Some studies [<xref ref-type="bibr" rid="scirp.128072-ref30">30</xref>] - [<xref ref-type="bibr" rid="scirp.128072-ref35">35</xref>] introduced new techniques to discover dynamical map features. Medrano and Solis extended and improved the existing characterization of general quadratic actual polynomial maps dynamics with coefficients [<xref ref-type="bibr" rid="scirp.128072-ref36">36</xref>] . Bai et al. [<xref ref-type="bibr" rid="scirp.128072-ref37">37</xref>] analyze the invariant solutions of Coupled Burgers’ equations utilizing one-dimensional optimum systems. The ground-state energy and entropy for a one-dimensional Heisenberg chain with alternating D-terms are investigated by Xiang et al. [<xref ref-type="bibr" rid="scirp.128072-ref38">38</xref>] .</p><p>Moreover, analyzing dynamical behaviors, such as fixed-point, iteration, orbit under specific values, and the orbit’s fate, is challenging due to the complicated mathematical calculation and programming codes [<xref ref-type="bibr" rid="scirp.128072-ref39">39</xref>] [<xref ref-type="bibr" rid="scirp.128072-ref40">40</xref>] . Therefore, in the present study, one-dimensional real map-based techniques are proposed to determine their dynamical behavior without complicated programming, compressing a mathematician or physicist’s effort.</p><p>The progression of the current research work is as follows. The formulation is thoroughly described in Section 2. Section 3 offers a numerical and graphical discussion of the maps mentioned earlier. On top of that, we provide a detailed comparison between numerical, visual, and DPS.exe analysis. The final words are given in Section 4.</p></sec><sec id="s2"><title>2. Methodology</title><p>In any research, one of the best unspoken tools is arriving at reliant elucidations to the problems through systematic assortment and analysis. Firstly, the dynamical behavior of one-dimensional stated maps is discussed using different coding software [<xref ref-type="bibr" rid="scirp.128072-ref41">41</xref>] - [<xref ref-type="bibr" rid="scirp.128072-ref45">45</xref>] . Then, an executable FORTRAN coding system is used in the background of the newly proposed software. In this regard, the two algorithms are present. Finally, a comparison of graphical, numerical, and proposed software is illustrated for the said maps. The newly suggested MS-Dos software allows mathematicians to determine the above behavior of various one-dimensional real maps except for any complicated code.</p>Developing Dynamical Programming Software (DPS.exe)<p>The first requirement is to introduce the works, for example, a flowchart to classify functions’ essence to explore the dynamical simulation framework. The diagram (<xref ref-type="fig" rid="fig1"><xref ref-type="fig" rid="fig">Figure </xref>1</xref>) depicts the developing technique and application of the process of the newly proposed DPS.exe.</p><p>Case-I: One-dimensional first-degree map</p><p>The general form of the one-dimensional first-degree equation is y = f ( x ) = a x + b , where a and b are the real constants, and x is the variable.</p><p>Fixed point analysis</p><p>A specific value of a , b the FORTRAN command [<xref ref-type="bibr" rid="scirp.128072-ref44">44</xref>] gives the output <xref ref-type="fig" rid="fig2"><xref ref-type="fig" rid="fig">Figure </xref>2</xref>(i). But in this case, if b = 0 and x = 1, then f ( x ) = x and thus all the</p><p>points of f ( x ) will be the fixed points, which is x = b 1 − a ( a ≠ 1 ) . This condition will be ( a ≠ 1 ) overcome using the IF statement in the first line in FORTRAN command [<xref ref-type="bibr" rid="scirp.128072-ref44">44</xref>] . Again, when f ( x ) = y = x there is no fixed point, i.e., the parallel lines meet at infinity. If b 1 − a &gt; 10000000 or b 1 − a &lt; − 10000000 (this range can be changed for more reliable calculation), the line is eventually parallel, so no fixed point exists. Nevertheless, if − 1000000 &lt; b 1 − a &lt; 1000000 (under consideration), then there must be a fixed point presented by s c t = b 1 − a . The nature of the fixed point [<xref ref-type="bibr" rid="scirp.128072-ref16">16</xref>] (attracting, repelling, or neutral) will be determined using the following conditions:</p><p>x 0 = s c t = b 1 − a is { attracting   ;       if   | f ′ ( x 0 ) | &lt; 1 repelling   ;           if   | f ′ ( x 0 ) | &gt; 1 neutral   ;               if   | f ′ ( x 0 ) | = 1 .</p><p>The nature of the fixed point entirely depends on the value of a as f ′ ( x ) = a , f ′ ( x 0 ) = a . Therefore, the output of the FORTRAN executable (DPS.exe) interface <xref ref-type="fig" rid="fig2"><xref ref-type="fig" rid="fig">Figure </xref>2</xref>(ii).</p><p>Orbit analysis</p><p>Let x = x 0 be the initial seed. Then the orbit analysis of f ( x ) = a x + b is x 0 → f ( x 0 ) = a x 0 + b , f 2 ( x 0 ) = f ( f ( x 0 ) ) = a ( a x 0 + b ) + b = a 2 x 0 + a b + b and so on.</p><p>The output of this segment for f ( x ) = 2 x − 3 with the initial seed x 0 = 3.5 is presented in <xref ref-type="fig" rid="fig2"><xref ref-type="fig" rid="fig">Figure </xref>2</xref>(iii).</p><p>Fate of orbit</p><p>After continuing the iteration process sufficiently many more times, finally, the fate of the orbit is presented in <xref ref-type="fig" rid="fig2"><xref ref-type="fig" rid="fig">Figure </xref>2</xref>(iv). Now, users may need to repeat the process for any new function. This programming procedure automatically returns to the initial stage <xref ref-type="fig" rid="fig2"><xref ref-type="fig" rid="fig">Figure </xref>2</xref>(v). This section’s output proceeds the mathematician to the end of a program or the program’s initial phase.</p><p>Case-II: One-dimensional second-degree map</p><p>The general form of the one-dimensional second-degree equation is y = f ( x ) = a x 2 + b x + c , where a, b and c are the real constants and x is the variable.</p><p>Fixed point analysis</p><p>The specific values of a , b , c FORTRAN [<xref ref-type="bibr" rid="scirp.128072-ref46">46</xref>] give the output <xref ref-type="fig" rid="fig3"><xref ref-type="fig" rid="fig">Figure </xref>3</xref>(i).</p><p>Fixed point of f ( x ) = a x 2 + b x + c is x = − ( b − 1 ) &#177; ( b − 1 ) 2 − 4 a c 2 a .</p><p>When ( b − 1 ) 2 − 4 a c &gt; 0 then two fixed points exist, and those two fixed points are</p><p>− ( b − 1 ) + ( b − 1 ) 2 − 4 a c 2 a and − ( b − 1 ) − ( b − 1 ) 2 − 4 a c 2 a</p><p>Therefore, the output of the FORTRAN executable (DPS.exe) file is pictured in <xref ref-type="fig" rid="fig3"><xref ref-type="fig" rid="fig">Figure </xref>3</xref>(ii).</p><p>Orbit analysis</p><p>The output of this segment f ( x ) with the initial seed x 0 = 1.7 demonstrated in <xref ref-type="fig" rid="fig3"><xref ref-type="fig" rid="fig">Figure </xref>3</xref>(iii).</p><p>Case-III: One-dimensional third-degree maps</p><p>The general form of the one-dimensional third-degree equation is y = f ( x ) = a x 3 + b x 2 + c x + d where a, b, c, and d are the real constants and x is the variable.</p><p>Fixed point analysis</p><p>For a specific value of a , b , c , d FORTRAN command [<xref ref-type="bibr" rid="scirp.128072-ref46">46</xref>] generates the <xref ref-type="fig" rid="fig4"><xref ref-type="fig" rid="fig">Figure </xref>4</xref>(i).</p><p>Root process for finding fixed points</p><p>The fixed point is the point of intersection of y = f ( x ) and y = x . Using Mathematica or any other programming command [<xref ref-type="bibr" rid="scirp.128072-ref41">41</xref>] - [<xref ref-type="bibr" rid="scirp.128072-ref45">45</xref>] , one can find its fixed points. As it is complicated and lengthy, the numerical procedure may help to obtain the solution.</p><p>Numerical process for finding fixed points</p><p>The solution of finding the given equation’s solution is to set an initial value of x. This value maybe −10,000 or less. Now choose f ( x ) = a x 3 + b x 2 + c x + d , g ( x ) = x .</p><p>If f ( x ) = g ( x ) , then x is a fixed point, start checking with −10,000. If both f ( x ) and g ( x ) are not equal, then do the process for x = − 10000 + 0.00001</p><p>Similarly, if it is not equal yet, then do it again for x = − 10000.00001 + 0.00001</p><p>All these procedures can be quickly done using the FORTRAN command [<xref ref-type="bibr" rid="scirp.128072-ref46">46</xref>] . If any fixed point can be found, then the nature of the fixed point can be determined by the logic of | f ′ ( x 0 ) | .</p><p>For the specific function f ( x ) = x 3 + 2 x 2 − 3 x + 4 , the output is revealed in</p><p><xref ref-type="fig" rid="fig4"><xref ref-type="fig" rid="fig">Figure </xref>4</xref>(ii).</p><p>Orbit analysis and the fate of the orbit</p><p>If the iterative value for any specific function goes to &lt;−10<sup>12</sup> or goes to &gt;10<sup>12</sup> after some iterations, then the fate of the orbit goes to negative infinity or positive infinity, respectively. For the specific function f ( x ) with the initial seed x 0 = − 3.3 , the output is visualized in <xref ref-type="fig" rid="fig4"><xref ref-type="fig" rid="fig">Figure </xref>4</xref>(iii).</p><p>Case-IV: One-dimensional higher degree maps</p><p>One-dimensional higher degree equation can be expressed in the following form:</p><p>f ( x ) = a 0 + a 1 x + a 2 x 2 + a 3 x 3 + ⋯ + a n x n</p><p>where a 0 , a 1 , a 2 , ⋯ , a n are n numbers of coefficients and x is the variable. Developing DPS.exe for the one-dimensional higher degree function is more complicated, as described in the later section.</p><p>Equation generating technique</p><p>A glance at the development of DPS.exe for one-dimensional higher-degree maps has been described in <xref ref-type="fig" rid="fig5"><xref ref-type="fig" rid="fig">Figure </xref>5</xref>.</p><p>Suppose anyone is interested to know the dynamical behavior of a function of the 5<sup>th</sup> degree. Then executing this part of programming displays <xref ref-type="fig" rid="fig6"><xref ref-type="fig" rid="fig">Figure </xref>6</xref>(i) and needs to input the value of n as 5. Here the number of the variable associated with each term depends on the desire of any individual. So there needs to build an array of variables A <xref ref-type="fig" rid="fig6"><xref ref-type="fig" rid="fig">Figure </xref>6</xref>(ii) such as A (10,000). Now the focus is on the value of those variables. Anyone needs to input the values of variables for any specific degree function. Inputting the values of the associated variables, visualize the complete process and associated programming code stored in DPS.exe engine code. Then for particular values n = 5 and coefficients, A 0 = 2 , A 1 = − 1 , A 2 = 3 , A 3 = − 4 , A 4 = 1 , A 5 = 7 , the entire function is demonstrated in <xref ref-type="fig" rid="fig6"><xref ref-type="fig" rid="fig">Figure </xref>6</xref>(iii).</p><p>Numerical procedure obtaining fixed point</p><p>This procedure is identical to the numerical process of the third-degree equation, the generalized form of f ′ ( x ) is f ′ ( x ) = ∑ n = 0 n − 1 n a n x n − 1 . For the specific 5<sup>th</sup>-degree function, the output of the following programming segment is portrayed in <xref ref-type="fig" rid="fig6"><xref ref-type="fig" rid="fig">Figure </xref>6</xref>(iv).</p><p>Orbit analysis and the fate of the orbit</p><p>This procedure is equivalent to the third-degree equation, and the output of this segment of engine code unfolds in <xref ref-type="fig" rid="fig6"><xref ref-type="fig" rid="fig">Figure </xref>6</xref>(v).</p><p>Case-V: Experiment on Higher degree function</p><p>Here, the 6<sup>th</sup>-degree equation f ( x ) = 3 + 2 x − 7 x 2 − 1.5 x 3 + 2.2 x 4 − 3.7 x 5 + 1.1 x 6 has been considered. Then DPS.exe exhibits the dynamical info in <xref ref-type="fig" rid="fig7"><xref ref-type="fig" rid="fig">Figure </xref>7</xref>(i) an <xref ref-type="fig" rid="fig7"><xref ref-type="fig" rid="fig">Figure </xref>7</xref>(ii). To see the orbit for any specific initial seed, press 5, input the number of iterations (12, but it depends on the user’s desire), and the initial seed’s value ( x 0 = 3.1 ). The desired interface is in <xref ref-type="fig" rid="fig7"><xref ref-type="fig" rid="fig">Figure </xref>7</xref>(iii).</p><p>Case-VI: Exponential maps</p><p>The generalized form of an exponential map is f ( x ) = a e b x + c , where a , b , c are the arbitrary constants. The following source code asks the user of DPS.exe for specific values of a , b , c and finally expresses the function. The process of</p><p>finding the fixed point, nature of the fixed point, orbits, and fate of the orbit of f ( x ) under a specific initial seed is the same as mentioned previously. The dynamical behavior of f ( x ) = 3 e 1.1 x − 2 appears in <xref ref-type="fig" rid="fig8"><xref ref-type="fig" rid="fig">Figure </xref>8</xref>(i). After pressing 5, the system will represent the orbit analysis for any particular seed x 0 = − 0.9 in <xref ref-type="fig" rid="fig8"><xref ref-type="fig" rid="fig">Figure </xref>8</xref>(ii). Analogously, anyone can determine any exponential functions dynamical behavior by changing the value of coefficients.</p><p>Case-VII: Logarithmic maps</p><p>The generalized form of the exponential map is f ( x ) = a log ( b x ) + c , where a , b , c are arbitrary constants. The dynamical behavior of f ( x ) = 2 log ( 3 x ) + 4 appeared in <xref ref-type="fig" rid="fig9"><xref ref-type="fig" rid="fig">Figure </xref>9</xref>(i) an <xref ref-type="fig" rid="fig9"><xref ref-type="fig" rid="fig">Figure </xref>9</xref>(ii). Here, the number of iteration is 15, and the initial seed is 5.</p><p>Case-VIII: Absolute value maps</p><p>The generalized form of the absolute value map is f ( x ) = | a x b | + c , where a , b , c are arbitrary constants. The dynamical behavior of f ( x ) = | − 2 x 3 | − 4 demonstrated in <xref ref-type="fig" rid="fig1"><xref ref-type="fig" rid="fig">Figure </xref>1</xref>0.</p></sec><sec id="s3"><title>3. Result and Discussions</title><p>Exploring the exactness of the obtained result using DPS.exe has to compare it numerically and graphically. In numerical cases, the initial seed’s specific value</p><p>gives the following values f ( x ) , f 2 ( x ) = f ( f ( x ) ) , f 3 ( x ) = f ( f ( f ( x ) ) ) , ⋯ etc. On the other hand, graphical analysis shows the graph of a function, fixed point, and orbit of the fixed point under a specific initial seed. It is more apparent to determine the dynamical behavior from its graphical analysis. However, DPS.exe has been more straightforward for a mathematician to gather all the information about dynamical behavior without programming knowledge.</p><sec id="s3_1"><title>3.1. One-Dimensional First-Degree Map</title><p>Suppose the one-dimensional first-degree equation is f ( x ) = 2 x + 1 and the initial seed x 0 .</p><p>Numerical analysis</p><p>Therefore, f ( x ) = 2 x + 1 , f 2 ( x ) = 4 x + 3 , f 3 ( x ) = 8 x + 7 , f 4 ( x ) = 16 x + 5 , ⋯ , and so on. f ( x ) = x gives the desired fixed point, and the fixed point is x = − 1 . This fixed point is repelling because if x 0 = − 1.1 that is a nearby point of the initial seed, then the orbit of the function appears as follows:</p><p>− 1.1 → − 1.2 → − 1.4 → − 1.8 → ⋯ and so on.</p><p>Thus, the orbit of the function under the considered initial seed is −∞.</p><p>Graphical Analysis</p><p>The graphical representation of f ( x ) = 2 x + 1 and its dynamical behavior ensues in <xref ref-type="fig" rid="fig1"><xref ref-type="fig" rid="fig">Figure </xref>1</xref>1.</p><p>Here, the orbit of the point for the given function is repelling, represented by the blue staircase.</p><p>DPS.exe analysis</p><p>In this process, mathematicians need not apply any programming command, just run DPS.exe and insert the coefficients. As f ( x ) is a linear function, after clicking DPS.exe, press 1 for the linear function section, which is exhibited in the appendix (A-I <xref ref-type="fig" rid="fig">Figure </xref>A1.1).</p><p>Now for the function f ( x ) insert “2” as the value of “a”, and “1” as “b”. The computer will then do the rest of the job to determine all dynamical behavior, displayed in the appendix (A-I <xref ref-type="fig" rid="fig">Figure </xref>A1.2).</p><p>DPS.exe also offers to see the orbit of the function for any desired initial seed. For this, the user has to press “5” and enter. Then, insert the initial seed and the number of iterations. Finally, the appendix demonstrates the interface (A-I <xref ref-type="fig" rid="fig">Figure </xref>A1.3).</p><p>Finally, all comparisons of one-dimensional first-degree maps, namely, numerical, graphical, and DPS.exe are presented in <xref ref-type="table" rid="table1">Table 1</xref>.</p></sec><sec id="s3_2"><title>3.2. One-Dimensional Second-Degree Map</title><p>Suppose the one-dimensional second-degree equation is f ( x ) = 1 2 x 2 − 1 . Then,</p><p>as previously, the DPS.exe interface is displayed in the appendix (A-II Figures A2.1-A2.3), and all comparisons of one-dimensional second-degree maps, namely, numerical, graphical, and DPS.exe are presented in <xref ref-type="table" rid="table2">Table 2</xref>.</p></sec><sec id="s3_3"><title>3.3. Higher Degree Maps</title><p>Suppose the one-dimensional higher (4<sup>th</sup>) degree map is f ( x ) = 2 x 4 − 3 x 3 − 4 x 2 − 5 x − 7 . The DPS.exe interface is manifested in the appendix (A-III Figures A3.1-A3.4), and all comparisons of one-dimensional higher-degree maps are presented in <xref ref-type="table" rid="table3">Table 3</xref>.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Comparison between numerical, graphical, and DPS.exe analysis</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >Numerical analysis</th><th align="center" valign="middle" >Graphical analysis</th><th align="center" valign="middle" >DPS.exe analysis</th></tr></thead><tr><td align="center" valign="middle" >Fixed point</td><td align="center" valign="middle" >−1.0 (have to solve f ( x ) = x for x).</td><td align="center" valign="middle" >−1.0 (additional programming required).</td><td align="center" valign="middle" >−1.0 (system generated automatically).</td></tr><tr><td align="center" valign="middle" >Nature of Fixed point</td><td align="center" valign="middle" >Repelling (determined from assumption).</td><td align="center" valign="middle" >Repelling (it is revealed by taking an initial seed near the fixed point and using a graphical cobweb).</td><td align="center" valign="middle" >Repelling (the system determines this and automatically represents its nature to the user).</td></tr><tr><td align="center" valign="middle" >Orbit for initial seed x<sub>0</sub> = −1.1</td><td align="center" valign="middle" >− 1.1 → − 1.4 → − 1.8 → ⋯ difficult to continue this process.</td><td align="center" valign="middle" >Values inserted in the list variable can show using the programming, but a little bit complicated.</td><td align="center" valign="middle" >The system automatically generates the orbit of f ( x ) for a given initial seed inserted by the user. Furthermore, the result is the same as the numerical process.</td></tr><tr><td align="center" valign="middle" >The fate of orbit for initial seed x<sub>0</sub> = −1.1</td><td align="center" valign="middle" >A sufficient number of iterations is essential.</td><td align="center" valign="middle" >As the cobweb is heading towards the significant negative, assume that fate is −∞.</td><td align="center" valign="middle" >DPS.exe automatically calculates itself and generates the result as −∞.</td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Comparison between numerical, graphical, and DPS.exe analysis</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >Numerical Analysis</th><th align="center" valign="middle" >Graphical Analysis</th><th align="center" valign="middle" >DPS.exe Analysis</th></tr></thead><tr><td align="center" valign="middle" >Fixed point</td><td align="center" valign="middle" >−0.73205 and 2.73205 (have to solve f ( x ) = x for x).</td><td align="center" valign="middle" >−0.73205 and 2.73205 (additional programming required).</td><td align="center" valign="middle" >−0.73205 and 2.73205 (system generated automatically).</td></tr><tr><td align="center" valign="middle" >Nature of the fixed point</td><td align="center" valign="middle" >A large scale of calculation and assumption is required.</td><td align="center" valign="middle" >−0.73205 is attracting fixed points, whereas 2.73205 is repelling, which arises by taking an initial seed near the fixed point and additional programming required for graphical analysis.</td><td align="center" valign="middle" >−0.73205 is attracting a fixed point, whereas 2.73205 is repelling. The system itself determines this and automatically represents its nature to the user.</td></tr><tr><td align="center" valign="middle" >Orbit for initial seed x<sub>0</sub> = −0.7</td><td align="center" valign="middle" >Complicate to continue this process.</td><td align="center" valign="middle" >Same as previous</td><td align="center" valign="middle" >Same as previous</td></tr><tr><td align="center" valign="middle" >The fate of orbit for initial see x<sub>0</sub> = −0.7</td><td align="center" valign="middle" >Same as previous</td><td align="center" valign="middle" >Assuming from cobweb</td><td align="center" valign="middle" >Same as previous</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Comparison between numerical, graphical, and DPS.exe analysis</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >Numerical analysis</th><th align="center" valign="middle" >Graphical analysis</th><th align="center" valign="middle" >DPS.exe analysis</th></tr></thead><tr><td align="center" valign="middle" >Fixed point</td><td align="center" valign="middle" >As previously, it is complicated as well as time-consuming. Furthermore, various methods have to apply to determine.</td><td align="center" valign="middle" >−1 and 2.77447 (additional programming required).</td><td align="center" valign="middle" >−1.000064 and 2.774464 (system generated automatically).</td></tr><tr><td align="center" valign="middle" >Nature of the fixed point</td><td align="center" valign="middle" >Same as previous</td><td align="center" valign="middle" >Both the fixed points are repelling, which arises by taking an initial seed near the fixed-point and additional programming required for graphical analysis.</td><td align="center" valign="middle" >Both the fixed points are repelling. The system itself determines this and automatically represents its nature to the user.</td></tr><tr><td align="center" valign="middle" >Orbit for initial seed x<sub>0</sub> = 0.5</td><td align="center" valign="middle" >Difficult to continue this process.</td><td align="center" valign="middle" >Same as previous</td><td align="center" valign="middle" >Same as previous</td></tr><tr><td align="center" valign="middle" >The fate of orbit for initial seed x<sub>0</sub> = 0.5</td><td align="center" valign="middle" >Same as previous</td><td align="center" valign="middle" >Same as previous</td><td align="center" valign="middle" >Same as previous</td></tr></tbody></table></table-wrap><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Comparison between numerical, graphical, and DPS.exe analysis</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >Numerical Analysis</th><th align="center" valign="middle" >Graphical Analysis</th><th align="center" valign="middle" >DPS.exe</th></tr></thead><tr><td align="center" valign="middle" >Fixed point</td><td align="center" valign="middle" >−1.99067 and 0.874516</td><td align="center" valign="middle" >Additional programming provides −1.99067 and 0.874516 are the fixed points.</td><td align="center" valign="middle" >−1.991655 and 0.874379 Generated by the system automatically.</td></tr><tr><td align="center" valign="middle" >Nature of Fixed point</td><td align="center" valign="middle" >One is attracting, and the other is repelling (from assumption).</td><td align="center" valign="middle" >One is attracting, and the other is repelling, which arises by taking an initial seed near the fixed-point and additional programming required for graphical analysis.</td><td align="center" valign="middle" >One is attracting, and the other is repelling. The system itself determines this and automatically represents its nature to the user.</td></tr><tr><td align="center" valign="middle" >Orbit for initial seed x<sub>0</sub> = −0.5</td><td align="center" valign="middle" >Same as previous</td><td align="center" valign="middle" >Same as previous</td><td align="center" valign="middle" >Same as previous</td></tr><tr><td align="center" valign="middle" >The fate of orbit for initial seed x<sub>0</sub> = −0.5</td><td align="center" valign="middle" >Same as previous</td><td align="center" valign="middle" >As the cobweb is heading towards the significant negative, the fate is −1.99067.</td><td align="center" valign="middle" >DPS.exe automatically calculates itself and generates the result as −1.99067.</td></tr></tbody></table></table-wrap></sec><sec id="s3_4"><title>3.4. Exponential Maps</title><p>Suppose the exponential function is f ( x ) = 1 2 e 2 x − 2 . Therefore, the DPS.exe interface is demonstrated in the appendix (A-IV Figures A4.1-A4.3), and all comparisons of the exponential map are presented in <xref ref-type="table" rid="table4">Table 4</xref>.</p></sec></sec><sec id="s4"><title>4. Conclusions</title><p>The one-dimensional real map is perceived as difference equations, iterated maps, or recursion relations in mathematical systems that model a single variable due to evolving over discrete steps. It has a remarkable significance in modeling natural phenomena, for example, population dynamics, electronics, and economics. However, this study has profoundly elaborated a possible one-dimensional real maps coding system to know the dynamical behavior and proposed a new technique, executable dynamical programming software in short DPS.exe. The appropriateness of the proposed DPS.exe is then systematically investigated graphically and numerically.</p><p>The present work conducted a theoretical, graphical, and extensive numerical analysis to comprehensively explore one-dimensional real maps of dynamical behavior: first-degree, second-degree, third-degree, nth-degree, exponential, logarithmic, and absolute. The main focus is on one-dimensional real maps to demonstrate dynamic behavior in the system. A sensible relationship between the graphical, numerical, and DPS.exe has drowned. Furthermore, DPS.exe is an effective software for determining the dynamical behavior of one-dimensional real maps rather than general calculating, Mathematica, or other programming languages. This analytical research suggests that the newly proposed MS-Dos software allows mathematicians and physicists to determine various one-dimensional real maps’ dynamical behavior without complicating programming code. We plan to analyze the chaotic maps using the current mechanism.</p></sec><sec id="s5"><title>Acknowledgements</title><p>The authors thank the anonymous reviewers for their suggestions and invaluable comments.</p></sec><sec id="s6"><title>Author Contributions</title><p>Each author equally contributed to this paper and read and approved the final manuscript.</p></sec><sec id="s7"><title>Declaration of Competing Interest</title><p>The authors declare that they have no known competing financial interests or personal relationships that could have influenced the work reported in this paper.</p></sec><sec id="s8"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s9"><title>Cite this paper</title><p>Ullah, M.S., Akter, M. and Ariful Kabir, K.M. (2023) Analysis of Dynamical Behavior of One-Dimensional Real Maps: An Executable Dynamical Programming Software Approach. Applied Mathematics, 14, 652-672. https://doi.org/10.4236/am.2023.149038</p></sec><sec id="s10"><title>Appendix A</title><p>A-I: One-dimensional first-degree map</p><p>A-II: One-dimensional second-degree map</p><p>A-III: Higher degree maps</p><p>A-IV: Exponential maps</p></sec></body><back><ref-list><title>References</title><ref id="scirp.128072-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Newton, I. (1729) The Mathematical Principles of Natural Philosophy (Vol. 2). 3rd Edition, William Dawson and Sons, London. https://www.google.com.hk/books/edition/The_Mathematical_Principles_of_Natural_P/Tm0FAAAAQAAJ?hl=en&amp;gbpv=0</mixed-citation></ref><ref id="scirp.128072-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Poincaré, H. (1967) New Methods of Celestial Mechanics. American Institute of Physics, Melville.</mixed-citation></ref><ref id="scirp.128072-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Poincaré, H. (1890) Chapitre I. Propriétésgénérales des équationsdifférentielles. Acta Mathematica, 13, 8-45. https://doi.org/10.1007/bf02392507</mixed-citation></ref><ref id="scirp.128072-ref4"><label>4</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Julia</surname><given-names> G. </given-names></name>,<etal>et al</etal>. (<year>1918</year>)<article-title>Mémoiresurl’iteration des fonctionsrationnelles</article-title><source> Journal de Mathématiques Pureset Appliquées</source><volume> 8</volume>,<fpage> 47</fpage>-<lpage>245</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.128072-ref5"><label>5</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Fatou</surname><given-names> P. </given-names></name>,<etal>et al</etal>. (<year>1917</year>)<article-title>Sur les substitutions rationnelles</article-title><source> Comptes Rendusdel Académie des Sciences de Paris</source><volume> 164</volume>,<fpage> 806</fpage>-<lpage>808</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.128072-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Birkhoff, G.D. (1920) Recent Advances in Dynamics. Science, 51, 51-55. https://doi.org/10.1126/science.51.1307.51</mixed-citation></ref><ref id="scirp.128072-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Smale, S. (1999) Smale Horseshoes and Symbolic Dynamics in Perturbed Nonlinear Schr&amp;#246;dinger Equations. Journal of Nonlinear Science, 9, 363-415. https://doi.org/10.1007/s003329900074</mixed-citation></ref><ref id="scirp.128072-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Lorenz, E.N. (1972) Predictability: Does the Flap of a Butterfly’s Wings in Brazil Set off a Tornado in Texas. American Association for the Advancement of Science, Washington DC.</mixed-citation></ref><ref id="scirp.128072-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Robert, M.M. (1976) Simple Mathematical Models with Very Complicated Dynamics. Nature, 261, 459-467. https://doi.org/10.1038/261459a0</mixed-citation></ref><ref id="scirp.128072-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Mandelbrot, B.B. (1983) The Fractal Geometry of Nature. American Journal of Physics, 51, 286-287. https://doi.org/10.1119/1.13295</mixed-citation></ref><ref id="scirp.128072-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Douady, A. and Hubbard, J.H. (1984) &amp;#201;tude dynamique des polyn&amp;#244;mes complexes. Université de Paris-Sud, Paris. https://pi.math.cornell.edu/~hubbard/OrsayFrench.pdf</mixed-citation></ref><ref id="scirp.128072-ref12"><label>12</label><mixed-citation publication-type="book" xlink:type="simple">Sullivan, J.M. (2012) Mathematical Pictures: Visualization, Art and Outreach. In: Behrends, E., Crato, N. and Rodrigues, J., Eds., Raising Public Awareness of Mathematics, Springer, Berlin, 279-293. https://doi.org/10.1007/978-3-642-25710-0_21</mixed-citation></ref><ref id="scirp.128072-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">de Melo, W. and van Strien, S. (1992) One-Dimensional Dynamics. Springer eBooks. https://doi.org/10.1007/978-3-642-78043-1</mixed-citation></ref><ref id="scirp.128072-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Sharkovsky, A.N., Kolyada, S.F., Sivak, A.G. and Fedorenko, V.V. (1997) Dynamics of One-Dimensional Maps. Springer, New York. https://doi.org/10.1007/978-94-015-8897-3</mixed-citation></ref><ref id="scirp.128072-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Rodrigues, A. (2021) One-Dimensional Dynamical Systems: An Example-Led Approach. Chapman and Hall/CRC Press, New York. https://doi.org/10.1201/9781003144618</mixed-citation></ref><ref id="scirp.128072-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Grandmont, J.M. (2021) Cycles and Chaos in Economic Equilibrium: Periodic and Aperiodic Behaviour in Discrete One-Dimensional Dynamical Systems. Princeton University Press eBooks, 44-63. https://doi.org/10.2307/j.ctv19fvxt1.5</mixed-citation></ref><ref id="scirp.128072-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Sharkovsky, A.N. (2014) On the History of One-Dimensional Dynamics. ESAIM: Proceedings and Surveys, 46, 83-85. https://doi.org/10.1051/proc/201446007</mixed-citation></ref><ref id="scirp.128072-ref18"><label>18</label><mixed-citation publication-type="book" xlink:type="simple">Elaydi, S. (2019) Global Dynamics of Discrete Dynamical Systems and Difference Equations. In: Elaydi, S., P&amp;#246;tzsche, C. and Sasu, A., Eds., ICDEA 2017: Difference Equations, Discrete Dynamical Systems and Applications, Springer, Cham, 51-81. https://doi.org/10.1007/978-3-030-20016-9_3</mixed-citation></ref><ref id="scirp.128072-ref19"><label>19</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Tanaka</surname><given-names> S. </given-names></name>,<etal>et al</etal>. (<year>2019</year>)<article-title>Brain as a Dynamical System</article-title><source> Journal of Brain and Nerves</source><volume> 71</volume>,<fpage> 657</fpage>-<lpage>664</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.128072-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Mehran, N.Z., Panahi, S., Hosseini, Z., Golpayegani, S.M.R.H. and Jafari, S. (2020) One Dimensional Map-Based Neuron Model: A Phase Space Interpretation. Chaos, Solitons &amp; Fractals, 132, Article ID: 109558. https://doi.org/10.1016/j.chaos.2019.109558</mixed-citation></ref><ref id="scirp.128072-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Clark, T., Strien, S.V. and Trejo, S. (2013) Complex Box Bounds for Real Maps. Communications in Mathematical Physics, 355, 1001-1119. https://doi.org/10.1007/s00220-017-2958-y</mixed-citation></ref><ref id="scirp.128072-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Iwanaga, S. and Namatame, A. (2016) Contagion of Evacuation Decision Making on Real Map. Mobile Networks and Applications, 21, 206-214. https://doi.org/10.1007/s11036-016-0704-x</mixed-citation></ref><ref id="scirp.128072-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Jia, X., Yujun, K., Enzhan, Z. and Weili, J. (2011) Research on Mobility Model for VANETs Based on Real Map. International Conference on Computational Problem-Solving (ICCP), Chengdu, 21-23 October 2011, 1-6. https://doi.org/10.1109/ICCPS.2011.6089939</mixed-citation></ref><ref id="scirp.128072-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Joshi, Y. and Blackmore, D. (2012) Exponentially Decaying Discrete Dynamical Systems. Recent Patents on Space Technology, 2, 37-48. https://doi.org/10.2174/1877611611202010037</mixed-citation></ref><ref id="scirp.128072-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Kozlovski, O. (2013) Periodic Attractors of Perturbed One-Dimensional Maps. Ergodic Theory and Dynamical Systems, 33, 1519-1541. https://doi.org/10.1017/etds.2013.28</mixed-citation></ref><ref id="scirp.128072-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Rivera-Letelier, J.R. and Shen, W. (2014) Statistical Properties of One-Dimensional Maps under Weak Hyperbolicity Assumptions. arXiv: 1004.0230.</mixed-citation></ref><ref id="scirp.128072-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Bruin, H. and Vejnar, B. (2020) Classification of One Dimensional Dynamical Systems by Countable Structures. arXiv: 2006.14926.</mixed-citation></ref><ref id="scirp.128072-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Ohmori, S. and Yamazaki, Y. (2020) Ultradiscrete Bifurcations for One Dimensional Dynamical Systems. Journal of Mathematical Physics, 61, Article ID: 122702. https://doi.org/10.1063/5.0012772</mixed-citation></ref><ref id="scirp.128072-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Sushko, I., Gardini, L. and Avrutin, V. (2016) Non-Smooth One-Dimensional Maps: Some Basic Concepts and Definitions. Journal of Difference Equations and Applications, 22, 1816-1870. https://doi.org/10.1080/10236198.2016.1248426</mixed-citation></ref><ref id="scirp.128072-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Li, Q., Tang, S. and Feng, X. (2012) Computing 1D Discontinuous Boundaries of Dynamical Systems. 2012 24th Chinese Control and Decision Conference (CCDC), Taiyuan, 23-25 May 2012, 1427-1430. https://doi.org/10.1109/CCDC.2012.6244229</mixed-citation></ref><ref id="scirp.128072-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Ufuktepe, &amp;#220;. (2014) Applications of Discrete Dynamical Systems with Mathematic&amp;#228;. Conference: RIMSAt, Kyoto, February 2014.</mixed-citation></ref><ref id="scirp.128072-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Jorba, &amp;#192;., Rabassa, P. and Tatjer, J.C. (2016) Local Study of a Renormalization Operator for 1D Maps under Quasiperiodic Forcing. Discrete &amp; Continuous Dynamical Systems, 9, 1171-1188. https://doi.org/10.3934/dcdss.2016047</mixed-citation></ref><ref id="scirp.128072-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">Khmou, Y., Said, S. and Frikel, M. (2018) A Comparison of Entropy Metrics in 1D Discrete Dynamical System. 1st International Conference on Signals, Automation and Telecommunications, Beni Mellal, 2-4 May 2018.</mixed-citation></ref><ref id="scirp.128072-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">Bashkirtseva, I. and Tsvetkov, I. (2018) Impact of the Parametric Noise on Map-Based Dynamical Systems. AIP Conference Proceedings, 2025, Article ID: 040004. https://doi.org/10.1063/1.5064888</mixed-citation></ref><ref id="scirp.128072-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">Ballard, T., Palada, H., Griffin, M. and Neal, A. (2019) An Integrated Approach to Testing Dynamic, Multilevel Theory: Using Computational Models to Connect Theory, Model, and Data. Organizational Research Methods, 24, 251-284. https://doi.org/10.1177/1094428119881209</mixed-citation></ref><ref id="scirp.128072-ref36"><label>36</label><mixed-citation publication-type="other" xlink:type="simple">Medrano, F.F. and Solis, F.J. (2015) Stability of Real Parametric Polynomial Discrete Dynamical Systems. Discrete Dynamics in Nature and Society, 2015, Article ID: 680970. https://doi.org/10.1155/2015/680970</mixed-citation></ref><ref id="scirp.128072-ref37"><label>37</label><mixed-citation publication-type="other" xlink:type="simple">Bai, Y., Bilige, S. and Chaolu, T. (2018) Potential Symmetries, One-Dimensional Optimal System and Invariant Solutions of the Coupled Burgers’ Equations. Journal of Applied Mathematics and Physics, 6, 1825-1839. https://doi.org/10.4236/jamp.2018.69156</mixed-citation></ref><ref id="scirp.128072-ref38"><label>38</label><mixed-citation publication-type="other" xlink:type="simple">Xiang, C. and Wang, H. (2019) Ground-State Energy and Entropy for One-Dimensional Heisenberg Chain with Alternating D-Term. Journal of Applied Mathematics and Physics, 7, 1220-1225. https://doi.org/10.4236/jamp.2019.75082</mixed-citation></ref><ref id="scirp.128072-ref39"><label>39</label><mixed-citation publication-type="other" xlink:type="simple">Devaney, R.L. and Choate, J. (2000) Chaos: A Tool Kit of Dynamics Activities. Key Curriculum Press, Emeryville.</mixed-citation></ref><ref id="scirp.128072-ref40"><label>40</label><mixed-citation publication-type="other" xlink:type="simple">Layek, G.C. (2015) An Introduction to Dynamical Systems and Chaos. Springer, New Delhi. https://doi.org/10.1007/978-81-322-2556-0</mixed-citation></ref><ref id="scirp.128072-ref41"><label>41</label><mixed-citation publication-type="other" xlink:type="simple">Lynch, S. (2007) Dynamical Systems with Applications Using Mathematica. Birkhauser, New York.</mixed-citation></ref><ref id="scirp.128072-ref42"><label>42</label><mixed-citation publication-type="other" xlink:type="simple">Kulenovic, M.R.S. and Merino, O. (2002) Discrete Dynamical Systems and Difference Equations with Mathematica. CRC Press, New York.</mixed-citation></ref><ref id="scirp.128072-ref43"><label>43</label><mixed-citation publication-type="other" xlink:type="simple">Lynch, S. (2004) Dynamical Systems with Applications Using MATLAB. Birkh&amp;#228;user, Boston. https://doi.org/10.1007/978-0-8176-8156-2</mixed-citation></ref><ref id="scirp.128072-ref44"><label>44</label><mixed-citation publication-type="other" xlink:type="simple">Lynch, S. (2010) Dynamical Systems with Applications Using Maple. Birkh&amp;#228;user, Boston. https://doi.org/10.1007/978-0-8176-4605-9</mixed-citation></ref><ref id="scirp.128072-ref45"><label>45</label><mixed-citation publication-type="other" xlink:type="simple">Nicola, B., Luigi, P. and Antonio, R. (2000) Mechanics and Dynamical Systems with Mathematica. Springer, New York.</mixed-citation></ref><ref id="scirp.128072-ref46"><label>46</label><mixed-citation publication-type="other" xlink:type="simple">Lipschutz, S. and Poe, A. (1982) Programming with FORTRAN. McGraw-Hill, Singapore.</mixed-citation></ref></ref-list></back></article>