# Noninteger iterates of functions: How to get ODE from flow at a given time?

Suppose you have the autonomous ordinary differential equation $dx(t)/dt = f(x(t))$ with $x: \mathbb{R} \to \mathbb{R}$ and the initial condition $x(0)=x_0$. Then, assuming some regularity conditions, you get as solution the flow $\Phi(x_0,t):=x(t)$. To give a trivial example: If $f(x)=x$, then $\Phi(x_0,t)=x_0 \exp(t)$.

Now, I'm not interested in the trajectories for a given initial condition, that is in $\Phi(x_0,t)$ with $x_0$ fixed and $t$ variable; but in the map $x_0 -> \Phi(x_0, t)$ for a fixed $t$ (say $t=1$).

Given the function $f$, you can easily (at least in principle, by solving the ODE) get the function $\Phi(\cdot, 1)$. There are a lot of theorems about existence and uniqueness of this problem and analytical and numerical algorithms.

But how can one get $f$ out of $\Phi(\cdot, 1)$? Is this a well posed problem? Are there any theorems?

This problem is closely related to "interpolating" the $n$-fold functional iterates of $g$ (with $g^{} = \mathrm{Id}, g^{} = g, g^{} = g \circ g, g^{[n+m]} = g^{[n]} \circ g^{[m]}$ for $n,m \in \mathbb{N}$) from $n \in \mathbb{N}$ to real values. If such an interpolation succeeds, on can get the ODE out of the flow $\Phi(\cdot, 1)$ by determining $\Phi(\cdot, 1)^{[\alpha]}$ for small $\alpha >0$. I have done some calculation, that give results, but lack in rigor.

For noninteger iterates of functions, a classical reference is http://www.math-inst.hu/~p_erdos/1960-07.pdf.

## 3 Answers

If g is a real-analytic function defined near x0 with g(x0) = x0 and 0 < λ ≠ 1 where λ := g'(x0), then Koenigs proved that there exists a real-analytic homeomorphism h defined near x0 such that hgh-1(x) = L(x) where

L(x) := x0 + λ(x-x0)

and h is unique up to a constant factor*.

This allows g to be embedded in a local flow: Φ(x,t) := h-1 Lt h(x), where

Lt(x) := x0 + λt (x-x0),

such that Φ(x,1) = g(x) where defined.

Then Φ satisfies the differential equation dx(t)/dt = V(x(t)) where the velocity V (denoted by f in the Question) is given by

V(x) := ∂Φ(x,t)/∂t |t=0.

Note, however, that the flow Φ(x,t) into which the original function g embeds as the time-1 map need not be unique when there exist more than one fixed point of g. As an example, for x > 0 consider the function

gc(x) := cx.

Then for 1 < c < e1/e the function gc has two distinct fixed points each satisfying the hypotheses of the Koenigs theorem, and these give two distinct flows into which gc embeds as the time-1 map.

Concretely, set c = √2, so that gc(x) = x for both x = 2 and x = 4, with derivatives

gc'(2) = ln(2) and

gc'(4) = ln(4).

Calculating the respective real-analytic flows Φ2(x,t) and Φ4(x,t), both are defined for (x,t) = (3, 1/2).

But Φ2(3, 1/2) and Φ4(3, 1/2) first differ in the 25th decimal place. Hence they are solutions of distinct differential equations. I.e., they have different velocity functions.

* This is actually true in greater generality; see J. Milnor's book Dynamics in One Complex Variable, 3rd ed., Princeton University Press, 2006.

Being the time one map of a flow, there is an invariant one dimensional foliation where the flow "evolves" (assuming no singularities). For example, this implies that if there is a periodic point of $\phi(.,1)$ then there is an invariant circle. It is not hard to construct diffeomorphisms with periodic points and not having an invariant circle.

Even if there exists an invariant foliation of dimension one, and one can construct a flow for which the diffeomorphism is the time one map, uniqueness is not guarantied, as you can see by taking the time one map of a flow in $S^1$ with velocity $2\pi$ (this example can be made more interesting).

Another thing, which is related to what you say is that if you can take "square roots" of $\phi(-,1)$ indefinitely then one can recreate the ODE, but this is not true in general.

Aldrovandi and Freitas' article Continuous iteration of dynamical maps begins to lay the groundwork for viewing systems in physics as continuously iterated functions, providing an alternative to ODEs and PDEs. The function $f(x)$ that is continuously iterated describes the "physics" of $\Phi(x,t)$. If the scale of time $t$ is known, then setting $t=1$ trivially gives $f(x)=\Phi(x,1)$. If the scale of time $t$ is not known, then finding $f(x)$ probably becomes a difficult, if not intractable problem.

Consider a similar problem, but using cellular automata. Stephen Wolfram has searched for 25 years for a cellular automata rule that accurately models the laws of physics, but he hasn't released any results indicating he is able to model with fidelity any interesting physical process.

For more references see continuous iteration.