Discretized Fast-Slow Systems near Pitchfork Singularities

Discretized Fast-Slow Systems near Pitchfork Singularities Luca Arcidiacono, Maximilian Engel, Christian Kuehn arxiv:90.065v [math.ds] 8 Feb 09 February 9, 09 Abstract Motivated by the normal form of a fast-slow ordinary differential equation exhibiting a pitchfork singularity we consider the discrete-time dynamical system that is obtained by an application of the explicit Euler method. Tracking trajectories in the vicinity of the singularity we show, how the slow manifold extends beyond the singular point and give an estimate on the contraction rate of a transition mapping. The proof relies on the blow-up method suitably adapted to the discrete setting where a key technical contribution are precise estimates for a cubic map in the central rescaling chart. Keywords: pitchfork bifurcation, slow manifolds, invariant manifolds, loss of normal hyperbolicity, blow-up method, discretization. Mathematics Subject Classification 00): 3E5, 37M99, 37G0, 3C5, 39A99. Acknowledgments This work was supported by the German Research Foundation DFG) via the SFB-TRR 09. Luca Arcidiacono acknowledges support from the graduate program TopMath of the Elite Network of Bavaria and the TopMath Graduate Center of TUM Graduate School at Technische Universität München. Introduction We study the dynamical system generated by the two-dimensional cubic polynomial map ) ) x x x + h xy x P : = ) ) + λε).) y ȳ y + hε for small ε, h > 0. The parameter ε implies a time scale separation between the fast Zentrum Mathematik, Fakultät für Mathematik, Technische Universität München, Boltzmannstr. 3, D-8578 Garching bei München

variable x and the slow variable y, while the parameter h is viewed as the stepsize in an explicit Euler method for the ordinary differential equation ODE) dx = ẋ = xy x ) + λε, dτ dy = ẏ = ε, dτ.) which is the normal form of a fast-slow system exhibiting a pitchfork singularity at the origin x, y) = 0, 0). Indeed, for ε = 0, we can view y as a bifurcation parameter for the flow in the x-variable: when y < 0, equation.) has a hyperbolic sink at 0, y) and, when y > 0, there are three equilibria with 0, y) unstable and the other two, y, y) and y, y), locally asymptotically stable. The origin x, y) = 0, 0) is called singular since hyperbolicity is lost at this point, and this is also the case for the map.). We will analyze the dynamics close to the origin for small ε, h > 0. Since we focus on the local behaviour around the singularity, we will neglect potential higher order nonlinearities in the majority of this work, but will show how to adapt the proof when including those. For our analysis we will make use of the blow-up method [, ], which has turned out to be a successful tool for treating singular points of fast-slow systems. It was first applied to fastslow systems by Dumortier and Roussarie [3] to gain insight in the dynamics around nonhyperbolic equilibria. The method uses a non-injective transformation that maps a higher dimensional object such as a sphere onto the non-hyperbolic equilibrium constituting the singularity. The dynamics on this larger, blown up version of the singularity may then be desingularized by an appropriate rescaling of time and exhibit partially) hyperbolic behaviour. Then one can use dynamical systems techniques to analyze the dynamics in blown-up space. Finally, a typical result allows one to extend invariant manifolds past the singular point in the blown down system; see e.g. [3, Chapter 7] for an introduction and [6, 5, 7, 7, 0,, ] for an, of course non-exhaustive, list of different applications to planar fast-slow systems. By means of the blow-up method Krupa and Szmolyan [0, ] analyze different kinds of singularities in fast-slow ODEs, i.e. fold, canard, transcritical and pitchfork singularities, and show how certain invariant manifolds, so-called slow manifolds, extend around the singular points for small ε > 0. In the case of fold points, Nipp and Stoffer [8] transform the blow-up technique to the corresponding explicit Runge-Kutta, in particular Euler, discretization and prove the extension of slow manifolds for the discrete time system around the singularity. Whereas they apply an abstract existence theory for invariant manifolds developed in [8], Engel and Kuehn [] use direct estimates in the blow-up charts to prove the extension of slow manifolds for transcritical singularities. In both cases, a crucial aspect of the discretized blow-up lies in finding the right rescaling of the step size h. In a similar spirit, we investigate, how trajectories of.) behave near the origin and show how the slow manifold may be continued beyond the pitchfork singularity in the discrete setting. We prove that, depending on the sign of λ, trajectories starting in the vicinity of the single slow branch near {0, y) R : y < 0} are attracted exponentially to one of the parabolic branches near {± y, y) R : y > 0}. Furthermore, for λ = 0, we show that canard-type orbits can track the unstable branch {0, y) R : y > 0}. Our analysis uses three charts that cover different parts of the blown up space around the singularity. We track trajectories through several checkpoints along a curve of fixed points of the cubic map and give estimates on the contraction of the transition mappings.

In this way, we also give an alternative way of proof to the result in [] for the ODE case, by letting the step size h 0 in.). This paper is structured as follows. After giving a short introduction to continuous time fast-slow systems and pitchfork singularities in Section, we formulate the setup and main results of this paper at the start of Section 3. The major part of Section 3 is dedicated to the proof of the main theorem which is divided into several steps. We start with λ 0. The relevant coordinate changes are discussed in Section 3.. Sections 3.3 and 3.5 describe the dynamics in the vicinity of the branches of the critical manifold, which allows to define the slow manifolds and control contractivity of the transition map. In Section 3. we describe the continuation of a slow manifold through the blown-up singularity by direct estimates on the trajectories. Finally Section 3.6 combines the results in a blown down version, which finishes the proof for λ 0. The required modifications to cover the canard case λ = 0 are outlined in Section 3.7. Section 3.8 shows how the previous results can be adapted to a more general setting. Continuous-time fast-slow systems with pitchfork singularity. Fast-slow ODEs Fast slow systems occur in various fields of science such as neurobiology or chemistry and are usually found as a system of ODEs with two time scales, this means, they are of the form dx = x = fx, y, ε), dt dy = y = εgx, y, ε), x R m, y R n, 0 < ε, dt.) where f, g, are C k -functions with k 3. The small parameter ε consitutes the separation between two time scales. The variables x and y are often called the fast variables) and the slow variables) respectively. The time variable in.), denoted by t, is termed the fast time scale. By a change of variables, one can also consider the slow time scale τ = εt and rewrite.) as ε dx dy = εẋ = fx, y, ε), dτ = ẏ = gx, y, ε). dτ.) The singular limit ε = 0 can be seen from two different perspectives corresponding with the two time scales. Setting ε = 0 in.) yields x = fx, y, 0), y = 0,.3) which is called the layer problem or fast subsystem), since we can view the equation layer-wise parametrized by the constant y. Setting ε = 0 in.) gives the differential algebraic equations 0 = fx, y, 0), y = gx, y, 0),.) 3

