Nicolas Grosjean

Thierry Huillet

SOME COMBINATORIAL ASPECTS OF DISCRETE NON-LINEAR POPULATION DYNAMICS

Keywords: non-linear population dynamics without and with immigration, Carleman transfer matrix, logistic map

Motivated by issues arising in population dynamics, we consider the problem of iterating a given analytic function a number of times. We use the celebrated technique known as Carleman linearization that turns (for a certain class of functions) this problem into simply taking the power of a real number. We expand this method, showing in particular that it can be used for population models with immigration, and we also apply it to the famous logistic map. We also are able to give a number of results for the invariant density of this map, some being related to the Carleman linearization.

Introduction

We consider simple 1-dimensional discrete-time dynamics: x n+1 = φ (x n ), x 0 = x, with the evolution mechanism φ (•) being an analytic function. Our main interest in these problems arises from population dynamics models describing the temporal evolution of some population with size x n ≥ 0. We first assume φ (0) = 0 (no immigration). Such non-linear models are amenable to a Carleman linearization giving x n from the initial condition x in terms of the n-th power of some uppertriangular infinite-dimensional transfer matrix which can be diagonalized. Equivalently, φ is h-conjugate to the linear map λx for some Carleman function h, would λ = φ (0) = {-1, 0, 1}. The coefficients of h, as a power series in x, are obtained from the left eigenvector of P with the eigenvalue λ. The Carleman linearization technique goes back the 60' ( [START_REF] Erdös | On analytic iteration[END_REF], [START_REF] Jabotinsky | Analytic iteration[END_REF], [START_REF] Kowalski | Nonlinear dynamical systems and Carleman linearization[END_REF], [START_REF] Gralewicz | Continuous time evolution from iterated maps and Carleman linearization[END_REF]). When λ = 1 (the critical case), we give the linear Carleman representation of x n , using a 'Jordanization' technique. Special such models arising in population dynamics are defined and investigated. We next consider the problem of computing the invariant density (and its support) of the dynamics in a chaotic population model regime, including quadratic maps. The study of the invariant measures of quadratic and related maps has a very long story starting in the 70' ( [START_REF] Bunimovich | A Transformation of the Circle[END_REF], [START_REF] Ruelle | Applications conservant une mesure absolument continue par rapport à dx sur [0, 1[END_REF], [START_REF] Bowen | Invariant measures for Markov maps of the interval[END_REF], [START_REF] Jakobson | Absolutely Continuous Invariant Measures for one-parameter Families of One-dimensional Maps[END_REF], [START_REF] Collet | Positive Liapunov exponents and absolute continuity for maps of the interval[END_REF]). We show that in some special cases, the h-conjugate representation of φ is useful for that purpose. We illustrate our point of view on the celebrated logistic population model φ (x) = rx (1 -x). Next we consider φ 0 (x) = c + φ (x), modelling some population dynamics with immigration c > 0. In the presence of a fixed point for φ 0 , such models are also Carleman linearizable; equivalently, φ 0 is shown to be g-conjugate now to an affine map for some explicit Carleman function g. As an illustration, we finally deal with the logistic population model with immigration. We develop its intimate relation to a family of companion logistic population models without immigration, the former being obtained from the latter through a suitable affine transformation. We exploit this deep connection to determine under which condition the logistic model with immigration is chaotic or not and, using this observation, we compute in some cases its invariant density.

The precise organization of the paper is as follows:

In Section 2, we recall and develop the Carleman transfer matrix linearization technique, including:

-its link with a conjugate representation of the map φ in the case φ (0) = 0 and φ (0) = {-1, 0, 1} .

-the application of this scheme to the specific critical case φ (0) = 1, leading to a method akin to the Jordanization of a matrix.

-the consequences of this construction in terms of the invariant measure of the dynamical system.

-the application of this general setting to a class of specific population evolution models.

In Section 3, using the above tools, we focus on the one-parameter logistic population model. Specifically, we characterize the loci and the types of the divergence of its invariant measure and we give a way to compute the disconnected components of its support for some parameter r range. For some values of the parameter r, we show how to compute explicitly the invariant density of the system.

In Section 4, expanding the tools introduced in Section 2 to include maps obeying φ 0 (0) = 0, we study the effect of adding immigration to population dynamics models, by relating it to a an affine conjugate equivalent of the new mechanism with immigration φ 0 (0) > 0. Once again, we apply these results to the logistic map. Our main result on this point is summarized in Figure 1 showing the values of the parameters (r, c) for which this topological conjugation is admissible. As a consequence, the chaoticity of the logistic model with immigration is revealed by the one of the corresponding model without immigration.

Carleman matrix in the triangular case α

0 = 0 With α k , k ≥ 1, real numbers, let φ (x) = k≥1 α k x k , α 1 = 0, be some smooth power series defined (convergent) in some neighborhood x - c < x < x + c , of the origin where -∞ ≤ x - c < 0 < x + c ≤ ∞ 1 .
We avoid the trivial linear case φ (x) = α 1 x. We shall let I c = (x - c , x + c ) be the interval of convergence. Note α 0 = 0. Consider the dynamical system (1)

x n+1 = φ (x n ) , x 0 = x.
Define the infinite-dimensional (Carleman) upper-triangular matrix

(2)

P (k, k ) = x k φ (x) k , k ≥ k ≥ 1.
By Faà di Bruno formula (see e.g. [START_REF] Comtet | Analyse combinatoire. Tomes 1 et 2[END_REF], Tome 1, p. 148), with B k,l (α 1 , α 2 , ...) the (ordinary) Bell polynomials in the coefficients

α k := x k φ (x) of φ (x), (3)
P (k, k ) = B k ,k (α 1 , α 2 , ..., α k -k+1 ) = k! * * c l k -k+1 l=1 α c l l c l ! ,
where the last double-star summation runs over the integers c l ≥ 0 such that

k -k+1 l=1 c l = k and k -k+1 l=1
lc l = k (there are p k,k terms in this sum, the number of partitions of k into k summands). In particular

P (k, k) = α k 1 and P (k, k + 1) = kα 2 α k-1 1 . P is called the Carleman 2 matrix of φ. If for example, φ (x) = x -x 2 , P (k, k ) = B k ,k (z, -1, 0, ...) | z=1
, the Hermite polynomials evaluated at z = 1. We conclude (see [START_REF] Jabotinsky | Analytic iteration[END_REF], [START_REF] Rabinovich | Solving nonlinear recursions[END_REF], [START_REF] Berkolaiko | Analysis of Carleman representation of analytical recursions[END_REF], [START_REF] Kowalski | Nonlinear dynamical systems and Carleman linearization[END_REF] and [START_REF] Gralewicz | Continuous time evolution from iterated maps and Carleman linearization[END_REF]):

Proposition: With e 1 = (1, 0, 0, ...) and x = x, x 2 , ... 3 , (4)

x n = e 1 P n x where P is an upper-triangular 'transfer' matrix with

P (k, k) = α k 1 =: λ k , k ≥ 1 (the eigenvalues of P ).
From (1), x n is also x n = φ n (x) where φ n is the n-th iterate of φ by composition and so ( 4) is an alternative linear representation of x n . Note

n≥0 λ n x n = e 1 (I -λP ) -1 x,
involving the resolvent of P .

Remark: We have 1 1 -uφ (x) = 1 + k≥1 u k φ (x) k = 1 + k≥1 u k k ≥k x k x k φ (x) k = 1 + k≥1 u k k ≥k x k P (k, k ) = 1 + k ≥1 x k k k=1 u k x k φ (x) k k k=1 u k P (k, k ) = x k 1 1 -uφ (x)
.

With φ n (x) = x n , the n-th iterate of φ and u := u, u 2 , ...

1 1 -uφ n (x) = 1+ k≥1 u k k ≥k x k P n (k, k ) = 1+ k ≥1 x k k k=1 u k P n (k, k ) = 1+u P n x.
Taking the derivative with respect to u at u = 0 gives φ n (x) = x n = e 1 P n x. Taking the k-th derivative with respect to u at u = 0 gives x k n = e k P n x. We conclude:

2 Carleman matrices are easily seen to be the transpose of Bell matrices. 3 Throughout, a boldface variable, say x, will represent a column-vector and its transpose, say

x , will be a row-vector.

Proposition:

If ψ (x) = k≥1 ψ k x k is some smooth observable, defining ψ = (ψ 1 , ψ 2 , ...), therefore (5) ψ (x n ) = k ψ k e k P n x = ψ P n x.
This generalizes [START_REF] Collet | Positive Liapunov exponents and absolute continuity for maps of the interval[END_REF].

By Cauchy formula, whenever φ is defined on the unit circle, we also have the Fourier representation

P (k, k ) = 1 2π
2π 0 e ik θ φ e -iθ k dθ.

Chaos for ( 1) is sometimes characterized by the positivity of its Lyapounov exponent defined by

λ (x) = lim N →∞ 1 N N -1 n=0 log φ (x n ) ,
for almost all x. Considering the sensitivity to the initial condition problem, we have

J n+1 := dx n+1 /dx = φ (x n ) dx n /dx = φ (x n ) J n . Therefore |J N | = N -1 n=0 φ (x n ) and λ N (x) := 1 N log |J N | → λ (x). Letting x = x, x 2
, ... , x := 1, x, x 2 , ... and D :=diag(1, 2, 3, ...) so that dx/dx = Dx, we observe

J N := dx N /dx = e 1 P N Dx.
v k (k ) = λ k -λ k -1 k -1 l=1 P (l, k ) v k (l)
gives the entries v k (k ), k > k ≥ 1, of v k by recurrence; and v k (k) can be left undeterminate. Developing, for k > k, we get

Proposition: Suppose λ = α 1 = 1, then v k (k ) = v k (k) k -k+1 j=2 * d1+..+dj-1=k -k j-1 l=1 P (d 1 + .. + d l-1 + k, d 1 + .. + d l + k) λ k -λ d1+..+d l +k , or v k (k ) = v k (k) k -k+1 j=2 k=k1<k2<..<kj-1<kj =k j-1 l=1 P (k l , k l+1 ) λ k -λ k l+1 ,
where the star-sum in the first identity (involving the integer partition of k -k) runs over the integers d l ≥ 1 summing to k -k.

Thus, with V = [v 1 , ..., v k , ...] , V P = D λ V where V is upper-triangular (so here invertible) and D λ =diag λ, λ 2 , ... , (6)

x n = e 1 V -1 D n λ V x.
This shows that a general non-linear dynamical system (1) generated by φ with φ (0) = α 0 = 0 and α 1 = {-1, 0, 1} is in fact a linear infinite-dimensional system with 'transfer matrix' P which can effectively be diagonalized.

Proposition: Defining h k (x) = v k x = k ≥k v k (k ) x k , this also means that for all k ≥ 1 x n = h -1 k λ nk h k (x) .
Proof: Indeed, (P x) l = φ (x) l and

v k P = λ k v k ⇒ v k P x = λ k v k x ⇒h k (φ (x)) = λ k h k (x) where h k (x) = v k x = k ≥k v k (k ) x k , for all k ≥ 1. Iterating, x 1 = φ (x) = h -1 k λ k h k (x) ⇒ x 2 = φ (x 1 ) = h -1 k λ 2k h k (x) ⇒ x n = φ (x n-1 ) = h -1 k λ nk h k (x) .
When k = 1, with v : = v 1 and ( 7)

h (x) = v x = k ≥1 v (k ) x k =: k ≥1 h k x k
and h -1 the analytic continuation of the inverse function of h, both obeying h (0) = h -1 (0) = 0, this in particular means (8)

x 1 := φ (x) = h -1 (λh (x)) and x n = h -1 (λ n h (x)) .
The importance of this formula relies on the fact that, in this way, x n can be evaluated directly, for any initial point x, without actually computing the intermediate values x 1 , ..., x n-1 . But this is at the expense of the computation of h and h -1 , which are simple special functions only in some exceptional situations.

Remarks:

(i)In a neighborhood of the origin, the analytic continuation of the inverse function and the inverse function itself defined on the range of h indeed coincide. By Lagrange inversion formula, h -1 , viewed as a power series, is locally defined in some neighborhood of the origin4 , with coefficients

g k := x k h -1 (x) = 1 k x k-1 (h (x) /x) -k .
By Faà di Bruno formula,

g 1 = h -1 (0) = 1/h 1 and for k ≥ 2 (9) g k = h -k 1 k-1 l=1 (-1) l (k + l -1) k-1 B k-1,l h 1 , h 2 , ... ,
where (k) l = k (k -1) ... (k -l + 1), h k = h k+1 /h 1 and B k,l the 'ordinary' Bell polynomials.

(ii)If 0 < λ < 1 (a subcritical case), (8) says x n → 0 whatever x as n → ∞. Furthermore, a first order Taylor development shows that

x n = h -1 (0)+ h -1 (0) λ n h (x)
with h -1 (0) = 0 and h -1 (0) = 1/h 1 . So x n goes to 0 geometrically fast at rate λ and λ -n x n → h (x) /h 1 .

Proposition:

If we choose v k (k) = v 1 (1) k , then h k (x) = h 1 (x) k =: h (x) k . Proof: Consider the ratio r (x) = h k (x) /h (x) k . It holds that r (φ (x)) = r (x) because r (φ (x)) = h k (φ (x)) /h (φ (x)) k = λ k h k (x) λ k h 1 (x) k = r (x) .
But then all the k-derivatives r (k) (0), k ≥ 1, of r at x = 0 vanish because φ (0) = 0 and φ (0

) = λ = 1. Therefore r is constant on x ∈ I c . Choosing v k (k) = v 1 (1) k , then r (0) = lim x→0 h k (x) /h (x) k = v k (k) /v 1 (1) k = 1, so r (x) = 1.
One possible way to achieve this is to assume without loss of gener-

ality that v 1 (1) = h 1 = 1 = v k (k) . Corollary: If v k (k) = v 1 (1)
k , the upper-triangular matrices V and V -1 in ( 6) are the Carleman transfer matrices of the power series h and h -1 in (8); similarly, the matrix D λ is the transfer matrix of the linear power series x → λx.

Proof: by definition V (k, k ) = v k (k ) = x k h k (x) = x k h (x)
k and by La-

grange inversion formula, V -1 (k, k ) = x k h -1 (x) k .
2.2. The critical case λ = α 1 = 1 and Jordanization of P . Let us consider the critical case, α 0 = 0 and α 1 = 1. The transfer matrix P has only 1s on its diagonal, and is subsequently not (infinite-)diagonalizable. Let us search instead for a family of column vectors w 1 , w 2 , . . . , such that

P w 1 = w 1
and, if k > 1, P w k = w k-1 + w k , which corresponds to the Jordanization of the infinite matrix P . A solution is given by w

1 = (1, 0, . . . ), w k (k ) = 0 if k > k, w k (1) = 0 if 5 k > 1 and otherwise w k (k ) = 1 P (k -1, k ) w k-1 (k -1) - k i=k +1 P (k -1, i)w k (i)
With W the matrix whose columns are the w k , I the infinite identity matrix and S the infinite (shift) matrix with 1s on its superdiagonal, we then have

P W = W (I + S)
and the powers of P satisfy

P n W = W n k=0 n k S k
Note that since P (k, k + 1) = kα 2 in the critical case, we have that

W (k, k) = α -(k-1) 2 (k -1)! ,
and because W (k, k) = λ k for some λ, W cannot be the transfer matrix of some power series. By construction, the first row of W is (1, 0, . . . ).

Define σ(k, k ) = 0 if k ≤ k together with σ(k, k ) = 1=i1<i2<•••<i k+1 =k k l=1 P (i l , i l+1 ) if k > k ≥ 1 and, if k = 0, σ(k, k ) = δ 1,k . For example, σ(2, 4) = P (1, 2)P (2, 4)+P (1, 3)P (3, 4), σ(2, 3) = P (1, 2)P (2, 3
). Then, we conjecture that

W -1 (k, k ) = σ(k -1, k ), k ≥ k ≥ 1.
We get:

Proposition: Since we have

x n = e 1 P n x, observing e 1 W = e 1 W -1 = e 1 and e 1 S k = e k+1 , (10)
x n = n k=0 n k e k+1 W -1 x = x + n k=1 n k k ≥k+1 σ(k, k )x k .
This constitutes the infinite-dimensional matrix representation of x n in the critical case.

Slow extinction in the critical case. With φ (0) = 0 and φ (0) = α 1 = 1, consider the critical dynamics,

x n+1 (x) = φ (x n (x)) , x 0 (x) = x.
Suppose x n (x) approaches 0 as n gets large. An order-two Taylor development of φ near x = 0 therefore gives

x n+1 (x) ≈ φ (0) x n (x) + 1 2 φ (0) x n (x) 2 = x n (x) + 1 2 φ (0) x n (x) 2 ,
leading (assuming α 2 := φ (0) < 0) to [START_REF] Jakobson | Absolutely Continuous Invariant Measures for one-parameter Families of One-dimensional Maps[END_REF] x n (x) ∼ x/ (1 -(nα 2 ) /2) , as n is large, indeed going to 0. We conclude:

Proposition: For a critical dynamics (1) for which

α 1 = 1, if α 2 := φ (0) < 0, nx n (x) → 2x/ (-α 2 )
, and x n vanishes at slow algebraic rate n -1 .

2.3. Special models arising in population dynamics. Because we are interested in population dynamics systems for which x n ≥ 0 is the size of some population at time n, we shall limit ourselves to dynamical systems of type (1) generated by φ with φ (0) = 0 and α 1 = φ (0) > 0, α 1 = 1 (recall however that the latter construction leading to ( 6) and ( 8) holds even if α 1 < 0). The initial condition x will be assumed to belong to the domain [0, x b ] where

x b = inf (x > 0 : φ (x) = 0) ≤ x + c , possibly with x b = ∞ if x + c = ∞.
We also need to assume that the maximal value φ * that the function φ can take on [0,

x b ] is ≤ x b so that φ maps I = [0, x b ] onto J ⊆ [0, x b ]
. If these conditions hold, the dynamics (1) will be said a population dynamics model.

Examples of population models:

(i) (logistic map) φ (x) = rx (b -x) with 4 ≥ α 1 = br > 0, α 2 = -r < 0, x b = b > 0, x + c = ∞. Here I = [0, b] and J = 0, rb 2 /4 . (ii) (homographic map) φ (x) = rx/ (1 + ax) with α 1 = r > 0, α 2 = -ar < 0, x b = x + c = ∞. Here I = [0, ∞] and J = [0, r/a] . (iii) (Ricker map) φ (x) = rx exp (-ax) with α 1 = r > 0, α 2 = -ar < 0, x b = x + c = ∞. Here I = [0, ∞] and J = [0, φ * = 1/a] ⊂ I. (iv) φ (x) = x * -a (x * -x) 2 + b with a, b, x * > 0, ax 2 * + b = x 2 * . Here I = 0, 2 b/ (1 -a) and J = 0, b/ (1 -a) 1 - √ 1 -a . Note φ (0) = α 1 = a < 1. For this map, φ (x) = q -1 (aq (x) + b) where q (x) = (x -x * ) 2 .
Note that models (i) and (iii) are single-humped maps, meaning a typically smooth and non-negative function with exactly one critical point (where φ vanishes), a maximum, and at most one point of inflection to the right of the maximum.

2.4. Invariant measure of specific maps. Consider a map φ from some interval I onto J ⊆ I. Then the invariant measure µ solves the Perron-Frobenius equation µ = φ -1 • µ, [START_REF] Ding | Statistical Properties of Deterministic Systems[END_REF] and [START_REF] Lasota | Fractals, and Noise. Stochastic Aspects of Dynamics[END_REF]. When φ is chaotic, a density solution f possibly exists, thereby solving f (x) = I δ (x -φ (y)) f (y) dy, [START_REF] Thunberg | Periodicity versus chaos in one-dimensional dynamics[END_REF]. This is also [START_REF] Kowalski | Nonlinear dynamical systems and Carleman linearization[END_REF]

f (x) = y:φ(y)=x f (y) φ (y) = K(x) k=1 f φ -1 k (x) φ -1 k (x) = K(x) k=1 f φ -1 k (x) φ φ -1 k (x)
, where K (x) is the number of φ-antecedents of x belonging to the set {x : f (x) > 0}.

Consider now a non-monotone map φ from some interval I onto itself (a surjection whose domain and range coincide with I). For almost all initial condition x (drawn at random uniformly on [0, 1]), for almost all φ (0) in a chaotic range of φ, a density solution f of ( 12) can exist and be made explicit. For instance, consider the case investigated by [START_REF] Luévano | The Schröder functional equation and its relation to the invariant measures of chaotic maps[END_REF] where in addition: ∀x ∈ I, φ -1 (x) is made of K antecedents (independently of x).

Given φ and κ ∈ R, consider then the Schröder functional equation, with unknown function s, s (φ (x)) = κs (x) .

Taking the derivative and then the modulus,

φ (x) |s (φ (x))| = |κ| |s (x)|. Thus for each branch k, φ φ -1 k (x) |s (x)| = |κ| s φ -1 k (x)
, showing from [START_REF] Kowalski | Nonlinear dynamical systems and Carleman linearization[END_REF], see [START_REF] Luévano | The Schröder functional equation and its relation to the invariant measures of chaotic maps[END_REF],

that if |κ| = K f (x) = |s (x)| .
This shows that for chaotic maps φ mapping some interval I into itself and for which ∀x ∈ I, φ -1 (x) is made of K antecedents (independent of x), if it has an invariant density, it can be explicitly obtained from the solution of a Schröder functional equation. But s can also be related to h, the Carleman function of φ.

Corollary: Let φ map I onto itself with K (x) = K and λ = φ (0) > 0. Recalling the representation (8) of φ: h (φ (x)) = λh (x), with α = log K λ, we get

s (x) = h (x) 1/α and so f (x) = 1 |α| h (x) 1/α-1 h (x) .
This shows that in some cases, (8) can also be useful for the computation of the invariant density of φ. Let us illustrate these facts for the logistic map:

-With I = [0, 1], consider the logistic map φ (x) = rx (1 -x) with b = 1 and r = 2. For this map, K = 2 but φ maps [0, 1] to [0, 1/2] ⊂ [0, 1]. Then, for every x ∈ (0, 1/2), φ (x) = h -1 (2h (x)) where h (x) = -log (1 -2x) /2 and h -1 (x) = 1 -e -2x
/2, the inverse of h. The invariant measure is a Dirac at the fixed point x = 1/2 of φ. If 2 < r < r 3.5699456 (the Feigenbaum constant), the invariant measure is a Dirac measure concentrating at the fixed points obtained in the period doubling process of the logistic map.

-Consider the logistic map φ (x) = rx (1 -x) with b = 1 and r = 4: φ maps [0, 1] into [0, 1] and K = 2 for all x. Then φ

(x) = h -1 (4h (x)) where h (x) = arcsin ( √ x) 2 and h -1 (x) = sin ( √ x) 2 6 . The invariant measure is µ (dx) = f (x) dx where f (x) = π -1 (x (1 -x)) -1/2 = h 1/2 (x) = |s (x)| where s (x) = arcsin ( √ x) solves s (φ (x)) = 2s (x) .
We now wish to further investigate some aspects of the invariant density problem for the general logistic map.

Invariant densities for φ(x) = rx(1 -x)

Let us further consider the logistic map φ (x) = rx (1 -x) with b = 1 and assume now 4 > r = λ > 2. This map is not into I but onto only. In this case the dynamics moves to the interval [φ (r/4) , r/4] in finite time and stays there 7 : the support of the invariant measure is expected to lie within this interval and, whenever 4 ≥ r > 3.5699456 (the Feigenbaum constant for the onset of chaos), the set of points with positive measure possibly consists in a finite union of disconnected subintervals embedded within this support. Within the interval r ∈ [3.5699456..., 4], it is widely known that there are subintervals where the map φ is not chaotic (has negative Lyapounov exponent) corresponding to x n reaching a cycle of any length. In these subintervals, the invariant measure is expected to be a Dirac measure.

The purpose of this section is to exhibit some features of the invariant density in the chaotic regimes, when this density exists. In the following, we will use the notation α n (r) := φ n (1/2).

To illustrate this section, we will often refer to a specific value of r, r 3.6785735, which is the one studied by Ruelle in [START_REF] Ruelle | Applications conservant une mesure absolument continue par rapport à dx sur [0, 1[END_REF]. This value is the smallest r for which there is a density with full support [1/r, r/4], and it solves (r -2) r 2 -4 = 16. 6 In this case, the inverse function of h (x) = arcsin √ x 2 is sin √ x 2 , restricted to the interval [0, π/2] whereas its analytic continuation h -1 is the same function now defined on the whole positive real line (its maximal domain of convergence). 7 x > r/4 has no antecedent and cannot be reached, the function is strictly increasing on ]0, φ(r/4)], and the interval [φ(r/4), r/4] is mapped onto itself.

