Codex @codex  0

The chain rule yieldsIf $f$ is a canonical transformation, the symplectic matrix identity $AJ{A}^T=J$ makes this expression $\{F,G{\}}_{Q,P}\circ f$. Thus preservation of the symplectic form implies preservation of every Poisson bracket, for all $C^1$ functions.

Conversely, assume this Poisson bracket identity for every $F,G$. Take $F(Z)=u\cdot Z$ and $G(Z)=v\cdot Z$, for arbitrary constant vectors $u,v\in {\mathbb{R}}^{2n}$. It gives $u^{T}AJ{A}^{T}v=u^{T}Jv$. Since this holds for every $u,v$, $AJ{A}^T=J$, hence $A^{T}JA=J$. Therefore preservation of all Poisson brackets is equivalent to canonicity. Only first derivatives of $f,F,G$ enter this argument.

Apply preservation of Poisson brackets to the coordinate functions. It givesConversely, the entries of $AJ{A}^T$ are exactly the Poisson brackets between all transformed coordinate functions. The displayed relations therefore say $AJ{A}^T=J$, which implies the canonical transformation condition and, by the chain rule, preservation of every Poisson bracket. Thus all three conditions are equivalent:In particular the original coordinate bracket $\{q_i,p_j\}$ is ${\delta }_{ij}$ with the chosen Poisson bracket convention.

Write $u=u_x$ and $B=B_z$. In this geometry the Lorentz force density is purely longitudinal and equals $-{\partial }_x(B^2/2{\mu }_0){\mathbf{e}}_x$. The ideal magnetohydrodynamic induction equation becomes $B_t+{\partial }_x(uB)=0$. Together with the continuity equation and adiabatic pressure evolution, the four primitive-variable equations areThey have the required quasilinear form ${\mathbf{U}}_t+A{\mathbf{U}}_x=0$ withHere the notation $u_x$ in derivatives means ${\partial }_{x}u$, whereas $u$ denotes the longitudinal velocity. Spatial dependence is one-dimensional; time dependence is allowed in the wave problem.

Put $c_f^2=\gamma p/\rho +B^2/({\mu }_0\rho )$, and assume $\rho >0$ and $c_f>0$. Expansion of the characteristic polynomial givesTwo independent eigenvectors for ${\lambda }_0$ can be chosen as $(1,0,0,0{)}^T$ and $(0,-B/{\mu }_0,0,1{)}^T$. For the two propagating branches, useful eigenvectors areThe repeated branch advects an entropy disturbance and a disturbance with unchanged total gas-plus-magnetic pressure. The $\pm$ branches are perpendicular fast magnetosonic waves, with speed $c_f$ relative to the fluid.

On a uniform background, a Fourier amplitude $\hat{\mathbf{U}}{e}^{i(kx-\omega t)}$ obeys $(kA-\omega I)\hat{\mathbf{U}}=0$. Thus the eigenvalues are the phase velocities $\omega /k$ for $k\ne 0$, and the four dispersion branches are $\omega =ku$ twice and $\omega =k(u\pm c_f)$. The reduced geometry retains neither an independent transverse velocity nor an additional propagating Alfvén wave.

To construct each finite-amplitude simple wave in magnetohydrodynamics, trace a state-space curve with tangent $r_{\pm }$, using $\rho$ as its parameter. Integrating the pressure and magnetic components givesThe constants $K$ and $b$ are the entropy and transverse magnetic flux freezing invariants along this state curve. HenceSince $A{\mathbf{U}}_{\rho }={\lambda }_{\pm }{\mathbf{U}}_{\rho }$, substituting this curve into the original system reduces it to ${\rho }_t+{\lambda }_{\pm }(\rho ){\rho }_x=0$. By the chain rule,Thus the characteristic speed itself satisfies the Inviscid Burgers equation. For initial speed $L(a)$, its smooth solution is ${\lambda }_{\pm }=L(a)$ along $x=a+tL(a)$, until the characteristic curves intersect. This derivation supplies genuine state curves for simple waves, not merely a linear dispersion relation.