called the reduced problem or slow subsystem). The flow of.), the so-called slow flow, is restricted to the set S 0 := {x, y) R m R n fx, y, 0) = 0}, which consists of equilibria of the layer problem.3). We refer to this set as the critical set or often also critical manifold, in case it is a manifold. A subset S S 0 is called normally hyperbolic if the matrix D x fx, y, 0) R m m has no eigenvalue with vanishing real part for all x, y) S. In the vicinity of normally hyperbolic submanifolds of S 0, the dynamics can be very well described for sufficiently small ε > 0: Fenichel Theory [6, 9, 3, 9] gives the existence of a locally invariant manifold, the slow manifold S ε, which lies close to S 0 and maintains the stability properties of the layer problem.3). Furthermore, the restriction of.) to S ε is a regular perturbation of the reduced problem.). However, points p R m R n, which do not satisfy normal hyperbolicity are called singularities in this context and are more delicate to handle. From the view point of the layer equation.3) singularities often correspond to bifurcations of the fast subsystem, and the breakdown of normal hyperbolicity is typically associated with the intersection of multiple parts of S 0 at the point p where degeneracy of D x fp) follows from the Implicit Function Theorem. In the case of a pitchfork singularity, which we will consider in the following, we are precisely in such a situation.. Pitchfork singularity in continuous time We consider a two-dimensional fast-slow system of the form.) where the critical manifold resembles the shape of a pitchfork, and we call the associated non-hyperbolic singularity a pitchfork singularity. Such a situation occurs when the vector field f satisfies f0, 0, 0) = 0, x f0, 0, 0) = 0, y f0, 0, 0) = 0, xx f = 0, xxx f0, 0, 0) 0, xy f0, 0, 0) 0. These conditions guarantee that, for ε = 0, there is a non-hyperbolic equilibrium at the origin where the critical manifold has a transversal self-intersection and one part of the branches crosses the other tangentially to the x-direction. In particular, we assume that xxx f0, 0, 0) < 0, xy f0, 0, 0) > 0, such that the singularity is supercritical. Furthermore, we assume g0, 0, 0) > 0 such that the slow dynamics pass through the origin in positive y-direction. In other words, we consider the problem of how the slow dynamics behave in the vicinity of a splitting into three critical branches see [, Figure ]). The case of a subcritical pitchfork singularity or the situation of g0, 0, 0) < 0 are less challenging, since the dynamics only heads into the direction of one critical branch, and will therefore not be treated in this paper. There is a linear change of coordinates see []) which brings the system into the normal form ẋ = xy x ) + λε + h x, y, ε), ẏ = ε + h x, y, ε)),.5)

where h and h satisfy h x, y, ε) = O x y, xy, εx, εy, ε ), h x, y, ε) = O x, y, ε). Since we are mainly interested in the local dynamics around the origin, we may initially ignore the higher order terms and only consider the system ẋ = xy x ) + λε, ẏ = ε..6) The critical manifold is given as S 0 = {x, y) R : y = x } {x, y) R : x = 0}. For negative y there is only one stable equilibrium at x = 0, while for positive y we have an unstable equilibrium at x = 0 and two locally asymptotically stable ones at x = ± y. 3 Pitchfork singularity in discrete time A time-discretization of equation.6) by the explicit Euler method with time step size h > 0 yields the map ) ) x x x + h xy x P : = ) ) + λε). 3.) y ȳ y + hε As in continuous time, the system induced by 3.) clearly possesses the critical manifold S 0 = {x, y) R : y = x } {x, y) R : x = 0}, consisting of fixed points of 3.) for ε = 0. We split the set S 0 into the four branches S 0 a = {x = 0, y < 0}, S + a = {y = x, x > 0}, S a = {y = x, x < 0}, S 0 r = {x = 0, y > 0}. 3.) By linearization we see that these four branches are normally hyperbolic as long as for x, y) S 0 \ {0, 0)} we have + hy. 3.3) Since we want to restrict the analysis to the non-hyperbolic singularity at the origin x, y) = 0, 0), we always assume that h is chosen small enough so that 3.3) holds as well as the same stability properties as in the time-continuous case. For example, for a fixed initial condition with y 0 < 0, we have to ensure + hy 0 > which yields the restriction h < / y 0, which then implies that Sa 0 is normally hyperbolic and locally attracting. Note that we shall still use the notation in 3.) in this context. In contrast to the continuous case the fixed points on Sa 0 are not globally stable, but only inside an interval around zero of size + y. Outside this interval solutions diverge in growing oscillations. Compare h also with the reasoning of Lemma 3.5 in Section 3.3. Due to normal hyperbolicity and according to [8, Theorem.], for ε sufficiently small, there exist corresponding invariant slow manifolds Sa,ε,h 0 S a,ε,h + S a,ε,h S0 r,ε,h. However, the four branches of the critical manifold S 0 intersect at the origin x, y) = 0, 0), where we have D x P 0, 0) =, i.e. we observe the loss of normal hyperbolicity as in the ODE case. 5

3. Main result We want to investigate where points around S 0 a,ε,h get mapped to by iterations of P in order to find the continuation of S 0 a,ε,h beyond the singularity. For that purpose, fix some ρ > 0, let J R be a small interval containing 0 and define the section in around the point 0, ρ ) on the critical branch S 0 a by in = {x, y) R y = ρ, x J}. In particular, we always assume that the initial condition is chosen on in and J is sufficiently small so that trajectories effectively start close to the attracting slow manifold Sa,ε,h. 0 We are going to follow trajectories of.) starting in in up to height y = ρ. Since the line {y = ρ } can only be reached in case ρ N, we introduce ρ as the closest reachable εh height, which then satisfies ρ ρ < εh. This allows us to define the sections 0 out = {x, y) R y = ρ, x J}, ± out = {x, y) R y = ρ, x ρ J} around the points ρ, ρ ), ρ, ρ ) or 0, ρ ) on the branches S a +, Sa or Sr 0 respectively. Depending on the sign of λ we will show that the transition maps Π ± : in ± out or Π 0 : in 0 out if λ = 0 are well defined. The transition maps are given by the N-fold of P where N = [ ρ ρ ] is the closest integer to. For the discrete setting, induced by the εh εh map 3.), we have the following main result see Figure for an illustration of the case λ > 0). Theorem 3.. Fix ρ > 0 and let λ R. Then there are ε 0, h 0 > 0, depending on λ, such that for all ε 0, ε 0 ] and all h 0, h 0 ] the following holds. T) If λ > 0, the set in including the point in Sa,ε,h) 0 is mapped by Π + into a subset ) of + out that contains the point + out S a,ε,h + and has a width of order O c h) K εh for some constants c, K > 0. T) If λ < 0, an analogous statement holds with out, S a,ε,h and Π instead of + out, S + a,ε,h and Π+. T3) If λ = 0 the slow manifolds Sa,ε,h 0 and Sr,ε,h 0 coincide with the critical branches Sa 0 and Sr 0 and are connected by a canard solution, i.e. Π 0 in Sa,ε,h) 0 Sr,ε,h. 0 The set in ) ρ gets mapped by Π 0 into 0 out and its image has a width less than h ρ ) εh J. Remark 3.. An analogous result for the continuous-time system.6) has been shown in []. The sections are chosen to be ± out = {x, y) R x = ρ, y ρ J} such that the transition maps Π ± : in ± out are well defined if signλ) = ±. The image Π ± in ) ± out is an interval about the corresponding point on the slow manifold S ± a,ε and has a size of O e C ε ) for some constant C > 0. Note that we will assume h 0 ε 0 ρ to obtain meaningful time lengths for the dynamical analysis. For further details on the choice of h 0 see Section 3.6, but notice that it immediately implies the stability restriction h < ρ from the discussion below 3.3). 6