3.1.

Divergence points of the invariant density. When r = 4, the invariant density diverges likes x -1/2 near 0 and (1 -x) -1/2 near 1, {0, 1} being the boundary of the support of the invariant density. When r = 4, φ maps 1/2 to 1 and then to 0, a fixed point of the dynamics: α 0 (r) = 1/2, α 1 (r) = 1 and α 2 (r) = 0. The density diverges at {0, 1}, the successive image of 1/2 under φ.

When r

3.6785735, φ maps 1/2 to r/4 and then to 1/r and finally to 1 -1/r, a fixed point of the dynamics: α 0 (r) = 1/2, α 1 (r) = r/4, α 2 (r) = 1/r and α 3 (r) = 1 -1/r. The density is expected to diverge at points {1/r, 1 -1/r, r/4}, the successive images of 1/2 under φ.

This seems to be a generic fact.

For a generic value of r, when x = α 1 (r) = r/4, there is only one antecedent, namely y = 1/2 and there, φ (y = 1/2) = 0, so, in view of ( 12), the invariant density f (x) of the mass at x = r/4 must diverge, if it exists. This is also true for x = φ 1 (r/4) = φ 2 (1/2) = α 2 (r): for this x indeed, y = r/4 is an antecedent and because f (y) diverges for this value of y, so does f (x). By induction, Proposition: For all x = α n (r), n ≥ 1, there is a divergence of the invariant density.