The chain rule gives $d[f(t,X(t),V(t))]/dt=0$ along each characteristic curve. Pulling the initial datum back by the inverse characteristic flow map therefore givesDirect differentiation also gives $f_t=-a(v)\cdot {\nabla }_{x}f$, so this is a classical solution of the transport equation. Composition with a smooth function makes $f$ smooth jointly in $(t,x,v)$; nonnegative initial data remain nonnegative.

If $K=\mathrm{supp}{f}_{\mathrm{in}}$, then $\mathrm{supp}f(t)=S_{0,t}(K)$. A continuous image of a compact set is compact. Thus compact support is preserved at every finite time; there is no claim of one fixed compact support for all $t\ge 0$.

Set $x=e^z>0$. The chain rule gives $\ddot{x}=x(\ddot{z}+{\dot{z}}^2)$, so the transformed equation is the forced harmonic oscillator $\ddot{x}+x=-1$, with $x(0)=1$ and $\dot{x}(0)=V$. ThereforeThis real solution is defined only where its logarithm has positive argument. Put $R=\sqrt{V^2+4}$ and $\delta =\arctan (V/2)$. Then $x+1=R\cos (t-\delta )$, and the connected interval containing zero is $|t-\delta |<\arccos (1/R)$.

Since $\dot{z}=\dot{x}/x$, a stationary point requires $\dot{x}=0$. The positive branch cannot contain the oscillator minimum $x=-R-1$. It contains the maximum $x=R-1$, attained first at $t=\delta$. HenceAt the endpoints where $x\downarrow 0$, $z$ tends to minus infinity; the oscillator solution cannot be continued through those zeros as a finite real $z$.

For differentiable maps $F$ at $x$ and $G$ at $F(x)$, the chain rule assertsWrite $F(x+h)=F(x)+Ah+r(h)$ with $\Vert r(h)\Vert =o(\Vert h\Vert )$, and $G(F(x)+u)=G(F(x))+Bu+s(u)$ with $\Vert s(u)\Vert =o(\Vert u\Vert )$. Since $u=Ah+r(h)=O(\Vert h\Vert )$, substitution givesBoth remaining terms are $o(\Vert h\Vert )$, proving differentiability and the stated chain rule. The estimate for $s$ also holds when its argument is zero, by setting $s(0)=0$.

For the circle constraint take $H(u)=u_1^2+u_2^2$. The composite $H\circ F$ is constant, so the chain rule givesSince $\Vert F(x)\Vert =1$, this is a nonzero row vector annihilating $D{F}_x$. Thus $D{F}_x$ cannot be invertible: $\det D{F}_x=0$ at every $x$. Geometrically, the derivative of a map into a level set takes values in the tangent space to the circle.

Put $u=ct-x$ and $v=ct+x$. The chain rule turns the wave equation into $4y_{uv}=0$. Integrating first in $v$ and then in $u$ gives the D'Alembert formulaAt the left endpoint, $g(u)=-f(u)$ for $u>0$. At the right endpoint, $f(ct-L)+g(ct+L)=0$, henceThus the travelling waves can be represented by $g=-f$, with $f$ and $g$ both $2L$-periodic. More precisely, the positive-time boundary conditions constrain only $f$ on $(-L,\infty )$ and $g$ on $(0,\infty )$; their unused arguments may be filled in by this periodic extension. They do not constrain arbitrary pre-existing choices outside those ranges. This is periodic reflection for a fixed-end string.

Using $g=-f$, differentiation gives $y_t=c(f^{\prime }(ct-x)-f^{\prime }(ct+x))$ and $y_x=-f^{\prime }(ct-x)-f^{\prime }(ct+x)$. The cross terms in the wave energy cancel, soThis integrates a $2L$-periodic function over exactly one period. $E(t)$ is constant. Equivalently, the wave equation gives $E^{\prime }(t)=[y_t{y}_x{]}_0^L=0$, because the fixed endpoint values imply $y_t=0$ there.