y S + a,ε,h + out x in S 0 a,ε,h Figure : Extension of the slow manifold Sa,ε,h 0 bold line) through the neighbourhood of the origin for λ > 0. Further the critical manifold S 0 dashed line) and a sample trajectory from in to out + dotted line). 3. Transformation to the chart coordinates The proof of Theorem 3. uses the blow-up method for the dynamical system induced by the map P : R R x x x + h xy x ) + λε) y ε ȳ ε = y + hε ε, 3.) h h h where the fast-slow separation parameter ε and the stepsize h are also seen as variables. The quasi-homogeneous blow-up transformation around the pitchfork singularity is given as Φ : B := S 0, ) 0, ) R, x = r x, y = r ỹ, ε = r ε, h = r h. x, ỹ, ε, h, r) x, y, ε, h), The transformation of the x, y, ε)-coordinates is the same as in the continuous-time case see []). The change of variables in h is chosen such that the map is desingularized in the relevant charts. We exclude 0 from the domain of h since at h = 0 every point is a neutral fixed point. Due to the transformation h = h/r we have to exclude 0 from the domain of r as well. The transformation Φ induces a map P := Φ P Φ on the manifold B. We analyse the dynamics of P by using the charts K i, i =,, 3, K : D := R R + R + 0 R + R R R + 0 R + x, r, ε, h ) x, y, ε, h), K : D := R R R + R + R R R + R + x, y, r, h ) x, y, ε, h), K 3 : D 3 := R R + R + 0 R + R R + R + 0 R + x 3, r 3, ε 3, h 3 ) x, y, ε, h), 7

which are given by K : x = r x y = r ε = rε h = r h, K : x = r x y = ry ε = r h = r h, K 3 : x = r 3 x 3 y = r3 ε = r3ε 3 h = r3 h 3. To switch between different chart coordinates we use the following coordinate changes and κ : x = x ε y = ε r = r ε h = h ε, κ : x = x y ) r = r y ) ε = y h = h y, κ 3 : x 3 = x y ) r 3 = r y ) ε 3 = y h 3 = h y, κ 3 : x = x 3 ε 3 y = ε 3 r = r 3 ε 3 h = h 3 ε 3. For the proof of Theorem 3. we will proceed as follows. Transforming 3.) using the coordinate changes K i i =,, 3) induces a dynamical system on D i R in chart coordinates. Trajectories of 3.) are analyzed via their corresponding transformed versions in each of the charts. In every chart K i we will define sets Σ in i and Σ out i, show that the transition maps Π i : Σ in i Σ out i are well defined and study their contractivity. The mappings Π ± are then built by connecting the three chart-wise transition maps, which are elaborated in Section 3.6 in more detail. We refer to K as the entering chart, as we start our analysis in this chart and trajectories are brought closer to the origin. Charts of the type of K are often called rescaling charts, since the transformation is basically a rescaling with suitable powers of the fast-slowseparation constant ε. In this chart, the dynamics arbitrarily close to the origin are analyzed. Finally the exiting chart K 3 is used to describe, how trajectories exit the vicinity of the origin and is crucial for the contractivity statement of Theorem 3.. 3.3 Dynamics in the entering chart Fix some ρ > 0 and also consider the case λ > 0. from now on until Section 3.7. The case λ < 0 can be treated analogously, see also Section 3.6 for more details. Further take ε, h > 0 sufficiently small. During the next sections we will specify what sufficiently small means for ε and h such that ε 0 and h 0 are determined. In the coordinates x, r, ε, h ) of the first chart K, the set in is given as Σ in = { ) x, ρ, δ, ν ρ : x J }, for which we define δ := ε ρ and ν := hρ. 8

We investigate the dynamics within the domain D := { x, r, ε, h ) R : r [0, ρ], ε [0, 6δ], h [0, ν] }. In order to find an expression for the map 3.) in terms of the entering chart K, we first rewrite ȳ = y + εh in K -coordinates as This yields r = r + r h r ε. r = r h ε ). 3.5) The remaining three equations of 3.) in K -coordinates read as r x = r x + r h r x r r x ) + λr ε ), r ε = r ε, r h = r h. Hence, by using 3.5), we obtain the maps x = h ε ) [x + h x x ) + λr ε )], r = h ε ) r, ε = h ε ) ε, h = h ε ) h. 3.6) The dynamics in r, ε and h can be calculated explicitly for the first chart. Lemma 3.3. For ξ 0 := h 0)ε 0) > 0, the trajectories of 3.6) in r, ε, h are given by ξ0 n ξ0 r n) = r 0), ε n) = ε 0) ξ 0 ξ 0 n), h n) = h 0) ξ 0 n, ξ 0 for n N, n < ξ 0. Proof. We define η := h ε. Multiplying the last two equations of 3.6) gives the relation η = h ε = h ε ) h ε = η) η. Given the initial condition η0) = ξ 0 > 0, we obtain the solution ηn) = ξ 0 n). Hence, we can compute n h n) = h 0) k=0 ηk)) = h 0) n k=0 Similarly, we obtain the formulas for r and ε. ηk) ηk + ) = h 0) η0) ηn) = h 0) ξ 0 n. ξ 0 9