For the Ruelle value of r, the sequence φ n (1/2) = α n (r), n ≥ 1, stops (converges) after three steps, so one expects a divergence of the density f in that point which is 1 -1/r (the fixed point of φ) but also at the intermediate points α 1 (r) = r/4, α 2 (r) = 1/r. Following this argument, when r 3.89087 (defined by α 5 (r) = α 3 (r) = α 4 (r)), one expects a divergence of the density at the four points α 1 (r) = r/4 = φ 1 (1/2), α 2 (r) = φ 2 (1/2), α 3 (r) = φ 3 (1/2) and α 4 (r) = φ 4 (1/2). The sequence stops there because for this value of r, α 5 (r) = α 3 (r), α 6 (r) = α 4 (r) and no new α n (r) is expected, this sequence oscillating between α 3 (r) and α 4 (r) with α 3 (r) > α 4 (r). This concerns special values of r for which f has full support and finitely many peaks. This poses the following general question: let C m be a cycle of length m of φ. Consider the set

R := {r in the chaotic regime: ∃n, m < ∞: α n (r) ∈ C m and α n-1 (r) / ∈ C m } .
If r ∈ R and if corresponding map φ admits a density f , then there must be n + m -1 peaks where f diverges, by the arguments above. If r / ∈ R and if f exists, then f has infinitely many peaks within its support. We don't know if R is dense within the subset of values of r in the chaotic regime, nor do we have a specific value of r / ∈ R. For an r / ∈ R, the sequence {α n (r)} does not enter a cycle in finite time and one expects a density diverging on this countable set values for x within [φ (r/4) , r/4].