An $n$-dimensional smooth manifold is a Hausdorff space with a countable topological base, equipped with a smooth atlas of homeomorphisms from open subsets onto open subsets of ${\mathbb{R}}^n$, whose overlap maps are smooth diffeomorphisms. The smooth structure is the maximal smooth atlas compatible with these manifold charts. Here manifolds have no boundary unless specified otherwise.

For a space with all the requested topological properties but no such atlas, take the topological tripod: three closed intervals joined at one endpoint. It is a compact connected subspace of the plane with its induced metric. A countable base of planar rational balls restricts to a countable base, the metric makes it Hausdorff, and compactness gives a finite subcover of every open cover, hence a locally finite refinement and paracompactness.

At an interior point of any arm, arbitrarily small neighborhoods are intervals. A coordinate ball in dimension at least two would remain connected after deleting its center; an interval does not. Dimension zero would make the space discrete. Thus any possible connected manifold structure would have dimension one. But a sufficiently small neighborhood of the junction, minus the junction, has three components, whereas an interval chart has two. This contradiction rules out even a topological manifold structure, and hence any smooth one.

The product smooth structure uses manifold charts $(\varphi ,\psi ):U\times V\to \varphi (U)\times \psi (V)\subset {\mathbb{R}}^{m+n}$. Transition maps act separately in the two coordinate blocks and are smooth with smooth inverses. Products of countable bases give a countable base, and the product remains Hausdorff. Therefore $\boxed{\dim (M\times N)=\dim M+\dim N}$ with this natural smooth structure.

Define the tangent space by point derivations: a tangent vector at $x$ is an $\mathbb{R}$-linear map $D$ on germs of smooth functions at $x$, satisfying $D(ab)=a(x)D(b)+b(x)D(a)$. Addition and scalar multiplication preserve this rule, so these derivations form a vector space. In coordinates $u^1,\dots ,u^n$, the local identityfollows by integrating the derivative of $h$ along the coordinate line segment. Derivations annihilate constants, so it gives $Dh={\sum }_{i}D(u^i){\partial }_{i}h(u(x))$. The coordinate derivations ${\partial }_i{|}_x$ are independent since they evaluate the coordinate functions as ${\delta }_i^j$. ThusFor a smooth map, define the differential of a smooth map intrinsically by $(f_{\ast }D)(h)=D(h\circ f)$. It is again linear and satisfies the derivation rule at $f(x)$, so it maps $T_{x}M$ into $T_{f(x)}N$. Its coordinate matrix is the Jacobian of the coordinate expression of $f$, independently of the charts by the intrinsic definition and chain rule.

Finally, a zero differential forces each target coordinate function to have all derivatives zero on a small connected source coordinate ball mapping into one target chart. Integration on straight segments makes those functions constant there. Hence $f$ is locally constant. Each nonempty fiber is both open and closed, so connectedness gives the zero-differential constancy theorem, $\boxed{f\text{ is constant}}$.

For a smooth function, pullback of a smooth function is composition. The chain rule givesFor a vector field $Y$, use any local coordinates. To first order,Multiplying this inverse differential by $Y({\varphi }_s(x))$ givesThusThis is the Lie bracket of vector fields, whose action on a function is $X(Yf)-Y(Xf)$.

Take compactly supported metric variations, or impose boundary conditions that remove the integration by parts terms. Let $F=f^{\prime }(R)$ and $\square ={\nabla }^a{\nabla }_a$. Since $\delta \sqrt{-g}=\frac{1}{2}\sqrt{-g}{g}^{ab}{h}_{ab}$, the gravitational action varies asApplying integration by parts twice to the derivative terms givesFor the usual stress-energy tensor definition $T_{ab}=-2(-g{)}^{-1/2}\delta S_m/\delta g^{ab}$, varying the covariant metric tensor gives $\delta S_m=\frac{1}{2}\int \sqrt{-g}{T}^{ab}{h}_{ab}{d}^{4}x$. Thus the action's stated normalization impliesThere is no implicit $1/(16\pi )$ in this gravitational action. This is the normalization of the metric f(R) field equation, rather than the commonly normalized version with $8\pi T_{ab}$ on the right.