Furthermore, we can observe from equations 3.6) that the set {ε = 0} is invariant for the dynamics and, for given r, h, consists of the two-parameter family of invariant one-dimensional lines {ε = 0, r = r, h = h }. Each of these lines has a fixed point located at x, r, ε, h ) = 0, r, 0, h ), which has a three-dimensional center eigenspace and a one-dimensional stable eigenspace in x - direction with eigenvalue h recall that h hρ < < ). In other words, the 3 two-dimensional plane S 0 a, = {x, r, ε, h ) D : x = 0, ε = 0} is an invariant manifold in D only consisting of fixed points, attracting in the x -direction and neutral in the other directions, and corresponding to the branch Sa 0 of the critical manifold. In particular, for each h 0 we have the fixed point We obtain the following statement: p 0 ah ) = 0, 0, 0, h ). Proposition 3.. The invariant manifold S 0 a, extends to a center-stable invariant manifold M 0 a, at p 0 a0)) which is given in D by a graph x = l r, ε, h ) for some mapping l. Furthermore, ε is increasing in D and thereby, in particular, on M 0 a,), whereas h, r are decreasing in D and thereby also on M 0 a,). Proof. This is an immediate consequence of the considerations above and classical center manifold theory cf. e.g. [8, Chapter 5A]). From Lemma 3.3, we can see immediately that as long as r n), h n) > 0, we have n < ξ 0. Hence, the claim follows from the formulas in Lemma 3.3. Note that, on {r > 0, ε > 0, h > 0}, the manifold Ma, 0 corresponds to the union of the slow manifolds Sa,ε,h. 0 Assume we iterate system 3.6) until r n) reaches a value less or equal to ρ. This means that the trajectory of 3.) has reached a point above the line y = ρ / in x, y-coordinates.) As already remarked before, when defining ± out, the level ρ might not be hit precisely. However, let us assume for simplicity that we are in the situation of reaching r = ρ after n N iterates. Also in the next sections we will assume in a similar way that specific values are hit since the small errors do clearly not change our results.) This means we have ξ 0 n ξ 0 =, and reversely, n = 3ξ 0 = 3. νδ Denote by Π : Σ in R 3 the transition map after iterations. We can deduce the νδ following Lemma: Lemma 3.5. For δ > 0 sufficiently small, we have { Π Σ in ) Σ out := x [, ], r = ρ, ε = 6δ, h = ν }, and M 0 a, Σ in as well as M 0 a, Σ out are non-empty sets. 0

Proof. From the explicit solutions in Lemma 3.3, we obtain { Π Σ in ) r = ρ, ε = 6δ, h = ν }. For ε = 0 the x -equation of 3.6) reads Consequently we have x = x + h x x ) ) = h)x hx 3. 3.7) x < h )x for x > 0 and x > h )x for x < 0. Additionally a direct computation yields that for x we have h x h )x This means that x lies in the cone bordered by ± h )x for all x that satisfy x h see also Figure ). Thus contraction towards the fixed point at 0 is guaranteed for all initial values of x in that domain and the speed is at least the linear rate h. x h )x x h h h )x Figure : For x [ h, h ] the image x = x +h x x )) lies inside the cone spanned by ± h )x. Note that h n) stays inside the interval [ ν, ν] for n n and the stable eigenvalue in x -direction is given by h along the manifold S 0 a,. We may choose ν small enough so that ρ J [ ν, ν ]. Hence it follows from standard perturbation arguments for some constant c < ν. Therefore we have Π Σ in ) Σ out for a sufficiently small choice of δ. Since for δ 0 the point Ma, 0 {ε = δ} approaches 0, ρ, 0, ν), we can deduce that Ma, 0 {ε = δ} lies in Σ in for δ small enough. Additionally, since the manifold Ma, 0 is invariant under the forward iterations of 3.6), we have that Ma, 0 Σ out is non-empty as well. that, for δ sufficiently small, the map Π is a contraction with rate c) 3 νδ

3. Dynamics in the rescaling chart We use κ to transfer the set Σ out to the second chart and define { [ Σ in := κ Σ out ) = x δ, δ ], y = } δ, r = ε = ρδ, h = νδ. For the transformation of 3.) via the chart K, first observe that since ε = ε we have r = r. 3.8) Using the coordinates of K, the remaining equations from 3.) become r x = r x + r h [ r x r y r x ) + λr ], r ȳ = ry + r h r, r h = r h, which can be simplified with 3.8) so that we get in total x = x + h [ x y x ) + λr ], ȳ = y + h, r = r, h = h. Since r and h stay constant in this chart, we can plug in the values from Σ in to write the maps as ) x = x + νδ x y x ) + λρδ =: fx, y ), 3.9) ȳ = y + νδ. We denote this two-dimensional map by F : R R, F x, y ) = fx, y ), y + νδ ). Furthermore, we abbreviate λ := λρδ, such that we have fx, y ) = x +νδ x y x ) + λ ). We will see that the appropriate exiting set in chart K is given by { [ ] Σ out := x minλρ, )δ, δ + µ, y = } δ, r = ε = ρδ, h = νδ. Recall from chart K the attracting center manifold Ma, 0 see Proposition 3.). Similarly to [], this manifold corresponds with the global manifold Ma 0 on the blow-up manifold ) B. In chart K we therefore have the attracting center manifold Ma, 0 = κ M 0 a, with, reversely, M 0 a, = κ M 0 a, ) ), whose behaviour is described in the following main result of this section recall that we are still in the setting λ > 0).

Theorem 3.6. The set I := [ δ, gets mapped by iteration of F into the set [ I := δ minλρ, )δ, ] δ { δ } R + µ ] { δ } R, where µ = λρδ. In particular, the transition mapping defined and maps Ma, 0 Σ in to Ma, 0 Σ out. Π : Σ in Σ out is well For the proof of this theorem we will analyze the evolution of our starting set at different heights see Figure 3), treated in the propositions below. Before we do that, we collect some properties of the function f in the following Lemmas 3.7 3.9. y I I 3 x I I Figure 3: Path from I to I Theorem 3.6) via I to I Proposition 3.0), I to I 3 Proposition 3.) and I 3 to I Proposition 3.). Firstly, we characterize the positive fixed points of f, y ) on y -fibres. Lemma 3.7. For any y R, the mapping x fx, y ) has precisely one positive fixed point x y ), satisying the equation y = x y ) λ x y ). 3.0) As long as y [ δ, δ ], the family of fixed points x y ) is monotonically increasing and satisfies x y ) < δ +. Proof. Fixed points of f, y ) are characterized by the equation x y x ) + λ = 0. 3

y x y ) x Figure : The curves of fixed points of f, y ) in the x, y )-plane. As x = 0 is not a solution, we can solve this for y and obtain y = x λ x. The graph of x x λ x consists of two branches. Notice that x λ x > 0 for x < 0. For x > 0 the graph is monotonically increasing since the summands x and λ x are). Thus we have a curve of positive fixed points, that we can also parametrize by y, call them x y ). It remains to verify the bound of x y ). Using λ <, we have ) δ + λ ) δ + > δ + = Consequently the point ) δ +, δ lies below the graph of x x λ x, so right of the curve {x y ), y )} of fixed points. Since the curve is increasing, all other fixed points x y ) with y δ also satisfy the estimate x y ) δ +. Secondly, we find f, y ) to be monotonically increasing on a suitable interval and, by using Lemma 3.7, we find invariant sets under f, y ). Lemma 3.8. Let y [ δ, δ ]. Then the following holds:. The function f, y ) is monotonically increasing on [ δ, δ + ] for all ν sufficiently small.. The set [0, δ + ] is positively) invariant under f, y ). Furthermore, when y 0, the set [ δ, δ ] is positively) invariant under f, y ). Proof. We compute the derivative x fx, y ) = + νδ y 3νδ x, δ.