It remains to determine the type of the divergence of f . We shall develop our arguments for the Ruelle value of r. We have

φ (y) = φ 1 2 + y - 1 2 = r 4 -r y - 1 2 2 , φ (y) = r (1 -2y) .
Put x = r/4 -ε and y = 1/2 -ε with , > 0. Because φ (1/2) = 0, we have φ (y) = 2rε and φ (y) = r/4 -rε 2 = x = r/4 -ε. Thus rε 2 = ε.

In view of

f (x) = y:φ(y)=x f (y) φ (y) , therefore f (r/4 -ε) ∼ f (1/2) 2 √ rε .
Similarly, with some linear relation between , and , the antecedents of 1/r + are r/4 -and 1 -r/4 + , the latter being outside of the support of f so its value there is zero. The antecedents of 1 -1/r + are 1/r + and 1 -1/r -. The antecedents of 1 -1/r -are 1 -1/r + and 1/r -, the latter being outside of the support of f . With these relations, we find that

(1/r + ) ∼ a/ √ , f (1 -1/r -) ∼ b/ √ , f (1 -1/r + ) ∼ c/ √ with a := f (1/2) φ (1/r) φ (1/r) 1/2 , b := a φ (1/r) φ (1/r) -1 , c := a φ (1/r) 3/2 φ (1/r) -1 .
Proposition: There is an algebraic divergence of f at the peaks of order -1/2 whereby f is integrable.