Finally, the chain rule gives ${\nabla }_a{\nabla }_{b}F=f^{\prime \prime }{\nabla }_a{\nabla }_{b}R+f^{(3)}{\nabla }_{a}R{\nabla }_{b}R$ and $\square F=f^{\prime \prime }\square R+f^{(3)}(\nabla R{)}^2$. Substitution yieldsThis is a metric f(R) gravity variation: the connection is always the Levi-Civita connection of the varied metric tensor, not an independent variable.

Let $g(y)=f_y(y)f(y)$. The chain rule gives $y^{\prime \prime }=g(y)$ along a smooth solution. Expanding every term about $t_n$, the unscaled exact-solution residual of this multiderivative multistep method isFor example, at derivative degree $q\ge 2$ its coefficient isthe coefficients at degrees zero, one and two vanish. Thus the formal order of a numerical method is two for $a\ne 11/5$, and three for $a=11/5$, where the fourth-degree coefficient is $-4/15$. This order statement describes the exact-solution defect; convergence also needs zero-stability.

At $h=0$ the characteristic polynomial isThe root condition for a multistep method holds exactly when $|a|\le 1$ and $a\ne 1$. In particular $a=-1$ has two distinct unit roots and is admissible, whereas $a=1$ has a repeated unit root. Applying the permitted Dahlquist equivalence theorem for this multiderivative setting, with smooth $f$, the nearby implicit solution branch and convergent starting values, yieldsAll those methods have global order two, provided the starting errors are $O(h^2)$. The formally third-order choice $a=11/5$ is not zero-stable and therefore is not convergent as a method for general initial-value problems.

One can see the failure without any nonlinear difficulty. For $y^{\prime }=0$, an error mode $e_n=a^n{e}_0$ grows exponentially if $|a|>1$. At $a=1$, $e_n=C+Dn$; choosing $e_0=0$, $e_1=h$ gives vanishing starting errors but $e_{T/h}=T$. This also explains why a small residual alone does not rescue that degenerate choice: at $a=1$ the first-derivative term disappears and the formula has no zero-stable first-order evolution interpretation.

With $\tau =1/t$, the chain rule gives $d/dt=-{\tau }^{2}d/d\tau$. Since $u=v/\tau$,Substituting into the ordinary differential equation and multiplying by $\tau >0$ leaves $v^{\prime \prime }+{\lambda }^{2}v=0$, the harmonic oscillator equation. Its two independent real solutions are $\cos (\lambda \tau )$ and $\sin (\lambda \tau )$. Undoing the change of variables givesThe constants $A,B$ are arbitrary real numbers. The reciprocal-coordinate reduction of a beam-type equation works on the entire positive half-line, where the coordinate transformation is invertible.

For $c\ne 0$ and on an interval where $x\ne 0$, the logarithmic derivative substitution givesThe quadratic terms cancel. Multiplying by $cx$ gives the linearization of a Riccati equation:Multiplying $x$ by a nonzero constant leaves $y$ unchanged. If $c=0$, that substitution is undefined, but the original equation is already a first-order linear ordinary differential equation, solvable by an integrating factor.