We have already shown that f, y) is increasing on [ δ, δ + ]. To prove the claimed invariance of [0, δ + ] and [ δ, δ ], it thus suffices to consider only the above that f δ +, y) < δ +. Since additionally f0, y) > 0, the claimed positive invariance of [0, δ + ] follows. Regarding the invariance of [ δ, δ ] in case of y 0, first note that the unique x y) x 0) = δ 3 λρ) < δ f δ, y) > δ and f δ, y) < δ, yielding the invariance of the interval [ δ, δ ]. x which shows that the cubic function f, y ) may have two stationary points located at x = ± 3 ν δ + y ) and then is monotonically increasing in between these. By choosing ν sufficiently small, we can achieve that the stationary points exist for all y [ δ, δ ] and that 3 ν δ + y ) δ +. The choice of ν can be made independently of δ, for all arbitrarily small δ. E. g. for δ the inequality is fulfilled for all ν. Hence, the first claim follows. 5 Since we have precisely one positive fixed point by Lemma 3.7, the graph of the continuous function f, y ) hits the diagonal precisely once. As f0, y ) = νδ 3 λρ > 0 and lim x fx, y ) = we can conclude that it crosses the diagonal at x and f, y ) lies above the diagonal i.e. fx, y ) > x ) for x [0, x y )) and below the diagonal i.e. fx, y ) < x ) for x x y ), ). If y 0 we can widen the estimate fx, y ) < x to x, x y )), because there is overall only one fixed point. Figure 5: The graph of f, y) together with the diagonal x x and the fixed point x y). endpoints of the intervals. Recall from Lemma 3.7 that x y) δ + and, hence, we can conclude with the fixed point x 0) is given by x 3 0) = λ = δ 3 λρ), which can be directly checked with the defining equation 3.0). As the curve of fixed points is increasing, we can estimate for sufficiently small δ and all y 0. With the considerations above, we thus have 5

Finally, we can show the following contraction property of f on relevant intervals. Lemma 3.9. Let y [ δ, 8 δ ]. The function f, y ) restricted to the interval [ δ, δ ] is a contraction with constant ν 8. Proof. Recall from the proof of Lemma 3.8 that the derivative reads x fx, y ) = + νδ y 3νδ x. As stated in Lemma 3.8, the map f, y ) is monotonically increasing on [ δ, δ ] for y [ δ, 8 δ ], and, hence, the derivative x fx, y ) is non-negative for all x, y ) E := [ δ, δ ] [ δ, 8 δ ]. Thus we obtain max x,y ) E x fx, y ) = max x,y ) E x fx, y ) = max x,y ) E + νδ y 3νδ x 8 νδ δ = ν 8. Hence, the claim follows by a standard application of the mean-value theorem. We now turn to showing the transitions from I i to I i+ for i =,, 3 see Figure 3). Each transition is formulated in one of the following Propositions. Proposition 3.0. The set I = [ ] δ, δ { δ } gets mapped by iterations of 3.9) into the set I := [ δ 3 λρ δ), 8δ 3 λρ + δ) ] { 8 δ }, where δ) := δ e 6 δ. Proof. It takes 8 δ νδ ) = 8 ν δ iterations to get from y = δ to y = 8 δ with steps of size νδ. Since [ δ, δ ] is invariant under f, y ) for all y [ δ, 8 δ ] by Lemma 3.8 and f, y ) is a contraction with constant ν by Lemma 8 ) 3.9, the image of I under iterations of 3.9) has a width of at most δ ν 8 ν δ. 8 Since the exponential limit lim ν 0 ) ν 8 ν δ = e 6 δ 8 is attained monotonically from below, we can bound this width via δ ) ν 8 ν δ 8 δ e 6 δ =: δ). Hence, any trajectory starting in I goes through the set [x,min δ), x,max + δ)] { } 8 δ, where x,min and x,max denote the minimal and the maximal value, respectively, that the fixed point x y ) attains while y varies over [ δ, 8 δ ]. Since, according to Lemma 3.7, the curve of fixed points is increasing, we have x,min = x δ ) and x,max = x 8 δ ). 6

It remains to verify that x δ ) δ 3 λρ and x 8 δ ) 8δ 3 λρ. To that purpose, we plug these bounds into equation 3.0), satisfied by the fixed points, which gives δ 3 λρ) λρδ δ 3 λρ = δ 3 λρ) δ where the last inequality holds for all sufficiently small δ, and 8δ 3 λρ) λρδ 8δ 3 λρ λρδ 8δ 3 λρ = 8 δ δ, Since the curve of fixed points is increasing see Figure ), the claim follows. Proposition 3.. The set I gets mapped by iterations of 3.9) into the set I 3 := [ 0, λ 3 ] {0}.. Proof. Recall that the fixed point x 0) is given by x 0) = λ 3 = δ λρ ) 3. Furthermore, we have [ ] [ 3 3 δ λρ δ), 8δ λρ + δ) 0, ) δ ] λρ 3 for all sufficiently small δ. Analogously to Lemma 3.9, we observe that, for all y [ 8 δ, 0), the map f, y ) is a contraction on [ 0, λ 3 ] to the fixed point x y ). Since x y ) [ 0, λ 3 ] for all y < 0, we can deduce that [ 0, λ 3 ] is invariant under f, y ) for all y [ 8 δ, 0 ). This implies the claim. The final transition is described in the following proposition. Proposition 3.. The set I 3 gets mapped by iterations of 3.9) into the set where µ = λρδ. [ I = minλρ, )δ, δ + µ ] Proof. We first consider the set S see Figure 6), defined by S := { } δ, { x, y ) 0 x x y ), 0 y } δ. Let x,0, y,0 ) be in S. Now, since f, y,0 ) is increasing and also the curve of fixed points is increasing, we have x,0 = fx,0, y,0 ) x y,0 ) x y,0 + νδ ) = x ȳ,0 ). Hence, the point x,0, ȳ,0 ) lies in S, as long as y,0 δ νδ. This means that trajectories leave S through the upper boundary at y = δ. Moreover, we consider the sets S := {x, y ) y x } S, S 3 := {x, y ) x < y λ x } S, 7

S S y y x x S 3 for λρ < ) S 3 for λρ ) y y x x Figure 6: The sets S, S and S 3 in the x, y )-plane. The shape of S 3 depends on the size of λρ. as depicted in Figure 6. For points x, y ) S 3, we have x = fx, y ) > x + νδ 3 λρ = x + νδ λ as well as ȳ = y + νδ. Thus they will be mapped into S S 3, as long as y δ νδ. Points on the left-hand boundary of S, given by {y = x }, clearly satisfy the same estimate such that they are also mapped into S S 3. By monotonicity of f, y ), see Lemma 3.7, this holds for all points in S. Hence, we deduce that trajectories leave S S 3 through the upper boundary at y = δ. In order to find a right-hand bound for the set S, we compute ) δ + µ λ ) δ + µ δ + µ λδ = δ, where we used that µ = λ δ 8 = λρ δ. Since the curve of fixed points is increasing see Figure ), this means that we can estimate x δ ) δ + µ. We observe from the proof of Proposition 3. that 0, λ ) 3 {0} is a subset of S S 3. Furthermore, it is now easy to see that S S 3 ) R { }) [ δ λρδ, δ + µ ] { } δ, when λρ <, 8