3.2.

Disconnected invariant measure support. As a polynomial in the variable r, α n (r) has degree 2 n -1. When r is larger than the Ruelle value r 3.6785735, whenever the regime is chaotic, the support of the invariant measure is made of one single piece, namely the full interval [α 2 (r) = φ (r/4) , α 1 (r) = r/4].

The above value of r is when α 3 (r) is a fixed point of φ, namely α 3 (r) = φ (α 3 (r)) = α 4 (r). Because α 3 (r) = 1 -1/r and φ (1/r) = 1 -1/r, this Ruelle value of r is also obtained when α 2 (r) = r 2 (1 -r/4) /4 = 1/r, indeed leading to r 3.6785735. When r becomes slightly less than the Ruelle value, the support of the invariant density of φ splits in the two pieces [α 2 (r) , α 4 (r)] ∪ [α 3 (r) , α 1 (r)], because α 4 (r) < α 3 (r). These two pieces each split again in two additional pieces for the value of r for which α 5 (r) = α 7 (r) and α 6 (r) = α 8 (r), corresponding to α 5 (r) and α 6 (r) being respectively one of the two known fixed points of φ 2 (which are not the fixed points {0, 1 -1/r} of φ 1 = φ). Slightly below this value of r, the support of the invariant density of φ splits into the four pieces

[α 2 (r) , α 8 (r)] ∪ [α 6 (r) , α 4 (r)] ∪ [α 3 (r) , α 7 (r)] ∪ [α 5 (r)
, α 1 (r)]: this transition is seen to occur at r 3.5925722. The next third step generates the 2 3 = 8 pieces

[α 2 (r) , α 16 (r)] ∪ [α 12 (r) , α 8 (r)] ∪ [α 6 (r) , α 14 (r)] ∪ [α 10 (r) , α 4 (r)] ∪ [α 3 (r) , α 15 (r)] ∪ [α 11 (r) , α 7 (r)] ∪ [α 5 (r) , α 13 (r)] ∪ [α 9 (r) , α 1 (r)]
whenever α 16 (r) = α 12 (r), α 14 (r) = α 10 (r), α 15 (r) = α 11 (r) and α 13 (r) = α 9 (r) corresponding to α 12 (r) , α 10 (r), α 11 (r) and α 9 (r) being the four fixed points of φ 4 which are neither the ones of φ 1 nor the ones of φ 2 . For some special value of r which can be computed, the four pieces split into these eight pieces.

This binary splitting process of the support can be iterated until one reaches the critical Feigenbaum value of r, namely r 3.5699456 where the invariant measure is expected to be singular (whose support consists in an uncountable number of points). For values of r less than this Feigenbaum value and larger than three, the invariant measure is a Dirac measure concentrated on the period two cyclic points appearing in the period doubling process onsetting after r = 3.

3.3.

Invariant measure for some specific values of r. Let us consider the Ruelle value of r. For this value of r, the quartic map

ψ (x) = (r -2) 2 x (1 -x) [1 + (r -2) x (1 -x)] maps [0, 1] into [0, 1]
and it has K (x) = 2 branches for each x ∈ [0, 1] \ {1/2}. Furthermore, ψ (x) = q 2 (x) := q • q (x) where q (x) = (r -2) x (x -1) and, with a (x) = (r -1 -rx) / (r -2)