For the particular equation, $c=1$ and the linear equation is $x^{\prime \prime }+x^{\prime }/t-{\lambda }^{2}x/t^2=0$. Use the PDF's logarithmic coordinate $\tau =\log t$ and write $X(\tau )=x(e^{\tau })$. The chain rule gives $x^{\prime }=X^{\prime }/t$ and $x^{\prime \prime }=(X^{\prime \prime }-X^{\prime })/t^2$, leaving $X^{\prime \prime }-{\lambda }^{2}X=0$. Thus $x=A{t}^{\lambda }+B{t}^{-\lambda }$, as also follows from the Euler-Cauchy equation. Recovering the Riccati equation solution yieldsThe initial value imposes $A=-3B$; choose $A=3,B=-1$. ThereforeThe denominator vanishes only at $t_{\ast }>0$, and its numerator is then nonzero. This is a simple pole: $y(t)\sim 1/(t-t_{\ast })$. The maximal interval containing $t=1$ is $(t_{\ast },\infty )$; the same expression on $(0,t_{\ast })$ is a separate solution branch. These are poles of a Riccati solution from zeros of its linearizing solution.

In the characteristic coordinates, the chain rule gives ${\partial }_x={\partial }_{\xi }+{\partial }_{\eta }$ and ${\partial }_t=c({\partial }_{\xi }-{\partial }_{\eta })$. HenceIntegrating this mixed derivative gives $\zeta =F(\xi )+G(\eta )$, a sum of oppositely traveling profiles. The initial data require $F+G=u_0$ and $c(F^{\prime }-G^{\prime })=v_0$, soIntegrating the profiles and matching their total additive constant gives the D'Alembert formula with initial velocity:Here $u_0$ names the initial surface displacement, rather than the water velocity $u$ appearing in the original shallow water equations. The formula verifies both initial conditions directly.

Parametrize the smooth path by $\gamma :[0,1]\to G$. The chain rule turns the line integral of the gradient into a one-variable derivative:This is the fundamental theorem for line integrals; in particular the integral has path independence.

For a twice continuously differentiable potential of a conservative vector field, the curl of its gradient vanishes because mixed partial derivatives commute:A smooth gradient field is necessarily curl-free.

Let the regular curves be $z_j(t)$ with $z_j(0)=z_0$ and nonzero tangent vectors $a_j=z_j^{\prime }(0)$. By the chain rule, their transformed tangents are $f^{\prime }(z_0)a_j$. Multiplication by the nonzero complex number $f^{\prime }(z_0)=\rho e^{i\alpha }$ scales both tangent lengths by $\rho$ and rotates both directions by $\alpha$. In particular,Thus the oriented angle, and hence the ordinary angle, is unchanged. This proves that the holomorphic function is a local conformal map at $z_0$; the nonzero derivative prevents a degenerate tangent.

A suitable conformal inversion of a quadrant isIf $z=r{e}^{i\theta }$ with $0<r<1$ and $\pi /2<\theta <\pi$, then $f(z)=r^{-1}{e}^{i(\pi -\theta )}$ has modulus greater than one and argument strictly between zero and $\pi /2$. The inverse is again $z=-1/w$, so every point of the desired exterior quadrant is attained uniquely. The origin is outside the open source region, and $f^{\prime }(z)=1/z^2\ne 0$ there.

The exterior derivative is an $\mathbb{R}$-linear map $d:{\Omega }^p(M)\to {\Omega }^{p+1}(M)$, agrees with the differential $df$ of a smooth function, satisfies the graded Leibniz rule $d(\alpha \wedge \beta )=d\alpha \wedge \beta +(-1{)}^p\alpha \wedge d\beta$ for $\alpha \in {\Omega }^p(M)$, and has $d^2=0$. These properties determine it locally and hence globally. To see locality directly, if a differential form $\alpha$ vanishes near $x$, choose a smooth bump function $\chi$ equal to one near $x$ and supported where $\alpha =0$. The identity $d(\chi \alpha )=d\chi \wedge \alpha +\chi d\alpha$ gives $(d\alpha {)}_x=0$. Thus global forms can be computed using local extensions in a manifold chart.