and S S 3 ) R This concludes the proof. { } [ δ δ, δ + µ ] { } δ, when λρ. Theorem 3.6 is now an immediate consequence of combining Proposition 3.0, Proposition 3. and Proposition 3.. 3.5 Dynamics in the exiting chart Transforming Σ out to the coordinates of the third chart gives κ 3 Σ out ) = { x 3 [ min λρ, ), + µδ ], r3 = ρ, ε 3 = 6δ, h 3 = ν }. Furthermore, we define Σ in 3 := { x 3 [ θ, + θ], r 3 = ρ, ε 3 = 6δ, h 3 = ν }, 3.) with θ = max{ λρ, }, such that κ 3Σ out ) Σ in 3 for sufficiently small δ. As already noted in section 3.3, we may assume, due to the controllably small error, that specific levels are hit by the trajectories and therefore, in particular, we assume ρ = ρ, where ρ was used to define + out 3.). Similarly to the situation in K, we consider the domain D 3 := {x 3, r 3, ε 3, h 3 ) R : r 3 [0, ρ], ε 3 [0, 6δ], h 3 [0, ν]} for the chart K 3, and we obtain the map x 3 = + h 3 ε 3 ) [x 3 + h 3 x 3 x 3) + λr 3 ε 3 )], r 3 = + h 3 ε 3 ) r 3, ε 3 = + h 3 ε 3 ) ε 3, h 3 = + h 3 ε 3 ) h 3, 3.) corresponding with 3.). Similarly to the system obtained in K, the special structure of 3.) again allows to explicitly determine solutions of the induced dynamical system in the components r 3, ε 3 and h 3. Lemma 3.3. For ζ 0 := h 3 0)ε 3 0) > 0, the trajectories of 3.) in r 3, ε 3, h 3 are given by for n N. h 3 n) = h 3 0) ζ 0 + n ζ0 + n ζ0, r 3 n) = r 3 0), ε 3 n) = ε 3 0) ζ 0 ζ 0 ζ 0 + n), Proof. Let ϑ := h 3 ε 3. Multiplying both sides of the equations for ε 3 and h 3 in 3.) yields ϑ = + ϑ) ϑ. 9

Solving this recursion for some initial condition ϑ0) = ζ 0 > 0 gives We use this observation to calculate n h 3 n) = h 3 0) k=0 + θk)) = h 3 0) ϑn) = ζ 0 + n. n k=0 θk) θk + ) = h 30) θ0) θn) = h 30) ζ 0 + n, ζ 0 and analogously ζ0 + n r 3 n) = r 3 0) ζ 0 as well as ζ0 ε 3 n) = ε 3 0) ζ 0 + n). This finishes the proof. We observe from 3.) that the hyperplane {ε 3 = 0} is an invariant set for the system 3.), foliating into the invariant lines {ε 3 = 0, r 3 = r3, h 3 = h 3} for all r3, h 3 0. Each of these lines has three fixed points, located at x 3, r 3, ε 3, h 3 ) = s, r3, 0, h 3), s {, 0, }. Linearizing around each of these, for h 3 > 0, we see that in x 3 -direction the fixed point at x 3 = 0 is unstable with eigenvalue + h 3 while those at x 3 = and x 3 = are stable with eigenvalue h 3 recall that h 3 hρ < < ). 3 Since, in our considerations, we enter K 3 via Σ in 3, our main interest lies in the family of stable fixed points at x 3, r 3, ε 3, h 3 ) =, r3, 0, h 3), corresponding with the branch S a + of the critical manifold for r3 > 0. Each of these fixed points has a three-dimensional center eigenspace as well as a one-dimensional stable eigenspace in x 3 -direction with eigenvalue h 3. The union of these fixed points forms an invariant manifold, which we call S + a,3 = {x 3, r 3, ε 3, h 3 ) D 3 : x 3 =, ε 3 = 0.} In particular, for each h 3 0 it contains the fixed point p + a h 3 ) =, 0, 0, h 3 ), which has gained hyperbolicity due to the desingularization of the origin. In analogy to Proposition 3. we get the following: Proposition 3.. The invariant manifold S a,3 + extends to a center-stable invariant manifold M a,3 + at p + a 0)) which is given in D 3 by a graph x 3 = l 3 r 3, ε 3, h 3 ) for a smooth mapping l 3. Furthermore, ε 3 is decreasing in D 3 and thereby, in particular, on M a,3), + whereas h 3, r 3 are increasing in D 3 and thereby on M a,3). + Proof. This is an immediate consequence of the considerations above and classical center manifold theory. Similarly to the proof of Proposition 3., the second claim follows from Lemma 3.3. Note that, on {r 3 > 0, ε 3 > 0, h 3 > 0}, the manifold M a,3 + corresponds to the union of the slow manifolds S a,ε,h +. We will follow the iterations of 3.) until the values r 3 = ρ, ε 3 = δ, h 3 = ν are reached. Lemma 3.3 tells us that this is the case when ζ 0+n =, which means that the number ζ 0 0

of iterations equals n = 3ζ 0 = 3. Let Π νδ 3 : Σ in 3 R be the transition map induced by n iterations of 3.) and consider the set Σ out 3 := { x 3 + ρ J, r 3 = ρ, ε 3 = δ, h 3 = ν }. Note that with this choice we have K 3 Σ out 3 ) = + out. The following Lemma summarizes the properties of the transition map Π 3. Lemma 3.5. For δ > 0 sufficiently small, we have Π 3 Σ in 3 ) Σ out 3, where Π 3 Σ in ) 3 νδ has a x 3 -width of at most C ν for some constant C intersections M a,3 + Σ in 3 as well as M a,3 + Π 3 Σ in 3 ) are non-empty sets. Proof. From the explicit solutions in Lemma 3.3, we directly see that Π 3 Σ in 3 ) {r 3 = ρ, ε 3 = δ, h 3 = ν}. 3 ) > 0. Furthermore, the Next, we consider the x 3 -equation of 3.) for ε 3 = 0 and h 3 > 0, which reads x 3 = x 3 + h 3 x3 x 3) ) =: gx 3, h 3 ). The cubic function g, h 3 ) is increasing between the stationary points at x 3 = ± +h3 h 3 =: ±κ. 3.3) Taking h 3 3, we achieve that κ > + θ. Since the fixed point at x 3 = 0 is unstable we can deduce that the set 0, κ] gets attracted to the fixed point at x 3 =, and so does [ θ, + θ] 0, κ]. We turn to giving estimates on the contraction rate towards x 3 = for different starting values in [ θ, + θ]. Note that the map g, h 3 ) is increasing and concave on 0, κ]. Since it is increasing, the subintervals [ θ, ) and, + θ] get mapped by g, h 3 ) into themselves. Due to its concavity, the function g, h 3 ) lies below the tangent at x 3 =, so that the contraction rate for values in, + θ] is at least as strong as the linear rate h 3 coming from the linearization around x 3 =. Concavity also yields that g, h 3 ) lies above the secant on the interval [ θ, ]. Thus the slope of the corresponding secant gives an estimate for the contraction rate of points in [ θ, ) to the fixed point at x 3 =. In more detail, linear interpolation between the points θ, g θ, h 3 )) and, ) yields the slope where g θ, h 3 ) θ) = [ θ) + h 3 θ) θ) ))] θ) θ) θ) ) = h 3 = h θ) θ θ ) 3 θ θ = h 3 θ) θ) = h 3 C θ), Note that 0 < C θ) < and C θ) 0 as θ. C θ) := θ) θ). 3.)