φ (x) = a -1 (q (a (x))) and φ 2 (x) := φ • φ (x) = a -1 (ψ (a (x))) .
As stated before, this value is also one for which α 3 (r) = α 4 (r) (α 3 (r) is a fixed point of φ).

Let now h be the Carleman function associated to ψ, defined by

ψ (x) = h -1 (λh (x)), λ = ψ (0) = (r -2)
2 > 2, as from ( 8). The dynamical system generated by ψ has an absolutely continuous invariant measure, given by f ψ (x) = |s (x)| where s (ψ (x)) = 2s (x) and s (x) = h 1/α (x), α = log 2 λ > 1. Therefore so does the dynamical system generated by φ 2 with invariant density given by f

φ 2 (x) = |a (x)| f ψ (a (x)), with support [1/r, 1 -1/r] .
And then so does the dynamical system generated by φ itself, with the density

f φ (x) = 1 2 f φ 2 (x) + φ -1 • f φ 2 (x) , where φ -1 • f φ 2 (x) = φ (x) f φ 2 (φ (x)) is the image density of f φ 2 under φ -1 . The support of f φ is [1/r, r/4] = [1/r, 1 -1/r] ∪ φ -1 [1/r, 1 -1/r].
We conclude that there is a Carleman inspired analytic expression of the invariant density for the logistic model when r takes on the Ruelle value. This is not the only one case.

Looking for an affine function a (x) = c -dx such that φ 4 (x) = a -1 (Ψ (a (x))) for some map Ψ = ψ • ψ with ψ (x) = ax (1 -x) [1 + bx (1 -x)], such that Ψ maps [0, 1] onto itself, leads to another solution to the above Ruelle scenario, namely: a -1.7214305, b 0.1592945, c -2.2648957, d -5.5297915 and r 3.5925722. The latter is the value of r where the support of the invariant measure splits from 2 pieces to 4. In that case, we have λ := Ψ (0) = a 2 , Ψ having two branches on [0, 1] \ {1/2}. Using the same Carleman matrix approach as previously stated, it is possible to compute the invariant measure associated to the mapping Ψ (now with h the Carleman function of Ψ and s (x) = h 1/α (x) for α = log 2 a 2 ). From that result and following the arguments developed for the Ruelle value, one can clearly compute the invariant density first of the function φ 4 on a given subset of its support, namely f φ 4 , and then get the invariant density of φ on its support by

f φ = 1 4 f φ 4 + φ -1 • f φ 4 + φ -2 • f φ 4 + φ -3 • f φ 4 .
Proposition: For the particular value of r 3.5925722, the invariant density of φ can be computed.

Although for r 3.6785735 and for r 3.5925722, the maps φ were not from [0, 1] onto [0, 1], the powers φ 2 and φ 4 of these maps were shown to be affine-conjugate to maps ψ and Ψ which now map [0, 1] onto [0, 1], thereby amenable to the Schröder and Carleman representation of the invariant densities.

4. The general case α 0 = 0 Let φ 0 (x) = k≥0 α k x k = α 0 + φ (x), now with α 0 = φ 0 (0) = 0. Consider now the dynamical system (13)

x n+1 = φ 0 (x n ) , x 0 = x. Let (14) P 0 (k, k ) = x k φ 0 (x) k , k, k ≥ 0.
Note P 0 no longer is triangular. Specifically, for k, k ≥ 1 Then, with e 1 = (0, 1, 0, 0, ...) and x = 1, x, x 2 , ... ,

P 0 (k, k ) = x k (α 0 + φ (x)) k = k l=0 k l α k-l 0 x k φ (x) l = k l=1 k l α k-l 0 P (l, k ) = k l=1 k l α k-l 0 B k ,l (α 1 , α 2 , ..., α k -l+1 ) , together with P 0 (0, k ) = δ 0,k , P 0 (k, 0) = α k 0 , for k ≥ 0. Hence P 0 = BP where B (k, l) = k l α k-l 0 , k ≥ l ≥ 0 and P =       1 0 • • • • • • 0 . . . P . . .       , so with P 0 =      1 0 • • • • • • α 0 α 2 0 BP . . .      and B (k, l) = k l α k-l 0 , k ≥ l ≥ 1,
x n = e 1 P n 0 x. This shows that a general non-linear dynamical system (13) generated by φ 0 with φ 0 (0) = 0 is also a linear infinite-dimensional system with complete 'transfer matrix' P 0 . Because P 0 no longer is triangular, its eigenvalues are not known nor is its diagonalization easy, if even possible. So P n 0 is very complex. Note in particular that x n (x = 0) = e 1 P n 0 e 0 = P n 0 (1, 0), the (1, 0) -entry of P n 0 . In some cases however, the general problem (13) can be taken back to a simpler problem of type [START_REF] Berkolaiko | Analysis of Carleman representation of analytical recursions[END_REF], see [START_REF] Gralewicz | Continuous time evolution from iterated maps and Carleman linearization[END_REF]. Let y n = x n -ρ for some real number ρ. Then

y n+1 = φ 0 (y n + ρ) -ρ, y 0 = x -ρ.
Suppose there is a real number ρ such that φ 0 (ρ) = ρ (the map φ 0 has a fixed point). Then ( 16)

y n+1 = φ (y n ) := φ 0 (y n + ρ) -ρ, y 0 = x -ρ,
where φ (0) = φ 0 (ρ) -ρ = 0. If this is the case, the new system (16) now generated by φ is of the form (1). It can be solved as a (1) model with the new

φ (y) = k≥0 α k (y + ρ) k -ρ =: l≥1 y l α l , where l!α l = k≥l ρ k-l (k-l)! (k!α k ) is of convolution type. And then, assuming λ := α 1 = 1, (17)
x

n = h -1 λ n h (x -ρ) + ρ,
where h is associated to the new generator φ (just like h in (8) was to φ in (1)).

Introducing h (x) = h (x -ρ), this is also

x n = h -1 λ n h (x) similar to (8), except
that h (0) = 0. Note also that φ (0) =: α 1 = φ 0 (ρ) is not necessarily > 0. And depending on φ 0 (ρ) < 1 (or > 1), ρ is a stable (unstable) fixed point of [START_REF] Lasota | Fractals, and Noise. Stochastic Aspects of Dynamics[END_REF].

Equivalently, (see e.g. [START_REF] Gralewicz | Continuous time evolution from iterated maps and Carleman linearization[END_REF])

Proposition: With Q the upper-triangular Carleman matrix of φ (easily diagonalizable with V Q = D λ V where D λ =diag 1, λ, λ 2 , ... ), the Carleman matrix P 0 of φ 0 therefore obeys (18)