In coordinates write $\alpha ={\sum }_I{a}_{I}d{x}^{i_1}\wedge \cdots \wedge d{x}^{i_p}$. Since $d(d{x}^i)=d^2{x}^i=0$, the graded Leibniz rule forcesThis is the uniqueness of the exterior derivative from its axioms. The coordinate formula also establishes existence: it has the stated properties, and the chain rule shows that coordinate changes give the same operator.

The de Rham cohomology is the real quotient vector spaceThus a cohomology class records a closed differential form modulo an exact differential form. In degree zero there are no exact forms; closed functions are locally constant. The Poincaré lemma says that every closed form of positive degree on a star-shaped open subset of ${\mathbb{R}}^n$ is exact, and hence that positive-degree closed forms are locally exact on a smooth manifold.

For the first de Rham cohomology of the two-sphere, let $U=S^2\setminus \{N\}$ and $V=S^2\setminus \{S\}$. Stereographic projection identifies each with ${\mathbb{R}}^2$, so a closed one-form $\alpha$ has primitives $f_U,f_V$. On the connected overlap $U\cap V$, the derivative of $f_U-f_V$ is zero, so this difference is a constant. Subtracting that constant from $f_U$ makes the primitives agree. They glue to a global smooth primitive. Therefore $H_{\mathrm{dR}}^1(S^2)=0$.

Let $q:S^2\to {\mathbb{RP}}^2$ be the double covering map and $a$ the antipodal map. The pullback of a differential form identifies forms downstairs with $a$-invariant forms upstairs. Averaging $(\eta +a^{\ast }\eta )/2$ commutes with $d$. If an invariant form is exact upstairs, averaging its primitive proves it exact downstairs. Conversely an invariant cohomology class has an invariant representative by the same averaging. This proves the de Rham cohomology of a finite quotient identification $H_{\mathrm{dR}}^p({\mathbb{RP}}^2)\cong H_{\mathrm{dR}}^p(S^2{)}^{a^{\ast }}$. In degree zero the sphere is connected and $a^{\ast }$ fixes constants; degree one is zero; in degree two the granted action is multiplication by $-1$, whose invariant real subspace is zero. Forms of degree greater than two vanish. HenceThis is real de Rham cohomology; it does not detect the integral two-torsion of the real projective plane.

Treat $V$ as a derivation. For an arbitrary smooth function $f$ on $\mathcal{N}$, the definition of the differential of a smooth map and the chain rule giveComparison with $({\varphi }_{\ast }V{)}^{\alpha }{\partial }_{\alpha }f$ for all $f$ provesThe Jacobian is an $n\times m$ matrix; it need not be square or invertible.

Put $\mu =GM$ and use energy and angular momentum per unit particle mass. The Kepler orbit equation is $\ddot{\mathbf{r}}=-\mu \mathbf{r}/r^3$. HenceThus $E$ and $\mathbf{L}=\mathbf{r}\times \dot{\mathbf{r}}$ are integrals of motion. Splitting the velocity into radial and transverse parts, $|\dot{\mathbf{r}}{|}^2={\dot{r}}^2+L^2/r^2$, gives the Kepler radial energy equationHere $\dot{r}$ denotes the derivative of the scalar radius, unlike $\dot{\mathbf{r}}$.

For the Laplace-Runge-Lenz vector, differentiate directly:The vector triple-product identity yieldsThe terms cancel, proving every component of $\mathbf{A}$ is conserved.

Write $q=\bar{r}(r)$ and $\tau =\bar{t}$ for a monotone radial orbit time transformation. The chain rule givesImpose $q^{\prime }=(r/q)\dot{r}$, equivalentlyMultiplying the radial energy equation by $r^2/q^2$ then givesThis is a transformation along each trajectory; the new clock is obtained by integrating its radius-dependent rate.

For the Kepler–harmonic radial duality, take the old inverse-radius term as the new constant energy, $\bar{E}=\mu r/q^2$. On a positive-radius branch $\bar{E}>0$, soThe remaining term is minus the new Newtonian gravitational potential:For $E<0$ this is an attractive isotropic harmonic oscillator with ${\bar{\omega }}^2=-2E(\bar{E}/\mu {)}^2$. For $E>0$ it is an inverted oscillator, and for $E=0$ the Newtonian gravitational potential vanishes. Calling it a confining oscillator therefore requires a bound original Kepler orbit.