x3 θ Figure 7: The map g, h3) together with the diagonal x3 x3 dashed), the tangent at and the secant over [ θ, ]. We conclude that, on the interval [ θ, +θ], the map g, h3) is contracting with constant h3c. The transition map Π3 defined on Σ in 3 = { x3 [ θ, + θ], r3 = ρ, ε3 = 6δ, h3 = } ν is in x3-direction a perturbation of the n -fold of g, h3), as we now have r3, ε3 > 0. Also note that during the n = 3 iterations that define Π3, the variable h3 lies in νδ the interval [ ν, ν]. Thus, if we choose δ sufficiently small, we can achieve that Π3 is ) a contraction in x3-direction) with a rate of 3 C ν νδ. Consequently the image Π3Σ in 3 ) has a x3-width of at most θ ) 3 C ν νδ C ν Since for δ 0 the point M a,3 + { r3 = ρ, ε3 = 6δ, h3 = } ν approaches, ρ, 0, ν ), so especially its x3 coordinate approaches, we can deduce that M a,3 { + r3 = ρ, ε3 = 6δ, h3 = } ν lies in Σ in 3 for δ small enough. Moreover, since the manifold M a,3 + is invariant under the forward iterations of 3.), the point M a,3 + {r3 = ρ, ε3 = δ, h3 = ν} lies within Π3Σ in 3 ). Because M a,3 + {r3 = ρ, ε3 = δ, h3 = ν} converges to, ρ, 0, ν) as δ 0, we infer that for small δ the point M a,3 + {r3 = ρ, ε3 = δ, h3 = ν} and the exponentially small set Π3Σ in 3 ) around this point lie within Σ out 3 = {x3 + ρ J, r3 = ρ, ε3 = δ, h3 = ν}. 3.6 Blown-down dynamics As a last step, we will connect the individual results, obtained in each of the three charts, to prove Theorem 3.. ) 3 νδ.

Proof of Theorem 3. T) and T). Throughout the previous proofs we needed to make δ sufficiently small; now we choose δ 0 > 0 such that all of the statements hold true for every δ 0, δ 0 ]. This gives the value ε 0 := δ 0 ρ. As a consequence of condition 3.3) in the proof of Lemma 3.5, we assumed ν 3. Even stronger, Lemma 3.8. gives ν. For Lemma 3.5 we required ν + J { 5 Thus we can take h 0 := min, }. 5ρ ρ + 8 J In the previous sections we have seen that ) Σ in κ Σ out and 8ρ ). ) Σ in 3 κ 3 Σ out, 3.5) so that we can concatenate the three transition maps Π, Π and Π 3, using the appropriate coordinate changes in between, and define the map Π from Σ in to Σ out 3 Π := Π 3 κ 3 Π κ Π. In particular, we have seen from the analysis in the charts K, K and K 3 see Lemma 3.5, Theorem 3.6 and Lemma 3.5) that Ma 0, equalling Ma, 0 in K and Ma, 0 in K, is continued to Σ out 3 in the vicinity of M a,3. + Additionally, we have in = K Σ in ) and + out = K 3 Σ out 3 ). Hence, by reverting to the original coordinates x, y, ε, h), we obtain the transition map Π + from in to + out given by Π + := K 3 Π K. We have shown in Lemma 3.5 that the set Π 3 Σ in 3 ) Σ out 3 has a x 3 -width of ) C ν 3 νδ. Hence, by transforming back to the blown-down coordinates, it is easy to see that the width in x- direction of K 3 Π 3 Σ in 3 )) + out is O C h) K εh ) when we take K = 3ρ and C = C θ) ρ with θ and C from 3.) and 3.) respectively). Furthermore, Lemma 3.5 implies that K Σ in ) contains a point of the slow manifold Sa,ε,h 0. On the other hand, Lemma 3.5 shows that the exponentially small set K 3 Π 3 Σ in 3 )) contains a point of the slow manifold S a,ε,h +, and, due to 3.5) and the definition of Π+, we also have Π + in ) K 3 Π 3 Σ in 3 )) + out. This completes the proof of the statement in Theorem 3. for the case λ > 0. Finally, when λ < 0, observe that under the change of variables x x the x-equation in system 3.) gets transformed into which is equivalent to x = x + h xy x ) + λε ), x = x + h xy x ) λε ). Hence, the analysis is the same as for positive λ, with the same outcome under symmetric change of variables. 3

3.7 Canard Case The analysis in the case of λ = 0 may be carried out without a blow-up transformation. Hence, we treat the proof of Theorem 3. T3) separately here. Proof of Theorem 3. T3). One observes, that the system ) ) x x x + h xy x P : = ) )) y ȳ y + hε keeps {x = 0} invariant for all ε 0, since there is no dependency on ε in the x-equation. This means that we have S 0 a,ε,h = S 0 a and S 0 r,ε,h = S 0 r on the domain of our analysis see discussion around 3.3)) so that S 0 a,ε,h and S 0 r,ε,h are connected. The connecting trajectory γ starting in in is explicitly given by γn) = xn), yn) ) = 0, ρ + nhε ). The linearization along the trajectory γ is characterized by the variational equation ) + hyn 0 vn + ) = vn), vn) R, for all n N. 3.6) 0 While the fixed point w = 0, ) of 3.6) corresponds with the centre-direction along γ, the solution of 3.6) starting at v0) =, 0) corresponds with the transversal hyperbolic direction and can be explicitly solved to be n vn) = v n), v n)) = + h ρ + khε)), 0). k=0 Let us for simplicity assume that we have ρ εh =: N N. In particular, we set 0 out = {x, ρ ) x J}, where clearly γn) 0 out. We have already seen in Section 3.3 that due to the cubic structure the contraction rates in x-direction towards the locally stable fixed point for y < 0 are at least as strong as the linear rates achieved by linearization around x = 0. Similarly one also observes that the linear rates give a bound for the expansion rate for positive values of y. Hence, the contraction and expansion of trajectories from in to a neighborhood of 0 out can be estimated from above by the linear rate µ along the trajectory γ which satisfies N N µ = v N) = + h ρ + khε) + h ρ + khε) k=0 k=0 = N/ k=0 N/ k=0 + h ρ + khε) ) + hρ khε) ) = h ρ khε) ) Hence, we can give the bound µ N/ k=0 N/ k=0 h ρ khε) ) h ρ )) = h ρ )) N/. h ρ ) ) ρ εh <, meaning that the transition map Π 0 is contractive for the canard case since the contraction rates along S 0 a prevail over the expansion rates along S 0 r. Hence, the claim follows.