P 0 = B -1 ρ QB ρ = (V B ρ ) -1 D λ V B ρ ,
where B ρ is the lower-triangular Carleman matrix of the shift function b ρ (x) = x-ρ: B ρ (k, l) = (-1)

k-l k l ρ k-l and B -1 ρ (k, l) = k l ρ k-l . Remarks:

-The only Carleman matrices which are lower-triangular are the ones associated to an affine map as the one above.

-With φ 0 (ρ) = 1, φ 0 (ρ) , φ 0 (ρ) 2 , ... and ρ = 1, ρ, ρ 2 , ... , for all fixed point ρ of φ 0 : P 0 ρ = φ 0 (ρ) = ρ showing that ρ is a right eigenvector of P 0 associated to its eigenvalue 1.

-If the map φ 0 has more than one real fixed point, the latter construction holds for any of these fixed points, showing that ( 17) is not unique.

-If the map φ 0 has no real fixed point, P 0 is not real-diagonalizable.

For example, the matrix P 0 associated to the map φ 0 (x) = x + α 0 / (x + 1) with φ 0 (0) = α 0 = 0 is not diagonalizable. However, this model has a fixed point at infinity. Exchanging ∞ and 0 can be done while using the transformation y = 1/x. The dynamics for the ys is thus y n+1 = φ (y n ), y 0 = y = 1/x, with φ (y) = 1/φ 0 (1/y) = y (1 + y) 1 + y + α 0 y 2 , now with a fixed point at y = 0 and of type [START_REF] Berkolaiko | Analysis of Carleman representation of analytical recursions[END_REF], with φ rational. We have φ (0) = α 1 = 1 and φ (0) = α 2 = 0, a critical model, therefore amenable to Jordanization.

Another example is φ 0 (x) = x+x 2 +α 0 with φ 0 (0) = α 0 = 0. This map has no real fixed point but it has two complex fixed points ρ = ±i √ c. The matrix P 0 associated to this map φ 0 is not real-diagonalizable but it is complex-diagonalizable. In such cases, only [START_REF] Rabinovich | Solving nonlinear recursions[END_REF] holds, but not [START_REF] Sagitov | A special family of Galton-Watson processes with explosions[END_REF] where h is real-valued. 4.1. An equivalent conjugation representation of φ 0 having a fixed point. Assuming λ := α 1 = φ 0 (ρ) = 1, with h (x) = h (x -ρ), we obtained

x 1 = φ 0 (x) = h -1 λ h (x -ρ) + ρ = h -1 λ h (x) ,
where h is associated to the generator φ. Define g (x) := h (x) -h (0) = h (x -ρ)h (-ρ), now obeying g (0) = 0 and g (ρ) = -h (-ρ) . Clearly then, with µ

:= λ -1 h (-ρ) = 1 -λ g (ρ) , (19)
x 1 = φ 0 (x) = g -1 λg (x) + µ ,
showing that φ 0 is g-conjugate to the affine function λx + µ, with c = φ 0 (0) = g -1 (µ) > 0. This can be iterated to give (20)

x n = g -1 λ n g (x) + µ n , with λ n = λ n and µ n = µ 1 + λ + ... + λ n-1 = g (ρ) 1 -λ n .
If λ < 1 (ρ is a stable fixed point of φ 0 ), x n → ρ.

4.2.

Population models with immigration. In population dynamics systems for which x n ≥ 0, we shall limit ourselves to dynamical systems of type [START_REF] Lasota | Fractals, and Noise. Stochastic Aspects of Dynamics[END_REF] generated by φ 0 with φ 0 (0) = c > 0. No need to require here anymore that φ 0 (0) = φ (0) = α 1 > 0. The initial condition x will be assumed to belong to the domain [0, x b ] where x b = inf (x > 0 : φ 0 (x) = 0), possibly with x b = ∞. Here φ 0 (0) > 0 interprets as an immigration rate. We also need to assume that the maximal value φ * 0 that φ 0 (x) can take on [0, x b ] is ≤ x b so that φ 0 maps I = [0, x b ] onto J ⊆ I. These φ 0 are amenable to the formalism [START_REF] Rabinovich | Solving nonlinear recursions[END_REF] and [START_REF] Thunberg | Periodicity versus chaos in one-dimensional dynamics[END_REF].

Examples of population models with immigration φ 0 (0) = c > 0:

