INVARIANT CURVES AND FOCAL POINTS IN A LYNESS ITERATIVE PROCESS
|
|
- Hollie Jordan
- 5 years ago
- Views:
Transcription
1 International Journal of Bifurcation and Chaos, Vol. 13, No. 7 (2003) c World Scientific Publishing Company INVARIANT CURVES AND FOCAL POINTS IN A LYNESS ITERATIVE PROCESS L. GARDINI and G. I. BISCHI University of Urbino, Italy gardini@econ.uniurb.it bischi@econ.uniurb.it C. MIRA 19 rue d Occitanie, Fonsegrives, Quint, Toulouse Cedex, France christian.mira@insa-tlse.fr Received January 30, 2001; Revised October 22, 2001 We investigate the properties of recurrence of the type x n+1 = (a+ k 2 i=0 x n i)/x n (k 1), known as Lyness iterations from [Lyness, 1942, 1945, 1961] and recently analyzed by several authors in the case a > 0, see e.g. [Kocic et al., 1993; Csornyei & Laczkovich, 2000]. We reconsider Lyness recurrences at the light of some recent results on iterated maps with denominator, given in [Bischi et al., 1999a], where new kinds of singularities, such as focal points and prefocal curves, have been defined. In this paper, in particular, we give an answer to one of the open problems proposed in [Kocic & Ladas, 1993, pp. 141] concerning the dynamic behavior of Lyness recurrences for a < 0. We also give some new results in the case a > 0, and we improve a previous result on Lyness periodic recurrences. Keywords: Recurrences; rational maps; Lyness equations; focal points. 1. Introduction The dynamic behavior of some recurrence relations of order k (k 2) of the type x n+1 = k 2 α + b i x n i i=0 x n (k 1) (1) has been recently investigated by several authors (see e.g. [Kocic et al., 1993; Li & Liu, 1999; Csornyei & Laczkovich, 2001]). After the early works by Lyness [1942, 1945, 1961] several results have been established for recurrences of the type (1), mainly in order to detect bounded oscillatory behaviors of positive trajectories for non-negative values of the parameters. A survey of these results can be found in the book by Kocic and Ladas [1993, Chap. 5], where several interesting references are also given. Csornyei and Laczkovich [2001] give some results concerning special parameters values for which (1) is a periodic sequence, that is, all the infinite sequences which are solutions of (1) are periodic. The interest in the recurrences (1) is not only theoretical, but is often related to the fact that they are met in several applications, for example in ecologic modeling, see e.g. [Kocic et al., 1993], and references therein. The authors quoted above focused their attention on the existence of a positive equilibrium with oscillatory behavior around it, whereas the question whether positive and bounded trajectories exist when some parameter is negative was left as an open problem. The main purpose of the present paper is to give an answer to this and other questions concerning recurrences (1). We shall make 1841
2 1842 L. Gardini et al. use of some recent results on iterated maps with a denominator which can vanish, given in [Bischi et al., 1999a, 2000], where new kinds of singularities, called focal points and prefocal curves, have been defined (see also the related works by Bischi and Gardini [1997] and Mira [1999]). As shown throughout the paper, some interesting features of the recurrences (1) can be evidenced at the light of these new concepts, which allow us to stress some global structures of the phase curves on which the sequences generated by (1) are confined. For sake of simplicity we shall mainly focus on the simplest case, obtained for k = 2 in (1). However, many of the techniques and properties that we present in this paper can be generalized to recurrences (1) with higher values of k, i.e. recurrences of order greater than 2. Some possible generalizations are explicitly shown in the paper. The plane of the work is as follows. In Sec. 2 we present the simplest case, we recall some known results given in the literature, and we state the new results concerning the case a < 0. Such new results are proved in Sec. 3, where some properties and the terminology introduced in [Bischi et al., 1999a] are briefly described and extensively used. In Sec. 4 we show that the results and methods contained in Secs. 2 and 3 can also be applied to the case of positive parameter, thus confirming and extending some results given in the literature. In Sec. 5 some extensions are given for higher order recurrences (1) with k 2. In Sec. 6 we shall consider the particular case with k = 3, we prove that an uncountable many periodic solutions exist for any value of the parameter a, and we improve a result concerning periodic recurrences given by Lyness (see [Kocic et al., 1993] or [Csornyei & Laczkovich, 2001]). 2. The Simplest Recurrence Let us consider the recurrence (1) with k = 2, i.e. x n+1 = α + bx n x n 1 (2) with b 0. By the change of variable x = by we get the recurrence y n+1 = a + y n y n 1 where a = α/b 2, so that we are led to the recurrence in one parameter, which we rewrite again as x n+1 = a + x n x n 1 (3) In the following we shall refer to the recurrence (3) as the simplest recurrence. For a > 0 it has been proved that (see e.g. [Kocic et al., 1993; Kocic & Ladas, 1993]): (i) the solutions of (3) satisfy the equation H(x n 1, x n ) = const. where H(x n 1, x n ) ( = (a + x n 1 + x n ) ) ( ) x n 1 x n (4) (ii) if the initial conditions (i.c. henceforth) are positive, then the solution of the recurrence (3) is given by a positive and bounded sequence, i.e. there exist m > 0 and M > m such that m < x n < M n; (iii) a unique positive equilibrium exists and all the positive solutions are oscillating around the positive equilibrium. Of particular interest are the cases called periodic recurrences, in which all the trajectories are periodic. The following result holds: (iv) every positive solution of (3) is periodic of period 5 iff a = 1; every positive solution of (3) is periodic of period 6 iff a = 0. Moreover, in [Csornyei & Laczkovich, 2001] is proved that these are the only possibilities for the recurrence (2) to be a periodic recurrence. In other words, the infinitely many solutions of (2) are all periodic iff a = 0 or a = 1, giving rise to periodic solutions of periods 6 and 5, respectively. Indeed, as we shall see in the following, for the recurrence (2) such values of the parameter a can be seen as bifurcation values when we consider such recurrence from the point of view of the singularities specific to two-dimensional iterated maps of the plane which are not defined in the whole plane due to a vanishing denominator, as given in [Bischi et al., 1999a, 1999b, 2000]. Our first result concerns the properties of the recurrence (3) when a is a negative parameter. As we shall see, the equilibria of (3) exist only for a > 1/4, so we are interested in the parameter range 0.25 < a < 0. Our first result is expressed by the following Theorem: Theorem 1. Let 0.25 < a < 0 in recurrence (3), then
3 Lyness Iterative Processes 1843 (j) the solutions satisfy H(x n 1, x n ) = const., where H is defined in (4); (jj) there are two fixed points of the recurrence (3) given by S = (x s, x s), with x s = ( a)/2 and P = (x p, x p) with x p = ( a)/2, of saddle and centre types respectively, such that a positive region in the plane (x n 1, x n ) exists, bounded by a closed invariant curve, locus of homoclinic orbits of the saddle S, and including the equilibrium P, inside which all the trajectories are bounded and oscillating around P ; (jjj) the trajectories of i.c. outside the area described in (jj) belong to invariant phase curves which cross through the four points (0, a), ( 1, 0), ( a, 0), (0, 1) and are made up of at least three unbounded branches; (jv) there exists a three-cycle C 3 of saddle type, given by: 1, 1, (1 a), 1, 1, (1 a),... (v) in the plane (x n 1, x n ), the three lines (r i ), i = 1, 2, 3, of equation x n = 0, x n 1 + x n +a = 0, and x n +1 = 0 respectively, constitute the stable set of the three-cycle C 3 given in (jv), and their points generate, through (3), sequences for which an explicit analytic expression exists, given by: 1, u 0, (u 0 + a), 1, u 1, (u 1 + a),... where u n = (a 1)n (u a) (a 1)(u 0 + 1) (a 1) n (u a) + (u 0 + 1) u 0 a and n 0 In order to prove this theorem we shall use the terminology and some properties of the iterated maps of the plane having a vanishing denominator. In the next section we first recall such properties and then we give a proof of Theorem Proof of Theorem 1 In order to prove the Theorem stated in Sec. 2, let us rewrite the second-order recurrence in (3) as a two-dimensional system of the first order, i.e. an iterated map of the plane. As usual, this is obtained by letting (x n 1, x n ) = (x, y) so that (3) can be written as T : (x n 1, x n ) (x n, x n+1 ), where T (x, y) (x, y ) is the two-dimensional map given by x = y T : y = a + y (5) x 3.1. Properties of the equivalent two-dimensional map The map T is not defined on the line δ s of equation x = 0, so the iteration of T generates uninterrupted sequences provided that the initial condition (x 0, y 0 ) (i.c. henceforth) belongs to the set E given by E = R 2 \Λ (6) where Λ is the union of the preimages of any rank of the line δ s Λ = T k (δ s ) k=0 In other words, the sequences generated by the recurrence (3) can be obtained by the iteration of the two-dimensional map T : E E which is not defined in the whole plane. 1 Following the terminology introduced in [Bischi et al., 1999a] the line δ s will be called set of nondefinition, and its point Q = (0, a), where the second component of the map T assumes the form 0/0 constituting a simple focal point of T, and the associated prefocal set δ(q) is the line of equation x = a. This means that a one-to-one correspondence exists between the slopes m of arcs γ through the focal point Q and the points ( a, y) where their images T (γ) cross the prefocal set δ(q). In this case such correspondence is very simple, given by m y(m) = m This implies that the image by T of an arc γ crossing through Q with a slope m in Q, is an arc T (γ) which crosses the prefocal line at the point ( a, m), and, conversely, the preimage T 1 (η) of an arbitrary arc η which crosses the line x = a at a point ( a, y) is an arc which crosses through Q with a slope m = y in Q, where the inverse map is given by ( ) a + x T 1 (x, y) =, x (7) y The latter property justifies the terms focal point and prefocal set, since all the segments crossing 1 Notice that the set Λ of points excluded from the phase space of the iterated map has zero Lebesgue measure in R 2.
4 1844 L. Gardini et al. (a) (b) Fig. 1. the prefocal line are focalized by T 1 into arcs through Q [see the qualitative pictures in Fig. 1(a)]. Instead, if an arc γ crosses the set of nondefinition δ s at a nonfocal point (0, ξ), then its image by T is an unbounded arc, doubly asymptotic to the line of equation x = ξ (see [Bischi et al., 1999a, 1999b]). The property that infinitely many invariant curves of the map may cross through focal points of T or of its inverse T 1, although such points are not fixed points or cycles of the map T, is now well known. We shall see that this also occurs in this example, and the properties of focal points and prefocal sets of T and T 1 will help us to understand the peculiar properties of the recurrence (3) Proof Let us first consider the invariant curves of the map (5). We already know the existence of a constant function H(x, y) in the case a > 0, and it is straightforward to see that the same function plays the same role also in the case a < 0. That is, the phase curves of T satisfy H(x, y) = const. where ( H(x, y) = (a + x + y) ) ( ) (8) x y being H(T (x, y)) = H(x, y). Moreover, the map T is invertible and its inverse, given by (7), has invariant phase curves which satisfy the same equation and such curves are symmetric with respect to the line x = y. In fact, the two maps, T and T 1, are topologically conjugate by the symmetry s(x, y) = (y, x), with s T s = T 1. So, both T and T 1 have the same phase curves. Moreover, looking at the explicit definition (7) of T 1, it is easy to see that the line y = 0 is the set of nondefinition for T 1, we denote it by δ s, and its point Q = ( a, 0) is a focal point for T, the related prefocal set δ(q ) being the line of equation y = a. This means that whenever an arc, say γ, crosses through the line y = a, then T 1 (T 1 (γ)) = T (γ) is an arc which crosses through the focal point Q. And if an arc crosses the set of nondefinition of T 1, y = 0, in a nonfocal point (ξ, 0), ξ 0, then its image by T 1 is an unbounded arc, doubly asymptotic to the line of equation y = ξ [see Fig. 1(b)]. All these facts lead us to the properties of the phase curves of T described in Theorem 1. By elementary calculus we can see that for 0.25 < a < 0 there exist two fixed points of T, both in the
5 Lyness Iterative Processes 1845 positive half-plane, a centre P and a saddle S : P = (x p, x p), x p = a 2 S = (x s, x s ), x s = a. 2 The two eigenvalues, and related eigenvectors, of the Jacobian matrix of T computed in the saddle S are given by λ 1 = 1 1 4x 2 s 2x, r 1 = (1, λ 1 ) s λ 2 = 1 λ 1, r 2 = (1, λ 2 ) Two branches of the invariant curves issuing from the saddle S must merge into a closed invariant curve, locus of homoclinic orbits of S, in the positive half-plane. This is a consequence of the fact that the phase curves, defined by H(x, y) = const., are symmetric with respect to the line x = y. In fact, considering H(x, y) = H(x s, x s) we show that the curve must cross the line x = y at a point (p s, p s ), with x s < x p < p s. This can be proved by analyzing the properties of the one-dimensional restriction of H to the line x = y, given by the function F (x) = H(x, x), which reads ( F (x) = (a + 2x) ) 2 x The qualitative graph of F (x) (for 0.25 < a < 0) is shown in Fig. 2, and from the properties of its local extrema we prove the existence of p s, which means that the phase curves of T through the saddle S have homoclinic orbits on a closed curve through the fixed point [see Fig. 3(a)]. But in the graph of F (x), shown in Fig. 2, besides the existence of the saddle S (the local maximum in x s ) and of the center P (the local minimum in x p), another critical point of the function H(x, y) can be seen, given by the local maximum at the point x = 1. This implies that the two-dimensional T must have a saddle also at the point ( 1, 1). However, that point is not a fixed one of T, but a point of a three-cycle of T of saddle type, as described below. So, we have proved that homoclinic orbits of S exist on a closed curve through S, which crosses the line x = y at the point (p s, p s ) and consequently bounds a finite region in the positive quadrant of the plane. And it is clear that any i.c. inside this area belongs to a phase curve which is a closed curve Fig. 2. around the center P. This concludes the proof of parts (j) and (jj) of Theorem 1. To prove part (jjj) let us consider the two other branches of phase curve issuing from the saddle S. Being x s > a the saddle S is located above the two prefocal lines δ(q) and δ(q ) of T and T 1 respectively. Thus, one branch issuing from S must cross δ(q ) in a point [see the arc γ in Fig. 3(a)] and consequently its image T (γ) must cross Q, its image T 2 (γ) must cross through T (Q ) = (0, 1) δ s and, as a consequence, T 3 (γ) must be an unbounded arc, doubly asymptotic to the line of equation y = 1. A similar reasoning applies to the other branch of phase curve issuing from S (symmetric of the previous ones). That branch includes an arc η which crosses δ(q), so that T 1 (η) crosses through the focal point Q, T 2 (η) crosses the point T 1 (Q) = ( 1, 0) δ s, which means that T 3 (η) is unbounded and doubly asymptotic to the line of equation x = 1. These unbounded branches of the phase curve through S are shown in Fig. 3(b). Of course, other phase curves must cross the lines δ(q) and δ(q ), thus giving those with unbounded arcs the same properties as the phase curves through S as described above. This is true for the sequences generated starting from any point of the positive half-plane external to the closed curve of homoclinic orbits of S. For example, let us consider an initial condition taken at the points (q, q) with q > p s ; an invariant curve through (q, q) must cross the prefocal δ(q ) at a point, say (ζ, a), and the prefocal δ(q) in the symmetric point ( a, ζ). Indeed, such curves fill up the plane
6 1846 L. Gardini et al. Fig. 3. Fig. 4. (x, y) (just a few of them are shown in Fig. 4). By the above arguments we have proved that all these phase curves cross through the four points Q = (0, a), T 1 (Q) = ( 1, 0), (9) Q = ( a, 0), T (Q ) = (0, 1) (10) and are asymptotic to the lines of equation x = 1 and y = 1. But more, we now show that these phase curves are also asymptotic to another line, of equation x + y + c = 0, where c a, so that these invariant curves must include at least three unbounded branches, thus proving part (jjj) of Theorem 1. This is a consequence of the fact that each
7 Lyness Iterative Processes 1847 of the following lines: (r 1 ) : x = 1, (r 2 ) : x + y + a = 0, (r 3 ) : y = 1 (11) is mapped into itself by the third iterate of the map, T 3. In fact, let ( 1, u) be a point belonging to (r 1 ), then T ( 1, u) = (u, a u) (r 2 ), T 2 ( 1, u) = ( a u, 1) (r 3 ), and T 3 ( 1, u) = [ 1, ((1 a)/(a+u))] (r 1 ). Moreover, the restriction of T 3 to any one of these lines is represented by the one-dimensional map: u = f(u), f(u) = 1 a (12) a + u The function f(u) has very simple properties (the graph of f(u) is an hyperbola with vertical asymptote in u = a). There are two fixed points of f(u), u = 1 (attracting for a < 0 and repelling for a > 0, with f ( 1) = 1/(a 1)), and u = 1 a (attracting for a > 0, repelling for a < 0, with f (1 a) = a 1). The stable fixed point is globally attracting for f(u). Moreover, by using the function h(z) = (z + 1)/(z (1 a)) we have that f(u) is conjugate to the function g(y), where g(y) = h f h 1 (y) = ( ) 1 y a 1 is linear, so that we can easily write the analytical solutions of the iterates of g(y), and thus of f(u), obtaining, u 0 a, the following trajectory: u n =f n (u 0 )= (a 1)n (u 0 1+a) (a 1)(u 0 +1) (a 1) n (u 0 1+a)+(u 0 +1) It is also straightforward to see that the fixed points of f(u) belong to a cycle of period 3 of T, say C 3, given by: ( 1, 1 a), (1 a, 1), ( 1, 1) (13) and the eigenvalues of T 3 on each periodic point are given by λ 1 = (a 1), λ 2 = 1 λ 1 so that the three-cycle of T, C 3, is of saddle type, its stable set being constituted by the three lines, cyclical for T (and invariant for T 3 ), given in (11). These three lines constitute the phase curve associated with H( 1, 1) (or with any other point on the lines), so that they cannot be crossed by any other invariant phase curve of the plane, except at the points of nondefinition of T and T 1, that is, at the intersections of the lines δ s, of equation x = 0, and δ s, of equation y = 0, with the three cyclical lines in (11), i.e. the four points given in (9). As a consequence, the unbounded arcs which we have proved to be asymptotic to the lines x = 1 and y = 1 are also doubly asymptotic to x + y + a = 0. Fig. 5.
8 1848 L. Gardini et al. Fig. 6. or to one line parallel to it, equation x + y + c = 0, with c a. In order to see the qualitative shape of all the phase curves outside the homoclinic orbit of S, it is enough to consider an i.c. (x, x) on varying x (x < x s and x > p s ) or i.c. ( a, ξ) on δ(q), or (ξ, a) on δ(q ), and varying ξ. As these invariant curves must intersect the prefocal lines δ(q) and δ(q ), the same reasoning as above apply, and we can conclude that they all cross the four points given in (9), are unbounded, and made up of at least three pieces, doubly asymptotic to the lines x = 1, y = 1 and x + y + c = 0 for some c a (their qualitative shape is shown in Figs. 5 and 6). This completes the proof of Theorem 1. We close this section by stressing that the sequences generated by points belonging to unbounded phase curves are unbounded but not diverging. This means that if we consider an i.c. {x 0, y 0 } outside the homoclinic loop of S, then the sequence {x n, n 0}, is such that sup{ x n, n 0} = + but lim n x n does not exist. In other words, these sequences arbitrarily approach the Poincaré Equator (i.e. the points at infinite distance) along the cyclical lines given in (11), but further iteration of T gives points which are at a smaller distance from the origin, a kind of motion which is frequently observed in iterated maps with denominator, as shown in [Bischi et al., 2000]. 4. Some Results on the Simplest Recurrence with a > 0 In the previous sections we have shown that for 0.25 < a < 0 the phase curves of (3) are invariant closed curves around P only in a suitable region, and we can refer to this region by using the points (x, x), with x (x s, p s ), as defined in the previous section. As the parameter a 0, the fixed point S (as well as Q and Q ) approaches the origin, so that the portion of the positive orthant made up of trajectories oscillating around the center (bounded by the homoclinic orbit of S ) become wider and wider as a 0. A qualitative picture is shown in Fig. 7. It is easy to see that the properties of T, and of the phase curves, as described in Sec. 3 are not due to the negativity of the parameter a. That is, those properties continue to hold also for a > 0, with obvious changes in part (jj) of Theorem 1. So, starting from our Theorem 1 we can extend the results given in [Kocic et al., 1993; Kocic & Ladas, 1993], in the sense that we can assume the statements (ii) and (iii) that were recalled at the beginning of Sec. 2 are not only sufficient but also necessary. Indeed, we have the following Theorem 2. Consider (5) with a > 0, then (j) an i.c. (x, y) gives rise to a positive and bounded trajectory (oscillating around P ) iff (x, y) R 2 +.
9 Lyness Iterative Processes 1849 Fig. 7. (jj) the trajectory of an i.c. (x, y) / R+ 2 belongs to an invariant phase curve which crosses through the four points Q = (0, a), T 1 (Q) = ( 1, 0), Q = ( a, 0), T (Q ) = (0, 1) and includes at least three unbounded branches doubly asymptotic to the lines x = 1 ; x + y + c = 0 (c a) ; y = 1 (jjj) for each a 1 there exists a three-cycle C 3 of saddle type, given by: 1, 1, (1 a), 1, 1, (1 a),... (jv) in the plane (x n 1, x n ), the three lines (r i ), i = 1, 2, 3, of equation x n = 0, x n 1 + x n +a = 0, and x n +1 = 0 respectively, constitute the stable set of the three-cycle C 3 given in Theorem 1, and their points generate, through (3), sequences for which an explicit analytic expression exists, given by: 1, u 0, (u 0 + a), 1, u 1, (u 1 + a),... (14) where u n = (a 1)n (u a) (a 1)(u 0 + 1) (a 1) n (u a) + (u 0 + 1) u 0 a and n 0., From the structure of the phase curves and the properties seen above and in the previous sections, it follows that the cases a = 0 and a = 1, which according to [Kocic et al., 1993; Csornyei & Laczkovich, 2001] are associated with periodic recurrences, correspond to bifurcation cases related to contacts between singularities of different nature of the map T. In fact, for a = 0 the focal points Q and Q merge, and both merge with the saddle fixed point S. Moreover, the lines of nondefinition of T and T 1, given by δ s and δ s respectively, merge with the two prefocal lines δ(q) and δ(q ), respectively. The related changes in the shape of the phase curves can be seen by comparing Fig. 7 (obtained for a < 0) and Fig. 8 (obtained for a > 0). Instead, for a = 1 the focal points Q and Q merge with the periodic points of the degenerated three-cycle (14), and the prefocal lines δ(q) and δ(q ) merge with the two degenerated cyclical lines given in (11). We note that for a = 1 the i.c. which belong to the lines (11) are not included in the set E defined in (6) whose points generate uninterrupted sequences, because such points are mapped into the point Q = (0, 1) after a finite number of iterations of T. So, these points do not generate periodic solutions of period 5. We also note that in [Kocic et al., 1993; Kocic & Ladas, 1993] it was proved that for a = 0 the recurrence is periodic, because all the i.c. in the set E generate periodic
10 1850 L. Gardini et al. Fig. 8. sequences of period 6. However, in order to avoid misunderstanding, it is better to specify that this statement was not to be intended as minimum period 6, because, as we have seen in (14), a cycle of period 3 exists, given by 1, 1, 1, 1, 1, 1, Extensions In the previous sections we have considered the simplest case of the family of recurrences given in (1), obtained for k = 2. However, as stated in the introduction, it is natural to conjecture that similar properties continue to hold, with obvious changes, even for any k > 2. Let us consider the recurrence (1) of order k when all the parameters b i are equal, say b = b 0 = b 1 = = b k 2, b 0, so that, after a change of variable, the recurrence may be rewritten as x n+1 = a + x n (k 2) + + x n 1 + x n x n (k 1) (15) where a = α/b 2. Clearly the recurrence of order k in (15) is equivalent to a system of k equations of the first order. Identifying (x n (k 1),..., x n 1, x n ) = (y 1,..., y k 1, y k ) we get a k- dimensional map T (k) : R k R k, (y 1,..., y k 1, y k ) = T (k)(y 1,..., y k 1, y k ) where: y 1 = y 2. y k 1 = y k y k = a + y y k y 1 (16) Some general results on the generic recurrence (15) are already given in [Kocic et al., 1993; Kocic & Ladas, 1993] for the case a > 0, such as: the function H(x n (k 1),..., x n 1, x n ) k 1 = a + j=0 x n j k 1 j=0 ( ) x n j (17) is constant along the solutions of (15); every positive solution of (15) is oscillating about the (unique) positive equilibrium P. Moreover, concerning the special case of periodic recurrences, in [Kocic & Ladas, 1993] it is shown that for k = 3 and a = 1 the recurrence (15) is periodic of period 8, and in [Csornyei & Laczkovich, 2001] it is proved that for k = 3 the same recurrence is periodic iff a = 1 whereas for k > 3 it cannot be periodic.
11 Lyness Iterative Processes 1851 But more can be proved on the global dynamics of the recurrence (15) in the regions of the phase space having some negative components, or when the parameter a is negative. We state the following Theorem, which holds for the recurrence (15) independent of the sign of a, and which will allow us to obtain some new insight into the special cases of periodic recurrences associated with k = 3. Theorem 3. Consider the recurrence (15), or equivalently the map T (k) in (16), for any k 2. Then (a) The function given in (17) is constant along the solutions of the recurrence; (b) for any a k 1 there exists a periodic solution of (15) of period (k + 1) given by 1,..., 1, (k 1 a),... }{{} (k times)... 1,..., 1, (k 1 a),... }{{} (k times) (18) (c) there exist (k+1) cyclical hyperplanes Π i of R k, mapped one into the other by the map T (k), invariant for the map T k+1 (k), given, in their cyclical order, by: (Π 1 ) : y = 0 (Π 2 ) : a + y y k = 0 (Π 3 ) : y k + 1 = 0 (Π 4 ) : y k = 0. (Π k+1 ) : y = 0 Proof. The statement in (a) is straightforward, being H(T (k) (y 1,..., y k 1, y k )) = H(y 1,..., y k 1, y k ). To prove the existence of the periodic orbit of period (k + 1) let us start with the i.c. made up of k times the value 1 : (x n (k 1),..., x n 1, x n ) = ( 1,..., 1). Then a (k 1) x n+1 = = k 1 a 1 and a + (k 1 a) (k 1 j) (j 1) x n+1+j = 1 = 1 for j = 1,..., k, so that x n+1+(k+1) = k 1 a and we have the (k + 1)-cycle in (18). To prove (c) let us consider a point p Π 1, then p = ( 1, z 2, z 3,..., z k ) Π 1 T (k) (p) = (z 2, z 3,..., z k, (a+z 2 +z 3 + +z k )) Π 2 its image T 1+1 (k) (p) = (z 3,..., z k, (a + z 2 + z z k ), 1) Π 3 has a 1 in the last position, so that T 1+j (k) (p) has a 1 in the position (k (j 1)), and T 1+k (k) (p) = ( 1,... ) Π 1 This complete the proof of Theorem 3. We observe that from the construction given above it is possible to also get the explicit formulation of the (k 1)-dimensional map which represents the restriction of T 1+k (k) to each of the cyclical hyperplanes. The case k = 2 was already given in the previous sections, with T = T (2), so, as a further example, let us consider now the case k = The Recurrence with k = 3 and the Equivalent 3-D Map In this section we consider the Lyness recurrence (1) with k = 3, which reads x n+1 = a + x n 1 + x n x n 2 (19) and can be rewritten as an iterated threedimensional map T (3) : (x n 2, x n 1, x n ) (x n 1, x n, x n+1 ), where T (3) (y 1, y 2, y 3 ) (y 1, y 2, y 3 ) is given by: y 1 = y 2 y 2 = y 3 T (3) : y 3 = a + y 2 + y 3 y 1 (20) where, as usual, the variables (y 1, y 2, y 3 ) are defined by (y 1, y 2, y 3 ) = (x n 2, x n 1, x n ). As stated in Theorem 3, four cyclical planes of R 3 exist on which the fourth-iterate of the map (20), T(3) 4, can be identified with a two-dimensional map. After a straightforward computation, it is easy to see that the
12 1852 L. Gardini et al. explicit formulation of this two-dimensional map Z : (u, v) (u, v ) = Z(u, v) is given by: Z : v = u = 1 + u + v v v(1 a) + (1 + u + v) v(a + u + v) (21) This explicit relation allows us to see that the periodic points of T (3) of period 4 are much more than those described in Theorem 3 (b). Indeed, the twodimensional map Z in (21) has a line of fixed points, which are necessarily periodic points of period 4 for T (3). They are given, u 1 a, by: ( 1, u, 1), (u, 1, 1 a u), ( 1, 1 a u, 1), (1 a u, 1, u) Stated with other words, for any value u 1 a the sequence 1, u, 1, 1 a u, 1, u, 1, 1 a u,... is one of the infinitely many (uncountable) solutions of period 4 of the recurrence (15) with k = 3. Notice that this holds for any value of a, in particular for a = 1, at which the recurrence (19) is periodic with period 8, as proved in [Kocic & Ladas, 1993; Csornyei & Laczkovich, 2001]. So, like in the previous section, our results help to clarify that the term periodic of period 8 is not to be intended as minimum period 8 since infinitely many periodic sequences of period 4 exist. This confirms our previous remark given at the end of Sec. 4. We close this section by stating the results given above in the form of the following Theorem: Theorem 4. Consider the recurrence (19). For any value of a infinitely many periodic sequences of period 4 can be generated, given by 1, u, 1, 1 a u, 1, u, 1, 1 a u, 1,... for each u 1 a. Acknowledgments We thank the seminar participants of the European Conference on Iteration Theory, ECIT-2000, for their comments. This work has been performed under the auspices of CNR, Italy, GNFM, INDAM, Italy, and under the activity of the national research project Dynamic Models in Economics and Finance: Evolution, Uncertainty and Forecasting, MURST, Italy. References Bischi, G. I. & Gardini, L. [1997] Basin fractalization due to focal points in a class of triangular maps, Int. J. Bifurcation and Chaos 7, Bischi, G. I., Gardini, L. & Mira, C. [1999a] Maps with denominator. Part 1: Some generic properties, Int. J. Bifurcation and Chaos 9, Bischi, G. I., Gardini, L. & Mira, C. [1999b] New phenomena related to the presence of focal points in twodimensional maps, Ann. Math. Silesianae (special issue Proc. ECIT98) 13, Bischi, G. I., Mira, C. & Gardini, L. [2000] Unbounded sets of attraction, Int. J. Bifurcation and Chaos 10, Csornyei, M. & Laczkovich, M. [2001] Some periodic and non-periodic recursions, Monatshefte für Mathematik 132, Kocic, V. I. & Ladas, G. [1993] Global Behaviour of Nonlinear Difference Equations of Higher Order with Applications (Kluwer Academic Publishers, the Netherlands). Kocic, V. I., Ladas, G. & Rodrigues, I. W. [1993] On rational recursive sequences, J. Math. Anal. Appl. 173, Lin, J. & Liu, Z. [1999] Invariant curves of the generalized Lynass equations, Int. J. Bifurcation and Chaos 9, Lyness, R. C. [1942] Note 1581, The Mathematical Gazette 26, p. 62. Lyness, R. C. [1945] Note 1847, The Mathematical Gazette 29, Lyness, R. C. [1961] Note 2952, The Mathematical Gazette 45, Mira, C. [1999] Some properties of two-dimensional (Z 0 Z 2 ) maps not defined in the whole plane, Grazer Math. Berichte (special issue Proc. ECIT96) 339,
UNBOUNDED SETS OF ATTRACTION
International Journal of Bifurcation and Chaos, Vol. 10, No. 6 (2000) 1437 1469 c World Scientific Publishing Company UNBOUNDED SETS OF ATTRACTION GIAN-ITALO BISCHI Istituto di Scienze Economiche, University
More informationRoutes to Complexity in a Macroeconomic Model Described by a Noninvertible Triangular Map
Cubo A Mathematical Journal Vol.05/N ō 03 OCTOBER 2003 Routes to Complexity in a Macroeconomic Model Described by a Noninvertible Triangular Map Roberto Dieci Dipartimento di Studi Economici e Quantitativi,
More informationMaps with vanishing denominator explained through applications in Economics
Journal of Physics: Conference Series PAPER OPEN ACCESS Maps with vanishing denominator explained through applications in Economics To cite this article: Fabio Tramontana 2016 J. Phys.: Conf. Ser. 692
More informationMultistability and cyclic attractors in duopoly games
Chaos, Solitons and Fractals 11 (2000) 543±564 www.elsevier.nl/locate/chaos Multistability and cyclic attractors in duopoly games Gian Italo Bischi a, Cristiana Mammana b, Laura Gardini a,c, * a Istituto
More informationHomoclinic tangles associated with closed invariant curves in families of 2D maps
ITERATION THEORY (ECIT 04) W. Förg-Rob, L. Gardini, D. Gronau, L. Reich, J. Smítal (Eds.) Grazer Math. Ber., ISSN 1016 7692 Bericht Nr. 350 (2006), 1 14 Homoclinic tangles associated with closed invariant
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More informationThe Fibonacci sequence modulo π, chaos and some rational recursive equations
J. Math. Anal. Appl. 310 (2005) 506 517 www.elsevier.com/locate/jmaa The Fibonacci sequence modulo π chaos and some rational recursive equations Mohamed Ben H. Rhouma Department of Mathematics and Statistics
More informationSTABILITY OF THE kth ORDER LYNESS EQUATION WITH A PERIOD-k COEFFICIENT
International Journal of Bifurcation Chaos Vol 7 No 2007 43 52 c World Scientific Publishing Company STABILITY OF THE kth ORDER LYNESS EQUATION WITH A PERIOD-k COEFFICIENT E J JANOWSKI M R S KULENOVIĆ
More informationCenter Bifurcation of a point on the Poincaré Equator
ITERATION THEORY (ECIT 08) A.N. Sharkovsky, I.M. Sushko (Eds.) Grazer Math. Ber., ISSN 1016 7692 Bericht Nr. 354 (2009), 254-279 Center Bifurcation of a point on the Poincaré Equator Iryna Sushko and Laura
More informationTransitions to Chaos in a 3-Dimensional Map
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 12, Number 5 (2016), pp. 4597 4611 Research India Publications http://www.ripublication.com/gjpam.htm Transitions to Chaos in a 3-Dimensional
More informationINTRICATE ASSET PRICE
Chapter 1 INTRICATE ASSET PRICE DYNAMICS AND ONE-DIMENSIONAL DISCONTINUOUS MAPS F. Tramontana, L. Gardini and F. Westerhoff * Department of Economics and Quantitative Methods, University of Urbino, Via
More informationBASINS OF ATTRACTION FOR TRIPLE LOGISTIC MAPS. 1. Introduction
BASINS OF ATTRACTION FOR TRIPLE LOGISTIC MAPS JAN ANDRES AND JIŘÍ FIŠER Abstract. Attractors and the associated basins of attraction for triple logistic maps are detected and visualized numerically (see
More informationCoupled Chaotic Fluctuations in a Model of International Trade and Innovation: Some Preliminary Results
Coupled Chaotic Fluctuations in a Model of International Trade and Innovation: Some Preliminary Results Iryna Sushko a,b, Laura Gardini c, Kiminori Matsuyama d a Institute of Mathematics, NASU, b Kyiv
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationA dynamic marketing model with best reply and inertia
A dynamic marketing model with best reply and inertia Gian Italo Bischi - Lorenzo Cerboni Baiardi University of Urbino (Italy) URL: www.mdef.it/gian-italo-bischi/ Urbino, September 18, 14 Urbino, September
More informationDynamics of the equation complex plane
APPLIED & INTERDISCIPLINARY MATHEMATICS RESEARCH ARTICLE Dynamics of the equation in the complex plane Sk. Sarif Hassan 1 and Esha Chatterjee * Received: 0 April 015 Accepted: 08 November 015 Published:
More informationGlobal Attractivity of a Higher-Order Nonlinear Difference Equation
International Journal of Difference Equations ISSN 0973-6069, Volume 5, Number 1, pp. 95 101 (010) http://campus.mst.edu/ijde Global Attractivity of a Higher-Order Nonlinear Difference Equation Xiu-Mei
More informationNonhyperbolic Dynamics for Competitive Systems in the Plane and Global Period-doubling Bifurcations
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 3, Number 2, pp. 229 249 (2008) http://campus.mst.edu/adsa Nonhyperbolic Dynamics for Competitive Systems in the Plane and Global Period-doubling
More informationGlobal Attractivity in a Higher Order Nonlinear Difference Equation
Applied Mathematics E-Notes, (00), 51-58 c ISSN 1607-510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Global Attractivity in a Higher Order Nonlinear Difference Equation Xing-Xue
More informationGlobal Attractivity in a Nonlinear Difference Equation and Applications to a Biological Model
International Journal of Difference Equations ISSN 0973-6069, Volume 9, Number 2, pp. 233 242 (204) http://campus.mst.edu/ijde Global Attractivity in a Nonlinear Difference Equation and Applications to
More informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More informationSUBCONTINUA OF FIBONACCI-LIKE INVERSE LIMIT SPACES.
SUBCONTINUA OF FIBONACCI-LIKE INVERSE LIMIT SPACES. H. BRUIN Abstract. We study the subcontinua of inverse limit spaces of Fibonacci-like unimodal maps. Under certain combinatorial constraints no other
More informationBifurcations of Traveling Wave Solutions for a Generalized Camassa-Holm Equation
Computational and Applied Mathematics Journal 2017; 3(6): 52-59 http://www.aascit.org/journal/camj ISSN: 2381-1218 (Print); ISSN: 2381-1226 (Online) Bifurcations of Traveling Wave Solutions for a Generalized
More informationIntegrability and dynamical degree of periodic non-autonomous Lyness recurrences. Anna Cima
Proceedings of the International Workshop Future Directions in Difference Equations. June 13-17, 2011, Vigo, Spain. PAGES 77 84 Integrability and dynamical degree of periodic non-autonomous Lyness recurrences.
More informationMATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation
MATH 415, WEEK 11: Bifurcations in Multiple Dimensions, Hopf Bifurcation 1 Bifurcations in Multiple Dimensions When we were considering one-dimensional systems, we saw that subtle changes in parameter
More informationOn the smoothness of the conjugacy between circle maps with a break
On the smoothness of the conjugacy between circle maps with a break Konstantin Khanin and Saša Kocić 2 Department of Mathematics, University of Toronto, Toronto, ON, Canada M5S 2E4 2 Department of Mathematics,
More informationSingular Perturbations in the McMullen Domain
Singular Perturbations in the McMullen Domain Robert L. Devaney Sebastian M. Marotta Department of Mathematics Boston University January 5, 2008 Abstract In this paper we study the dynamics of the family
More informationarxiv: v1 [nlin.cd] 20 Jul 2010
Invariant manifolds of the Bonhoeffer-van der Pol oscillator arxiv:1007.3375v1 [nlin.cd] 20 Jul 2010 R. Benítez 1, V. J. Bolós 2 1 Departamento de Matemáticas, Centro Universitario de Plasencia, Universidad
More informationA MINIMAL 2-D QUADRATIC MAP WITH QUASI-PERIODIC ROUTE TO CHAOS
International Journal of Bifurcation and Chaos, Vol. 18, No. 5 (2008) 1567 1577 c World Scientific Publishing Company A MINIMAL 2-D QUADRATIC MAP WITH QUASI-PERIODIC ROUTE TO CHAOS ZERAOULIA ELHADJ Department
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationPeriodicities in Linear Fractional Recurrences: Degree growth of birational surface maps. Eric Bedford* and Kyounghee Kim
Periodicities in Linear Fractional Recurrences: Degree growth of birational surface maps Eric Bedford* and Kyounghee Kim revision: August 23, 2006 0. Introduction Given complex numbers α 0,...,α p and
More informationBIFURCATION ANALYSIS OF A CIRCUIT-RELATED GENERALIZATION OF THE SHIPMAP
International Journal of Bifurcation and Chaos, Vol. 16, No. 8 (26) 2435 2452 c World Scientific Publishing Company BIFURCATION ANALYSIS OF A CIRCUIT-RELATED GENERALIZATION OF THE SHIPMAP FEDERICO BIZZARRI
More informationM.ARCIERO, G.LADAS AND S.W.SCHULTZ
Georgian Mathematical Journal 1(1994), No. 3, 229-233 SOME OPEN PROBLEMS ABOUT THE SOLUTIONS OF THE DELAY DIFFERENCE EQUATION x n+1 = A/x 2 n + 1/x p n k M.ARCIERO, G.LADAS AND S.W.SCHULTZ Abstract. We
More informationInvariant manifolds of the Bonhoeffer-van der Pol oscillator
Invariant manifolds of the Bonhoeffer-van der Pol oscillator R. Benítez 1, V. J. Bolós 2 1 Dpto. Matemáticas, Centro Universitario de Plasencia, Universidad de Extremadura. Avda. Virgen del Puerto 2. 10600,
More informationOn Periodic points of area preserving torus homeomorphisms
On Periodic points of area preserving torus homeomorphisms Fábio Armando Tal and Salvador Addas-Zanata Instituto de Matemática e Estatística Universidade de São Paulo Rua do Matão 11, Cidade Universitária,
More informationSYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 21, 2003, 211 226 SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS Massimo Grosi Filomena Pacella S.
More informationAttractivity of the Recursive Sequence x n+1 = (α βx n 1 )F (x n )
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.7(2009) No.2,pp.201-206 Attractivity of the Recursive Sequence (α βx n 1 )F (x n ) A. M. Ahmed 1,, Alaa E. Hamza
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More informationLinear Fractional Recurrences: Periodicities and Integrability. Eric Bedford and Kyounghee Kim
Linear Fractional Recurrences: Periodicities and Integrability Eric Bedford and Kyounghee Kim 0. Introduction. Let k 3, and let α = (α 0,..., α k ), β = (β 0,..., β k ) be (k + 1)-tuples of complex numbers.
More informationDynamics of a Rational Recursive Sequence
International Journal of Difference Equations ISSN 0973-6069, Volume 4, Number 2, pp 185 200 (2009) http://campusmstedu/ijde Dynamics of a Rational Recursive Sequence E M Elsayed Mansoura University Department
More informationBehaviour of Lipschitz functions on negligible sets. Non-differentiability in R. Outline
Behaviour of Lipschitz functions on negligible sets G. Alberti 1 M. Csörnyei 2 D. Preiss 3 1 Università di Pisa 2 University College London 3 University of Warwick Lars Ahlfors Centennial Celebration Helsinki,
More informationFundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
More informationTHE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL
THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL HAL L. SMITH* SCHOOL OF MATHEMATICAL AND STATISTICAL SCIENCES ARIZONA STATE UNIVERSITY TEMPE, AZ, USA 8587 Abstract. This is intended as lecture notes for nd
More informationFOURTH ORDER CONSERVATIVE TWIST SYSTEMS: SIMPLE CLOSED CHARACTERISTICS
FOURTH ORDER CONSERVATIVE TWIST SYSTEMS: SIMPLE CLOSED CHARACTERISTICS J.B. VAN DEN BERG AND R.C.A.M. VANDERVORST ABSTRACT. On the energy manifolds of fourth order conservative systems closed characteristics
More informationTWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2
More informationA Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term
ETASR - Engineering, Technology & Applied Science Research Vol., o.,, 9-5 9 A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term Fei Yu College of Information Science
More informationDynamical Systems: Ecological Modeling
Dynamical Systems: Ecological Modeling G Söderbacka Abstract Ecological modeling is becoming increasingly more important for modern engineers. The mathematical language of dynamical systems has been applied
More informationConnection between bifurcations on the Poincaré Equator and the dangerous bifurcations
ITERATION THEORY (ECIT 08) A.N. Sharkovsky, I.M. Sushko (Eds.) Grazer Math. Ber., ISSN 1016 7692 Bericht Nr. 354 (2009), 53-72 Connection between bifurcations on the Poincaré Equator and the dangerous
More informationSufficient conditions for a period incrementing big bang bifurcation in one-dimensional maps.
Sufficient conditions for a period incrementing big bang bifurcation in one-dimensional maps. V. Avrutin, A. Granados and M. Schanz Abstract Typically, big bang bifurcation occur for one (or higher)-dimensional
More informationCenter Bifurcation for Two-Dimensional Border-Collision Normal Form
Center Bifurcation for Two-Dimensional Border-Collision Normal Form Iryna Sushko Institute of Mathematics National Academy of Sciences of Ukraine, Kiev, Ukraine, sushko@imath.kiev.ua Laura Gardini Department
More informationModule 2: Reflecting on One s Problems
MATH55 Module : Reflecting on One s Problems Main Math concepts: Translations, Reflections, Graphs of Equations, Symmetry Auxiliary ideas: Working with quadratics, Mobius maps, Calculus, Inverses I. Transformations
More informationCHAOTIC BEHAVIOR IN A TWO-DIMENSIONAL BUSINESS CYCLE MODEL
CHAOTIC BEHAVIOR IN A TWO-DIMENSIONAL BUSINESS CYCLE MODEL CRISTINA JANUÁRIO Department of Chemistry, Mathematics Unit, Instituto Superior de Engenharia de Lisboa, Rua Conselheiro Emídio Navarro, 1949-014
More informationAn Analysis of Nondifferentiable Models of and DPCM Systems From the Perspective of Noninvertible Map Theory
An Analysis of Nondifferentiable Models of and DPCM Systems From the Perspective of Noninvertible Map Theory INA TARALOVA-ROUX AND ORLA FEELY Department of Electronic and Electrical Engineering University
More information16 1 Basic Facts from Functional Analysis and Banach Lattices
16 1 Basic Facts from Functional Analysis and Banach Lattices 1.2.3 Banach Steinhaus Theorem Another fundamental theorem of functional analysis is the Banach Steinhaus theorem, or the Uniform Boundedness
More informationSimplest Chaotic Flows with Involutional Symmetries
International Journal of Bifurcation and Chaos, Vol. 24, No. 1 (2014) 1450009 (9 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414500096 Simplest Chaotic Flows with Involutional Symmetries
More information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
More informationConnecting Orbits with Bifurcating End Points
Connecting Orbits with Bifurcating End Points Thorsten Hüls Fakultät für Mathematik Universität Bielefeld Postfach 100131, 33501 Bielefeld Germany huels@mathematik.uni-bielefeld.de June 24, 2004 Abstract
More informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More informationReceived January 27, 2008; Revised March 24, 2009
Papers International Journal of Bifurcation and Chaos, Vol. 19, No. 10 (2009) 3235 3282 c World Scientific Publishing Company FROM THE BOX-WITHIN-A-BOX BIFURCATION STRUCTURE TO THE JULIA SET. PART II:
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: January 08, 2018, at 08 30 12 30 Johanneberg Kristian
More informationDYNAMICS OF A DISCRETE BRUSSELATOR MODEL: ESCAPE TO INFINITY AND JULIA SET
DYNAMICS OF A DISCETE BUSSELATO MODEL: ESCAPE TO INFINITY AND JULIA SET HUNSEOK KANG AND YAKOV PESIN Abstract. We consider a discrete version of the Brusselator Model of the famous Belousov-Zhabotinsky
More informationAPPPHYS217 Tuesday 25 May 2010
APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag
More information11 Chaos in Continuous Dynamical Systems.
11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization
More informationOn the fixed points set of differential systems reversibilities arxiv: v1 [math.ds] 6 Oct 2015
On the fixed points set of differential systems reversibilities arxiv:1510.01464v1 [math.ds] 6 Oct 2015 Marco Sabatini October 5, 2015 Abstract We extend a result proved in [7] for mirror symmetries of
More informationShapley Polygons in 4 4 Games
Games 2010, 1, 189-220; doi:10.3390/g1030189 OPEN ACCESS games ISSN 2073-4336 www.mdpi.com/journal/games Article Shapley Polygons in 4 4 Games Martin Hahn Department of Mathematics, University Vienna,
More informationSpeculative behaviour and complex asset price dynamics: a global analysis
Journal of Economic Behavior & Organization Vol. 49 (2002) 173 197 Speculative behaviour and complex asset price dynamics: a global analysis Carl Chiarella a,, Roberto Dieci b, Laura Gardini c a School
More informationHomoclinic bifurcations in Chua s circuit
Physica A 262 (1999) 144 152 Homoclinic bifurcations in Chua s circuit Sandra Kahan, Anibal C. Sicardi-Schino Instituto de Fsica, Universidad de la Republica, C.C. 30, C.P. 11 000, Montevideo, Uruguay
More informationPERIODIC POINTS OF THE FAMILY OF TENT MAPS
PERIODIC POINTS OF THE FAMILY OF TENT MAPS ROBERTO HASFURA-B. AND PHILLIP LYNCH 1. INTRODUCTION. Of interest in this article is the dynamical behavior of the one-parameter family of maps T (x) = (1/2 x
More informationPropagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R
Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract We consider semilinear
More informationGlobal Attractivity, Oscillations and Chaos in A Class of Nonlinear, Second Order Difference Equations
Global Attractivity, Oscillations and Chaos in A Class of Nonlinear, Second Order Difference Equations Hassan Sedaghat Department of Mathematics, Virginia Commonwealth University Richmond, Virginia, 23284-2014,
More informationCanards at Folded Nodes
Canards at Folded Nodes John Guckenheimer and Radu Haiduc Mathematics Department, Ithaca, NY 14853 For Yulij Il yashenko with admiration and affection on the occasion of his 60th birthday March 18, 2003
More informationResearch Article Global Dynamics of a Competitive System of Rational Difference Equations in the Plane
Hindawi Publishing Corporation Advances in Difference Equations Volume 009 Article ID 1380 30 pages doi:101155/009/1380 Research Article Global Dynamics of a Competitive System of Rational Difference Equations
More informationRational Maps with Generalized Sierpinski Gasket Julia Sets
Rational Maps with Generalized Sierpinski Gasket Julia Sets Robert L. Devaney Department of Mathematics Boston University Boston, MA 02215, USA Mónica Moreno Rocha The Fields Institute University of Toronto
More informationAutonomous systems. Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous.
Autonomous equations Autonomous systems Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous. i f i(x 1, x 2,..., x n ) for i 1,..., n As you
More informationOn localized solutions of chains of oscillators with cubic nonlinearity
On localized solutions of chains of oscillators with cubic nonlinearity Francesco Romeo, Giuseppe Rega Dipartimento di Ingegneria Strutturale e Geotecnica, SAPIENZA Università di Roma, Italia E-mail: francesco.romeo@uniroma1.it,
More informationCournot Duopoly when the Competitors Operate Multiple Production Plants
ISSN 1974-4110 WP-EMS Working Papers Series in Economics, Mathematics and Statistics Cournot Duopoly when the Competitors Operate Multiple Production Plants Fabio Tramontana, (U. Ancona & U. Urbino) Laura
More informationSupersymmetric Quantum Mechanics and Geometry by Nicholas Mee of Trinity College, Cambridge. October 1989
Supersymmetric Quantum Mechanics and Geometry by Nicholas Mee of Trinity College Cambridge October 1989 Preface This dissertation is the result of my own individual effort except where reference is explicitly
More informationA System of Difference Equations with Solutions Associated to Fibonacci Numbers
International Journal of Difference Equations ISSN 0973-6069 Volume Number pp 6 77 06) http://campusmstedu/ijde A System of Difference Equations with Solutions Associated to Fibonacci Numbers Yacine Halim
More informationCONSTRUCTING SOLUTIONS TO THE ULTRA-DISCRETE PAINLEVE EQUATIONS
CONSTRUCTING SOLUTIONS TO THE ULTRA-DISCRETE PAINLEVE EQUATIONS D. Takahashi Department of Applied Mathematics and Informatics Ryukoku University Seta, Ohtsu 50-1, Japan T. Tokihiro Department of Mathematical
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationIntroduction to Applied Nonlinear Dynamical Systems and Chaos
Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium
More informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More informationCourse 212: Academic Year Section 1: Metric Spaces
Course 212: Academic Year 1991-2 Section 1: Metric Spaces D. R. Wilkins Contents 1 Metric Spaces 3 1.1 Distance Functions and Metric Spaces............. 3 1.2 Convergence and Continuity in Metric Spaces.........
More informationINTRODUCTION TO FRACTAL GEOMETRY
Every mathematical theory, however abstract, is inspired by some idea coming in our mind from the observation of nature, and has some application to our world, even if very unexpected ones and lying centuries
More informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 203 Outline
More informationORDINARY DIFFERENTIAL EQUATIONS
STEPHEN WIGGINS ORDINARY DIFFERENTIAL EQUATIONS SCHOOL OF MATHEMATICS UNIVERSITY OF BRISTOL Contents Preface 9 1 Getting Started: The Language of ODEs 13 2 Special Structure and Solutions of ODEs 23 3
More informationChapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics
Chapter 23 Predicting Chaos We have discussed methods for diagnosing chaos, but what about predicting the existence of chaos in a dynamical system. This is a much harder problem, and it seems that the
More informationSynchronization, Chaos, and the Dynamics of Coupled Oscillators. Supplemental 1. Winter Zachary Adams Undergraduate in Mathematics and Biology
Synchronization, Chaos, and the Dynamics of Coupled Oscillators Supplemental 1 Winter 2017 Zachary Adams Undergraduate in Mathematics and Biology Outline: The shift map is discussed, and a rigorous proof
More informationMonotone and oscillatory solutions of a rational difference equation containing quadratic terms
Monotone and oscillatory solutions of a rational difference equation containing quadratic terms M. Dehghan, C.M. Kent, R. Mazrooei-Sebdani N.L. Ortiz, H. Sedaghat * Department of Mathematics, Virginia
More informationINTRODUCTION TO FURSTENBERG S 2 3 CONJECTURE
INTRODUCTION TO FURSTENBERG S 2 3 CONJECTURE BEN CALL Abstract. In this paper, we introduce the rudiments of ergodic theory and entropy necessary to study Rudolph s partial solution to the 2 3 problem
More informationPeriodicity and Solution of Rational Recurrence Relation of Order Six
Applied Mathematics 3 79-733 http://d.doi.org/.436/am..377 Published Online July (http://www.scirp.org/journal/am) Periodicity and Solution of Rational Recurrence Relation of Order Si Tarek F. Ibrahim
More informationUniversity of Utrecht
University of Utrecht Bachelor thesis The bifurcation diagram of the second nontrivial normal form of an axially symmetric perturbation of the isotropic harmonic oscillator Author: T. Welker Supervisor:
More informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
More informationExample of a Blue Sky Catastrophe
PUB:[SXG.TEMP]TRANS2913EL.PS 16-OCT-2001 11:08:53.21 SXG Page: 99 (1) Amer. Math. Soc. Transl. (2) Vol. 200, 2000 Example of a Blue Sky Catastrophe Nikolaĭ Gavrilov and Andrey Shilnikov To the memory of
More informationA Two-dimensional Mapping with a Strange Attractor
Commun. math. Phys. 50, 69 77 (1976) Communications in Mathematical Physics by Springer-Verlag 1976 A Two-dimensional Mapping with a Strange Attractor M. Henon Observatoire de Nice, F-06300 Nice, France
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 2013 Outline
More informationSection 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d
Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d July 6, 2009 Today s Session Today s Session A Summary of This Session: Today s Session A Summary of This
More informationHamiltonian Chaos and the standard map
Hamiltonian Chaos and the standard map Outline: What happens for small perturbation? Questions of long time stability? Poincare section and twist maps. Area preserving mappings. Standard map as time sections
More information