However, note that ) lim ρ h ρ h 0 ) εh =. Thus, as expected, in the limit h 0 one obtains the stability behaviour of the corresponding continuous-time system where contraction and expansion exactly compensate each other. It is still remarkable that the Euler method not only preserves the stability behaviour for trajectories close to the canard but even enhances stability as compared to the continuous-time case for sufficiently small h > 0. We already observed this surprising effect in the case of transcritical canards cf. []) but emphasize that in other similar situations, like the folded canard cf. [5]), the Euler method has clearly unfavourable stability properties. 3.8 Higher Order Terms We briefly discuss how our results can be generalized when higher order terms h x, y, ε) = O x y, xy, εx, εy, ε ) and h x, y, ε) = O x, y, ε) from.5) are included. The corresponding discretized dynamical system reads x x x + h xy x ) + λε + O x y, xy, εx, εy, ε ) ) y P : ε ȳ ε = y + hε + O x, y, ε)) ε. 3.7) h h h Note that due to the dependence of ȳ on x, points in the image of in under iteration of P will not share the same y-coordinate. Thus we cannot define the transition mappings Π ±,0 by just a fixed number of iterations of P, but instead pointwise for each initial value x 0, ρ ) in by Π ±,0 x 0, ρ ) = P n x 0 ) x 0, ρ ), where n x 0 ) = arg min distp n x 0, ρ ), {y = ρ }). n N We see that for a sufficiently small choice of ρ every trajectory will get εh close to {y = ρ }. In [] this had already to be taken into account for the case of a normal form without higher order terms. Firstly, let us discuss the problem for fixed λ 0. It is an important benefit of the blow-up method that in entering and exiting charts higher order terms have no significant impact. In more detail, we transform system 3.7) by K, proceeding as in Section 3.3, to obtain r = r h ε + h ε O r x ε, rε, rε)), This relation yields r x = r x + r h r x r r x ) + λr ε + O x r, x r 5, x r 5 ε, r 6 ε )). which simplifies to x = r r [ x + h x x ) + λr ε + O x r, x r, x r ε, r 3 ε ))]. 5

Consequently the transformed system in K can be written as x = h ε + O r ) ) [ x + h x x ) + λr ε + O r ) )], r = h ε + O r ) ) r, ε = h ε + O r ) ) ε, h = h ε + O r ) ) h. 3.8) For r = 0, this system is identical to 3.6). Hence, for sufficiently small r, we still obtain the existence of a center-stable manifold M 0 a, at the point p 0 a0) and the consequences thereof. A small choice of r means that we have to restrict ρ to sufficiently small values. For the exiting chart K 3, the situation is similar. In the rescaling chart, however, the higher order terms may not be bypassed that easily, but the strategy from Section 3. can be adapted. As in Section 3., we still have r = r and h = h. The remaining equations of 3.7) transform to and r x = r x + r h r x r y r x ) + λr + O r x y, r 5 x y, r 5 x, r 6 y, r 8 )), r ȳ = r y + r h r + O r x, r y, r ) ), which can be simplified and desingularized into x = x + h x y x ) + λr + O r x y, r x y, r x, r 3 y, r 5 )), 3.9) ȳ = y + h + O r x, r y, r ) ). The following arguments will not only require small ε but also sufficiently small ρ, so that the impact of normal form higher order terms can be controlled and the dynamics are determined by the remaining terms. Since the small parameter δ incorporates ρ and ε, it is more apparent in the original not blown up coordinates x, y, ε, h) how the choice of ρ determines the considered neighbourhood of the origin. Note that in the following small ρ and ε mean that the statements hold for sufficiently small fixed ρ and for all positive ε below some sufficiently small threshold. Restricting ρ and ε, we can assure upwarded movement in y -direction taking O ν δ ) steps to travel through the domain considered in the second chart. Our approach in Section 3. relied heavily on the curve of fixed points x y ) introduced in Lemma 3.7. In the more general setting involving the higher order terms, a curve corresponding to x y ) persists. In other words, one can show that for fixed, negative values of y there is exactly one positive fixed point for equation 3.9), given that ρ and ε and thus δ) are sufficiently small. Note that this can only be accomplished for y outside an interval of size O r), since otherwise the term x y is of order O x r) and therefore does not dominate additional terms of that order any longer. Using the curve of fixed points a result analogous to Proposition 3.0 can be shown, which implies that trajectories will enter the quadrant {x > 0, y < 0}. Equivalent statements to those of Lemma 3.8 can also be obtained for ρ and ε small enough. This means that for sufficiently small ν the mapping x x is monotone. Using the monotonicity one then easily checks that {x 0} is invariant under 3.9), simply by plugging in x = 0. In a similar manner one 6

can ensure that trajectories leave the rectangular set {x [0, δ ], y [ δ, δ ]} only through its upper boundary. Moreover, one can show that in the quadrant {x 0, y 0} trajectories are bounded away from the y -axis, independently from ε. The corresponding result is found in Proposition 3.. Hence, we may deduce that all trajectories will end up at a y -height close to δ with positive x -value, bounded away from 0 and smaller than δ. In other words, we obtain a result similar to Theorem 3.6. Transforming the exiting set into K 3 - coordinates allows to proceed similarly as in Section 3.5. As we have already seen in chart K, the higher order terms do not change the behaviour for sufficiently small ρ. Summarizing, we deduce that Theorem 3., T) and T), can be transferred to the general setting of 3.7). Furthermore, note that in the general case of the normal form.5) including higher order terms, the value of λ close to 0 giving a canard changes with the value of ε. For continuous time, this phenomenon is studied in detail for canards in folds in [0] and discussed for the transcritical and pitchfork case in [, Remark. and Remark.]. Using a Melnikov computation, one may show the existence of a function λ c ε / ) with λ c 0) = 0 such that for λ = λ c ε / ) the slow manifold Sa,ε extends to S r,ε + for sufficiently small ε. In Theorem 3. T3), we only treated the case λ c ε / ) 0 since we did not take into account perturbations from higher order terms. In order to obtain an analogous result to the ODE case, a treatment of the more general problem 3.7) about λ = 0 requires a discrete Melnikov computation, which is more complicated. Therefore, we are going to treat the general canard problem in the separate study [5]. 7