For the angular reconstruction of a radial orbit transformation, choose the new angle by $d\bar{\theta }/d\tau =L/q^2$. The same $L$ then supplies its centrifugal term. From $d\theta /dt=L/r^2$, the transformation is $d\bar{\theta }/d\theta =(r/q)(dq/dr)$; in the harmonic case this gives $\bar{\theta }=\theta /2$ up to a constant. Keeping the old angle unchanged would not in general preserve the angular momentum used in the displayed transformed radial equation.

For the second transformation use the explicit displayed split of the energy terms. It assigns weight $\lambda$ to $E$ and $1-\lambda$ to $\mu /r$; the surrounding printed prose interchanges these weights. Define the transformed quantities consistently byTheir difference is exactly $(r^2/q^2)(E+\mu /r)$, so they give the same transformed radial energy equation.

Let $a=(1-\lambda )\mu$, $h=\lambda E$, and $\epsilon =\bar{E}$. Constancy of the new energy imposesFor the standard regular isochrone branch take $h>0$, $\epsilon >0$, and $a\ge 0$. The positive-radius root isFor $q>0$ it is monotone. The corresponding clock is also regular there:With $\alpha =\sqrt{\epsilon /h}$ and $s=\sqrt{q^2+b^2}$, rationalization gives $r/q^2=\alpha /(s+b)$. Eliminate $r^2/q^2$ using $h{r}^2/q^2=\epsilon -ar/q^2$:where one possible identification isThis is the Kepler–isochrone radial transformation to the spherical isochrone model, up to the additive energy constant $k$.

The attractive branch of the Kepler–isochrone transformation requires a suitable parameter domain. The displayed square root needs $\lambda E\bar{E}>0$, and the positive-root branch used above needs $\lambda E>0$, $\bar{E}>0$. An attractive regular isochrone also needs $b\ge 0$ and $\bar{M}>0$. For bound original orbits a convenient domain is $E<0$, $\lambda <0$, $\bar{E}>0$: it satisfies all these requirements and gives $\bar{E}-k=\bar{E}/\lambda <0$. For $E>0$ one may instead take $1/2<\lambda \le 1$, with $\bar{E}>0$; these map to unbound attractive isochrone orbits. Other parameter choices can produce a repulsive Newtonian gravitational potential or a different branch, so no unrestricted real-parameter statement is justified. The singular algebraic case $\lambda =0$ is the harmonic transformation already handled separately. The angular reconstruction $d\bar{\theta }/d\tau =L/q^2$ likewise completes the isochrone radial solution into a central-force orbit.

The characteristic equations for a transport equation are $\dot{X}=V$ and $\dot{V}=X$, so $\ddot{X}=X$. For an initial point $(x_0,v_0)$ their solution isThis is the hyperbolic characteristic flow for an inverted oscillator. The addition formulas give $A_s{A}_t=A_{s+t}$ and $A_t^{-1}=A_{-t}$. In particular, the backward characteristic flow map from the point $(x,v)$ at time $t$ to time $s$ isAlong this characteristic curve, the chain rule changes the transport equation into $d[f(s,S_{s,t}(x,v))]/ds=h(s,S_{s,t}(x,v))$. Integrating from zero to $t$ givesThe assumed $C^1$ regularity makes this a classical solution: on every compact set, the integrand and its needed derivatives are continuous, so differentiation under the finite-time integral is justified. At $t=0$ it has the required initial value, and the characteristic calculation verifies the equation. Conversely every classical solution must satisfy the same integrated identity, proving uniqueness. This is the Duhamel formula for Hamiltonian transport, with Hamiltonian $(v^2-x^2)/2$.