(i) (logistic map) φ 0 (x) = c + rx (b -x), r > 0, x b = br + √ ∆ / (2r) > b > 0, ∆ = (br) 2 + 4rc, x + c = ∞. Here I = [0, x b ] and J = 0, c + rb 2 /4 ⊆ I provided φ * 0 = c + rb 2 /4 ≤ x b . (i ) (logistic map') φ 0 (x) = c + rx (b -x), r < 0, b > 0. In this setup, φ 0 (0) = φ (0) = α 1 = rb < 0 -If φ 0 (b/2) = c + r (b/2) 2 > 0, then x b = x + c = ∞ and I = J = [0, ∞] . -If φ 0 (b/2) = c + r (b/2) 2 ≤ 0, then x b < ∞. Here I = [0, x b ] and J = [0, c] ⊆ I provided c ≤ x b . (ii) (homographic map) φ 0 (x) = c + rx/ (1 + ax) with x b = ∞ = x +
c . Here I = [0, ∞] and J = [c, c + r/a] ⊂ I. For this map, g in ( 19) is readily seen to be an homographic function itself, using the matrix representation of an homography, matrix product translating into composition of homographic maps. In these 4 examples, there is a smallest ρ > 0 for which φ 0 (ρ) = ρ.

(v) Consider the quadratic population model with immigration φ 0 (x) = x + x 2 + c with φ 0 (0) = c > 0. In this example, φ 0 (x) has no real fixed point and clearly with

x n+1 = φ 0 (x n ), x 0 = x, x n ∼ λ (x, c) 2 n
, with λ (x, c) > 1 depending on both x and c. x n drifts to ∞ at doubly exponential speed with n.

(vi) Consider the quadratic population model with immigration φ 0 (x) = x + c/ (x + 1) with φ 0 (0) = c > 0. In this example, φ 0 (x) has no real fixed point either and clearly with

x n+1 = φ 0 (x n ), x 0 = x, x n ∼ √ 2cn → ∞, whatever x > 0.
x n drifts algebraically slowly to ∞.

4.3.

Logistic population models with or without immigration. Consider the logistic population dynamical system without immigration (21)

x n+1 = rx n (1 -x n ) , x 0 = x, where r ∈ (0, 4] . Let y n = ax n + b, with b = (1 -a) /2. Then, if a = r/R (22) y n+1 = c + Ry n (1 -y n ) , y 0 = y, where (23) c = (r -1) 2 -(R -1) 2 4R .
For the dynamics y n to be a logistic population dynamical system with immigration, mapping some interval I onto J ⊆ I, elementary algebra shows that the pair (r, R) must lie in one of the shaded zone Z of Figure 1 below, where |R| ≤ 4 and r ∈ (0, 4].

The detailed explanation of Figure 1 is as follows: for any value of R and r, there exists a c such that the dynamics rx(1 -x) and Ry(1 -y) + c are related. For the latter to be an admissible population dynamics including immigration, it has to satisfy several conditions. First, the number c has to be positive (which is related to the lines R = r and R = 2 -r on the graph). With this condition, the dynamics in y has to belong to one of three classes: R > 0 (corresponding to the vertically dashed area), R < 0 and the dynamics never reaches zero (the diagonally dashed area), or R < 0 and the dynamics reaches zero (the horizontally dashed area). The domain of the first-kind dynamics obviously corresponds to the superior half-plane, while the domains for the second and the third one are separated by R = r(2-r)/2. Another condition to be an admissible dynamics is that it has to map one finite (because the rx(1 -x) is defined on [0, 1]) domain into itself. This condition is related to the lines R = r(r -2)(r 2 -2r -4)/8 for the first-kind dynamics, R = -r for the second-kind one and r 2 -R 2 -2r -2R -4 = 0 for the third-kind one. Finally, to be consistent with the rx(1 -x) dynamics, the initial condition in y has to be mapped in an initial condition in x that belongs to [0, 1], which is in relation to the line r = R.

This means that, starting from a dynamical system of type (22) for some R chosen in the range |R| ≤ 4, for any value of r intersecting the shaded zone with the horizontal line of equation y = R, there is a positive immigration rate c given by ( 22) such that the dynamics of the ys is mapped into the dynamics (21) using the reverse affine transformation x n = (y n -b) /a. We conclude:

Proposition: With |R| ≤ 4, if r ∈ (0, 4] is in the shaded zone Z of the graph and also in the chaotic range for (21), then the logistic population dynamical system (22) with immigration rate c given by ( 23) is chaotic.

Starting from a logistic model without immigration and with reproduction rate R > 0 (not) in the chaotic region, adding immigration c can lead to a (c, R) model ( 22) either (chaotic) or non chaotic, depending on the chosen value of r: immigration can either (destabilize) stabilize a (non) chaotic system, [START_REF] Stone | Effects of immigration on the dynamics of simple population models[END_REF]. We also note that when R is negative, but not too negative, there is no r large enough to lie in the chaotic region of (21). Example: Let R = -4 and r = 4 (in Z), so that c = 1, corresponding to the logistic population dynamical system with immigration (24) y n+1 = 1 -4y n (1 -y n ) , y 0 = y.

Here a = -1 and b = 1. With h (x) = arcsin ( √ x) 2 , we have

y n = 1 -h -1 (4 n h ((1 -y))) ,
and this dynamics is chaotic.

It can be checked that the representation (20) takes the alternative form

y n = sin 2 (-2) n arcsin √ x + π 6 (1 -(-2) n ) ,
corresponding to g (x) = arcsin ( √ x), λ = -2 and ρ = 1/4 (with g (ρ) = π/6), in the notations of subsection 4.1.

(21) with r = 4 has the invariant density f (x) = π -1 (x (1 -x))

-1/2 and the logistic population dynamical system with immigration (24) has the invariant density f (y) = f ((y -1) / (-1)) = π -1 (y (1 -y))

-1/2 , so identical to f (with Lyapounov exponent log 2).

Conclusion

Using Carleman linearization techniques applied to 1-dimensional discrete-time population dynamical systems, we gave a technique to compute the current population state x n , for any initial point x, without actually computing the intermediate values x 1 , ..., x n-1 . This technique was shown to be related to the characterization of the invariant density measure, when it exists. But this is at the expense of the computation of h and h -1 , which are "simple" special functions only in some exceptional situations, such as for specific parameter values of the logistic map with or without immigration. What "simple" means and the class of models for which these functions are "simple" and/or lead to chaotic behavior are largely open problems. This methodology has also recently proved useful in the context of discrete-time branching process for which the map φ is absolutely monotone: for the family of so-called generalized linear-fractional branching processes first introduced in [START_REF] Sagitov | A special family of Galton-Watson processes with explosions[END_REF] and further studied in [START_REF] Grosjean | Additional aspects of the generalized linear-fractional branching process[END_REF], the function h has a simple structure and the iteration of φ does not lead of course to chaos being one-to-one on the unit interval.

2. 1 .

 1 The case |λ| = 1. Suppose λ := α 1 = ±1 and let v k P = λ k v k , define the left-row-eigenvector v k of P associated to the eigenvalue λ k := λ k . Then, with k > k ≥ 1,

  the incomplete lower-triangular Pascal (binomial) matrix. Note Bx k = (α 0 + x) k : the full Pascal matrix B is the Carleman matrix of the shift function b (x) = α 0 + x. It also holds B (k, k) = 1, k ≥ 0 and P 0 admits a LU factorization.

(

  iii) (Ricker map) φ 0 (x) = c + rx exp (-ax) with x b = x + c = ∞. Here I = [0, ∞] and J = [c, c + 1/a] ⊂ I.

(

  iv) φ 0 (x) = c+x * -a (x * -x) 2 + b. Here I = [0, x b ] and J = 0, c + x * -√ b ⊆ I .

Corollary: With |R| ≤ 4 ,

 4 suppose r ∈ (0, 4] is in the shaded zone Z of the graph.If x n = h -1 (r n h (x))is the conjugacy representation (8) of x n , then, with a = r/R and b = (1 -a) /2 y n = ah -1 (r n h ((y -b) /a)) + b is a conjugacy solution y n of (22), started at y. Corollary: If (22) is chaotic and corresponding (21) has an invariant density f , then (22) has the image invariant density: f (y) = f ((y -b) /a) /a.

Figure 1 .

 1 Figure 1. Map of the admissible correspondences between R and r. A point in the dashed areas corresponds to an admissible transformation. The miniatures of graphs indicate the kind of dynamics for Ry(1 -y) + c.

If h has a positive convergence radius, so does h -1 .

Acknowledgments:

T. Huillet acknowledges partial support from the "Chaire Modélisation mathématique et biodiversité". N. Grosjean and T. Huillet also acknowledge support from the labex MME-DII Center of Excellence (Modèles mathématiques et économiques de la dynamique, de l'incertitude et des interactions, ANR-11-LABX-0023-01 project). The authors would like to thank Geneviève Rollet for some fruitful discussions.