---
abstract: 'We show a general method allowing the solution calculation, in the form of a power series, for a very large class of nonlinear Ordinary Differential Equations (ODEs), namely the real analytic $\sigma\pi$-ODEs (and, more in general, the real analytic $\sigma\pi$-[*reducible*]{} ODEs) in many indeterminates, characterized by an ODE-function given by generalized polynomials of the indeterminates and their derivatives, i.e. functions formally polynomial with exponents, though the exponent can be any real number. The solution method consists in reducing the ODE to a certain canonical quadratic ODE, with a larger number of indeterminates, whose solutions include, as sub-solutions, the original solutions, and for which a recursive formula is shown, giving all coefficients of the solution power series directly from the ODE parameters. The reduction method is named [*exact quadratization*]{} and was formerly introduced in another our article, where we considered explicit ODEs only. In the present paper, which is self-contained to a large extent, we review and complete the theory of exact quadratization by solving issues, such as for instance the [*globality*]{} of the quadratization, that had remained open, and also extend it to the more general case of an implicitly defined ODE. Finally, we argue that the result can be seen as a partial solution of a differential version of the 22nd Hilbert’s problem.'
address: '$^{1}$ Istituto di Analisi dei Sistemi ed Informatica “Antonio Ruberti”, Consiglio Nazionale della Ricerche, Via dei Taurini 19, Roma, Italy. '
author:
- 'Francesco Carravetta$^{1}$'
title: On the Solution Calculation of Nonlinear Ordinary Differential Equations via Exact Quadratization
---

Introduction {#sc:int}
============

Let us consider the following quadratic ordinary differential equation (ODE) in $\re^m$, written component-wise ($i=1,\ldots,m$): x\_i=\_[j=1]{}\^m v\_[i,j]{}x\_ix\_j,\[Quad1\] where the ’time’ $t\in\cal T\subset \re$, $\cal T$ an open interval, (resp: $x\in\re^m$) is the independent (resp: dependent) variable, the dot denotes time-derivative, $x_i$ denotes the $i$th component of the vector $x$, and $v_{i,j}$ are real functions of $t$ analytic on $\cal T$. Whatever we fix a point $(x,s)\in \re^n\times \cal T$ we know that there exists an unique solution $t\mapsto x(t)$ of (\[Quad1\]) defined and analytic on some interval ${\bf I}_{x,s}\subset \cal T$, which in general depends of the point $(x,s)$, with $s\in {\bf I}_{x,s}$ and $x(s)=x$. Moreover, since the solution $x(t)$ is analytic, for any $t_0\in {\bf I}_{x,s}$ it can be written, component-wise, as the power series: x\_i(t)=\_[k=0]{}\^c\_k(i) [(t-t\_0)\^kk!]{},c\_k(i)=x\_i\^[(k)]{}(t\_0), \[SPSDint\] where $x^{(k)}$ denotes the $k$th time-derivative, and (\[SPSDint\]) holds at least whereby all series, for $i=1,\ldots,n$, converge, which surely occurs for all $t\in {\bf I}(t_0,r) \subset 
{\bf I}_{x,s}$, an open interval centered in $t_0$ with ray $r\in \re^+\cup \{+\infty\}$.\
Now, let us pose the following question:\
[*(q)*]{} is there a recursive (in $k$) formula giving the coefficients $c_k(i)$ of the power series (\[SPSDint\]) for all $k\in \na$ directly from the coefficients $v_{i,j}$?\
The above question has an interest going far beyond the single case study, insofar one is able to show that the quadratic ODE (\[Quad1\]) [*is not*]{} really a [*particular*]{} differential equation, and it is in fact a sort of [*nonlinear paradigm*]{}, a basic case to which a large class of analytic ODE can be always [*reduced*]{}, in the sense we are going to explain below.\
We say that a first order ODE of the type: x=f(x,t),\[S1\] where $f:\cal A\to \re^n$, is a function defined on the open set $\cal A\subset \re^{n+1}$ and locally Lipschitz therein, is [*globally, exactly*]{} (adjectives that hereinafter we omit whether clear from context) [*quadratizable*]{}, or that [*undergoes a (global, exact) quadratization*]{}, if there exist an integer $m\ge n$, and locally Lipschitz functions $v_{i,j}: \cal T\to \re$ – with $\cal T$ an open interval of $\re$ and $i,j=1,\ldots,m$ – such that every solution is a [*local sub-solution*]{}[^1] of a solution of a quadratic ODE of the type (\[Quad1\]), which is then accordingly referred to as [*a quadratization*]{} of (\[S1\]).

Before going more in depth, we clarify some points of the definition given above. We define the [*maximal domain*]{}, $\cal A\subset  \re^{n+1}$, of (\[S1\]) to be the largest (open) set such that for any $(x,t)\in\cal A$ there exist a unique solution of (\[S1\]) whose graph passes through $(x,t)$. We always understand that $\cal A$ is non void. If $(x,t)\in \cal A$, we denote by ${\bf I}_{x,t}\subset \re$ the largest (open) interval where the solution $s\mapsto x(s)$, such that $x(t)=x$ is defined, and identify the solution with the point $(x,t)$, or even with the sole $x$, when $t$ is understood, obvious, or un-essential [^2]. We say that $(x,t)\in \cal A$ – solution of the ODE (\[S1\]) – is a local (resp: global, or no adjective) sub-solution of a second ODE, of the same type (\[S1\]) but $m$-dimensional, with $m\ge n$, if it is equal, point-wise in ${\bf I}\subset {\bf I}_{x,t}$, with ${\bf I}$ an open sub-interval, (resp: point-wise in ${\bf I}_{x,t}$) to $n$ out of the $m$ components of a solution of the second ODE. If we collect the indices of the latter $n$ components in a set $\cal I^*\subset \{1,\ldots,m\}$ we can say that $(x,t)$ is a sub-solution of the second ODE [*with respect to*]{} $\cal I^*$.

Let us consider a map $\Phi:{\rm Dom}\{\Phi\}{\rm (open)}\subset {\cal A}\to \re^m$, satisfying the following property: for any $(x,t) \in {\rm Dom}\{\Phi\}$ there exists an open neighborhood $U\subset {\rm Dom}\{\Phi\}$ of $(x,t)$ such that, for any $(x',t')\in U$ belonging to the graph of the solution $(x,t)$ of (\[S1\]), $(\Phi(x',t'),t')$ belongs to the graph of the solution $(\Phi(x,t),t)$ of (\[Quad1\]). We call such a map: [*locally compatible* ]{} (with the ODEs (\[S1\]) and (\[Quad1\])), [^3] and the ODE (\[S1\]) will be said [*(globally, exactly) quadratizable under $\Phi$*]{} if every solution $(x,t)$ of (\[S1\]) is a local sub-solution of the type $(\Phi(x^*,t^*),t^*)$ – with $(x^*,t^*)\in {\rm Dom}\{\Phi\}$ – of a quadratic ODE of the type (\[Quad1\]). If (\[S1\]) is quadratizable under $\Phi$, we always understand in the following that $\Phi$ is a locally compatible map.

If (\[S1\]) is quadratizable under $\Phi$, and any of its solutions $(x,t)$ is an analytic function on ${\bf I}_{x,t}$ — which is the case as $f$ in $(\ref{S1})$ is analytic in all of its arguments — then any solution $(x,t)$ of (\[S1\]) whose graph crosses ${\rm Dom}\{\Phi\}$, say $(x(t^*),t^*)=(x^*,t^*)\in {\rm Dom}\{\Phi\}$, can be recovered (by analytic prolongation) by the solution $(\Phi(x^*,t^*),t^*)$ of (\[Quad1\]). In other words: we can [*calculate*]{} at least all solutions whose graph is not disjoint with the domain of $\Phi$, provided that we are able to calculate, even only locally, the solutions of a quadratic ODE of the type (\[Quad1\]).

In the paper [@Ca15] we have given a first answer to the question about which kind of ODEs can be exactly quadratized. The basic result (stated in a rather different setting than the one we use here) has been that all ODEs of the type (\[S1\]) with $f$ given component-wise by: f\_i(x,t)=\_[l=1]{}\^[\_i]{} v\_[i,l]{}(t) X\_[i,l]{}(x); X\_[i,l]{}(x)=\_[j=1]{}\^n x\_j\^[p\^l\_[i,j]{}]{},\[SP-form1\] where $p^l_{i,j}$ are real exponents, $\nu_i\in\na$ (if $\nu_i=0$ then $f_i=0$) and $v_{i,l}$ are analytic functions all defined on some open real interval, namely $\cal T$, undergo a [*dense*]{} quadratization (we clarify in a moment what ’dense quadratization’ means) under the following real analytic map[^4] $\Phi:{\rm Dom}\{\Phi\}\subset {\cal A}\to \re^m$ – $\cal A\subset \re^n\times \cal T$ being the (open) analyticity domain of (\[SP-form1\]), and $m=\nu_1+\ldots+\nu_n$, thus $m\ge n$ – which can be written component-wise as follows: \_[i,l]{}(x\_1,…,x\_n)=x\_i\^[-1]{} X\_[i,l]{}. \[ChV\]

In (\[SP-form1\]) the function $f$ is a [*formal polynomial*]{} of $x$, that is to say: a polynomial writing where the exponents are allowed to be [*any real number*]{}. In [@Ca15] we named this kind of functions [*$\sigma\pi$-functions*]{}, and $\sigma\pi$-[*ODE*]{} the associated differential equation.

A [*dense*]{} quadratization under a map $\Phi$ differs from the definition of (global) quadratization (under $\Phi$) given before in that it is not required that [*every*]{} solution-graph crosses ${\rm Dom}\{\Phi\}$, but rather that the set of all solution-graphs that do not cross it form a zero-measure set in $\re^{n+1}$, i.e. the set $\cal A\cap {\rm Dom}\{\Phi\}$ has zero-measure, or, which is the same, ${\rm Dom}\{\Phi\}$ is a [*dense*]{} subset in $\cal A$. Indeed, it can be at once inferred by (\[ChV\]) that ${\rm Dom}\{\Phi\}$ is $\cal A\setminus S$ where $S$ denotes the union of the $n$ hyperplanes in $\re^{n+1}$ described by the equation $x_i=0$, for $i=1,\ldots,n$, so, what has been proven in [@Ca15] can be stated in other words as follows: every solution of a $\sigma\pi$-system, except possibly a sub-set whose curves all lye in some (among the first $n$) [*coordinate hyperplane*]{} of $\re^{n+1}$, is a local sub-solution of an ODE of type (\[Quad1\]) for some $m\ge n$ (with different coefficients $v_{i,l}$).

We can similarly say that (\[S1\]) is $\sigma\pi$-[*reducible*]{} (resp: under a map $\Phi$), if any solution $(x,t)$ of it, is a sub-solution (resp: a sub-solution of the type $(\Phi(x,t),t)$) of the solution of some $\sigma\pi$-ODE (generally of a larger dimension). Obviously, the class of ODEs (\[S1\]) that are exactly quadratizable includes all those $\sigma\pi$-reducible, and as shown in [@Ca15] there are indeed [*many*]{} of such functions $f$, for instance: [*(i)*]{} all finite compositions of $\sigma\pi$-functions, [*(ii)*]{} for $x\in \re$, the functions $e^x,\ln x$, and all the trigonometric functions $\sin x,\cos x, \tan x$, [*(iii)*]{} all finite compositions of functions from (i) and (ii). The problem of better characterize the [*largest*]{} set of functions that are densely quadratizable (or, equivalently: $\sigma\pi$-reducible) in terms of the more usual concepts of [*algebraic*]{} or [*differentially algebraic*]{} function, was therein only sketched, and we’ll not be concerned here either with this important issue, but save it for another occasion. It should be stressed that, on account that the class of $\sigma\pi$-ODEs is by itself quite large, the class of all quadratizable ODEs, despite based on [@Ca15] has still vaguely defined outlines, appears to be [*huge*]{} as it is, because it includes all ODEs that one encounter in practice, typically defined through finite compositions of elementary, algebraic or transcendental functions of the type we have seen just above. With that being said, the property for an ODE of being exactly quadratizable appears not ensuing given some (substantially) [*restrictive*]{} condition, and in fact to be a somewhat universal property of ordinary differential equations [*of the type*]{} (\[S1\]), i.e. [*explicitable*]{} with respect to the first order derivative.

The aim of the present paper is, first, to give an answer to some issues that remained unsolved in [@Ca15]. The first of these issues is about whether an analytic $\sigma\pi$-ODE is [*globally*]{} quadratizable (and not just [*densely*]{} quadratizable). In other words, this amounts to check whether solutions $(x,t)\not\in {\rm Dom}\{\Phi\}$ are still local sub-solution of an ODE of type (\[Quad1\]). In the first part of the paper we’ll prove a different, and more general (and even shorter) proof of the basic theorem (Theorem 2.5 in [@Ca15]) on exact quadratization, showing that, any analytic $\sigma\pi$-ODE is indeed globally exactly quadratizable, and in particular all the solutions $(x,t)\not\in {\rm Dom}\{\Phi\}$ can be calculated as well from the solutions of a quadratic systems, different and smaller dimensional, but still of the same type (\[Quad1\]).

Another issue that remained open in [@Ca15] is about whether the exact quadratization can be carried over to [*implicitly defined*]{} ODEs, of the type $f(\dot x,x,t)=0$ with both $x$ and $f$ in $\re^n$. For $f$ continuously differentiable with respect to all of its arguments, and having full rank ($n$) with respect to the $\dot x$ indeterminates, the Dini’s theorem assures that, at least locally, one can write $\dot x=\phi(x,t)$ for some continuously differentiable $\phi$. We can then talk about a solution $(x,t)$ defined in some domain, in the same sense defined earlier for explicit ODE, carrying it over to implicit ODE as well. When $f$ is an analytic $\sigma\pi$-function, the local inverse $\phi$ is still not in general a $\sigma\pi$-function. Nevertheless, we show in §\[sc:imp\] that the local explicit ODE $\dot x=\phi(x,t)$ is $\sigma\pi$-reducible, and therefore quadratizable. We perform in this way the extension of the exact quadratization theory, at least up to an extent that includes the basic $\sigma\pi$ case.

Even though it relies on the results of [@Ca15], the present paper is self-contained, and in fact the definition given before for [*exact quadratization*]{} is new, and, in our view, is more general and suitable for the present purposes[^5]. In §II we review the basic results on exact quadratization, by focusing $\sigma\pi$-ODEs, and introduce some new definitions, among which those of [*singular, regular*]{} and [*critical*]{} solution. We do not consider more general classes of ODEs, for which we still refer readers to the results of [@Ca15]. As a matter of fact, the class of ODEs $\sigma\pi$-reducible considered in [@Ca15] are all reduced to a $\sigma\pi$-ODE through maps $\Phi$ that – unlike the quadratizing map (\[ChV\]) , used for $\sigma\pi$-ODEs – have all ${\rm Dom}\{\Phi\}\subset \cal A$ ($\cal A$ is the [*initial*]{} domain, of the ODE to be put, in this case, in $\sigma\pi$-form). In other words, all $\sigma\pi$-reducible ODEs considered in [@Ca15] are in fact [*globally*]{} reduced to a $\sigma\pi$-ODE.

After having shown the globality of the exact quadratization, and its extension to implicit $\sigma\pi$-ODEs, which represent the completion of the partial results of [@Ca15], in the final part of this work we come back to the question ([*q*]{}) we posed at the beginning, and show that the answer is positive. This is the central result of the paper, showing the importance of the exact quadratization theory, since a joint application of the latter with it allows to write the solution of any analytic $\sigma\pi$-reducible ODE in its own signature analytic form, a Taylor expansion, through an explicit recursive calculation of the coefficients from the coefficients of a quadratization of the original equation. A general formula is provided as well for a lower bound of the convergence ray of the solution Taylor expansion.

In the last chapter, we sketch a possible development of the present research as well as show the interplay of the results here presented with a differential version of the concept of [*uniformization*]{} of the 22rd Hilbert’s problem.

A review, and Further Results, on the Exact Quadratization for $\sigma\pi$-ODEs
===============================================================================

A foreword on terminology: systems, sub-systems, and sub-solutions {#ssc:ter}
------------------------------------------------------------------

Keeping all definitions and symbols given as yet, hereinafter we talk about a [*system defined by*]{} an ODE of the type (\[S1\]) (hereinafter: ’the ODE’) as a certain subset (and in particular the whole) of his solutions set. If $\cal A$ is the maximal domain of the ODE, a system is given as soon as a subset $\cal S\subset \cal A$ is given, so that we talk about a [*system generated by $\cal S$*]{} (associated to the ODE, a specification that often will be omitted) meaning this the set of all solutions $(x,t)\in \cal S$. We name the subset $\cal S$ [*a generating set*]{} of the system.

We say that $\cal S\subset \cal A$ is an [*invariance set*]{} for the ODE, if for any $(x,t)\in\cal S$, $(x(s),s)\in \cal S$, $\forall s\in {\bf I}_{x,t}$. An invariance set $\cal S$ will be also named a [*domain*]{} (of the ODE). If $\cal S$ is not an invariance set, we can always define the set $\cal S^*\supset \cal S$, said [*the completion*]{} of $\cal S$, collecting all points of solutions graph passing through a point of $\cal S$. Then, by construction, $\cal S^*$ is an invariance set. By reason of that, we often identify a system generated by a domain $\cal S$, with the domain $\cal S$ itself, and talk about the system as: the [*system $\cal S$*]{}. The largest invariant set is $\cal A$, and the ODE itself is a system (the system $\cal A$). Whether there isn’t any pre-defined $\cal S$, then the system associated to the ODE is understood having domain $\cal A$. We also call the system $\cal A$ the [*total*]{} system. Note that a generating set, or a domain, $\cal S\subset \cal A$ is not necessarily open. In the following we usually write the ODE component-wise: x\_i=f\_i(x,t),\[S1c\] understanding that $i\in \cal I=\{1,\ldots, n\}$, where $\cal I$ is named the [*set of indices*]{} (of the ODE). Let $\cal S\subset \cal A$ be a generating set, and $\cal I^*\subset \cal I$ be a subset of indices. We define the [*sub-system generated by $(\cal S,\cal I^*)$*]{} — of the system $(\cal A,\cal I)$, which possibly is understood — to be the set of all functions $\{x_j(\cdot), j\in \cal I^*\}$ which – following the definition given in §I – are local sub-solutions of solutions in $\cal S$ with respect to $\cal I^*$ (we’ll also refer to them as [*solutions of*]{} $(\cal S,\cal I^*)$). If $\cal S$ is a domain, than the sub-system generated by $(\cal S,\cal I^*)$ will be simply said: the sub-system $(\cal S,\cal I^*)$. Of course, for any given $\cal S_1\subset \cal S_2 \subset \cal A$, system $\cal S_1$ is a subsystem of system $\cal S_2$ and so on, if one understands that $\cal S_k$, $k=1,2$, is actually $(\cal S_k, \cal I)$, and $\cal A$ is actually $(\cal A,\cal I)$, in accordance with our definition. A sub-system $(\cal S_1,\cal I)$ of $(\cal S_2,\cal I)$ will be also said to be a [*full*]{} sub-system, whereas if $\cal I^*$ is a proper sub-set of $\cal I$, $(\cal S_1,\cal I^*)$ will be said a [*index-proper*]{} sub-system — of system $(\cal S_2, \cal I)$, which is also named [*the hull system*]{} (of system $(\cal S_1,\cal I^*)$).

Let $\cal S_1,\ldots,\cal S_p\subset \cal S$, and $\cal I_1,\ldots,\cal I_q\subset \cal I$ be disjoint partitions of the sets $\cal S$ and $\cal I$ respectively, then the collection $\{(\cal S_l,\cal I_m)\}_{l,m}$ will be said: a [*partition*]{}, and the sub-systems generated by $(\cal S_l,\cal I_m)$ will be said [*parts*]{}, of the system generated by $(\cal S,\cal I)$. If $\cal S$ is a domain, we say that the partition $\{(\cal S_l,\cal I_m)\}_{l,m}$, is [*disjoint*]{}, if the sub-systems generated by $({\cal S}_{l_1},{\cal I}_{m_1})$, $({\cal S}_{l_2},{\cal I}_{m_2})$, are disjoint for $l_1\ne l_2$, and $m_1=m_2$, where ’disjoint sub-systems’ means that they haven’t solutions in common. Note that a system partition in [*full*]{} sub-systems (namely a [*full system-partition*]{}) is necessarily disjoint if all the parts are complete (i.e. are defined on sub-domains), but this not necessarily occurs if the parts are index-proper sub-systems. Also note that, since every system with generating set $\cal S_l$ has a domain (the completion of $\cal S_l$: $\cal S_l^*$) it follows that from any full system-partition one can derive a disjoint partition, namely $\{(\cal S_l^*, \cal I)\}_{l}$, where $l=1,\ldots,p^*\le p$ (obviously completing all parts cannot increase the number of parts).

Autonomous sub-systems {#ssc:aut}
----------------------

A particular kind of sub-system of a given system $(\cal A,\cal I)$ occurs as the sub-system, say $(\cal S,\cal I^*)$, addresses (through $\cal I^*$) equations of the ODE that includes indeterminates $x_j$ with $j\in \cal I^*$ only. We name these sub-systems [*autonomous sub-systems*]{}. In such a case, the sub-system solutions are actually independent of the components $x_j$, for $j\not\in\cal I^*$, of the point $(x,t)\in \cal S$. This means that, denoting: $$\cal T=\{t\in \re \quad|\quad \exists x\in \re^n: (x,t)\in \cal A \}$$ $\cal S$ decomposes as follows $$\cal S=\bigcup_{t\in \cal T}\left({\cal P}_t^*\times {\cal Q}_t\times \{t\}\right),$$ where, if $\iota$ is the cardinality of $\cal I^*$, for any $t\in \cal T$ we have $\cal P_t^*\subset \re^{\iota}$, and ${\cal Q}_t\subset \re^{n-\iota}$. Thus, provided that the canonical injection $\re^{\iota+1}\to \re^{n+1}$ is applied, an autonomous sub-system $(\cal S,\cal I^*)$ (in $\re^{n+1}$) can be identified with the system $(\cal P^*, \cal I^*)$ (in $\re^{\iota+1}$), with $\cal P^*=\cup_{t\in \cal T}(\cal P^*_t\times \{t\})$.

An important case of autonomous subsystem, that we encounter in the following, is that given by all sub-solutions zero in ${\cal I}^*$, i.e. by solutions of a given system $\cal S$, such that $x_i\equiv 0$, $\forall i\in {\cal I}^*$. We call them sub-system [*projected on ${\cal I}\setminus {\cal I}^*$*]{} which – provided is non-void – is described by the scalar ODEs addressed by $\cal I\setminus \cal I^*$, modified by setting $x_i=0$ for any $i\in \cal I^*$. Let us denote by $S^0_i$, $i=1,\ldots,n$, the hyperplane $x_i=0$ in $\re^{n+1}$, and by $S^0$ their union: &&S\^0\_i={(x,t) \^[n+1]{}/ x\_i=0},\[S0i\]\
&&S\^0=\_i S\^0\_i.\[dfS0\] The sub-system projected on $\cal I\setminus \cal I^*$ is the index-proper sub-system $(\cal S^*,\cal I\setminus \cal I^*)$, with \^\*=S ( \_[jI\^\*]{}S\^0\_j),\[barS\] of the [*projected total system*]{} $(\cal A^*,\cal I)$, where $\cal A^*$ is defined by the same identity (\[barS\]) but the symbol $\cal A$ replaces the symbol $\cal S$. Note that $(\cal A^*,\cal I)$ is an index-proper subsystem of $(\cal A,\cal I)$ in $\re^{n+1}$. Also, we can consider $(\cal S^*,\cal I\setminus \cal I^*)$ as a [*full*]{} sub-system of the [*reduced*]{} system $(\cal A^*,\cal I\setminus \cal I^*)$ in $\re^{\iota+1}$, as soon as $\cal S^*$ and $\cal A^*$, which are defined as the traces on the hyperplanes $x_j=0$ $j\in 
{\cal I}^*$ of $\cal S$ and $\cal A$ respectively, are identified with subsets of $\re^{n-\iota+1}$ through the canonical injection $\re^{n-\iota+1}\to \re^{n+1}$.

Decomposition of analytic systems in regular and singular parts {#ssc:sinreg}
---------------------------------------------------------------

\[df:crs\] [*For a given system $\cal S$, we name ’critical set’ the set =[S]{}S\^0.\[Sst\] The points of ${\bf C}$ are accordingly named: critical points. Let us consider the smallest subset of indices $i_1,\ldots,i_L\in \{1,\ldots,n\}$ such that \_[j=1]{}\^L S\^0\_[i\_j]{}.\[senind\] Then, $i_1,\ldots,i_L$ will be called: criticality indices (of the system). The set of all criticality indices is denoted \^\*={i\_1,…,i\_L}. \[ist\]* ]{}

Stated in words, the critical set is the set of point of the ODE domain which are also points of some coordinate hyperplane. Along with critical points and criticality indices, we also give the following definition.\

\[df:reg\] [*A solution, $(x,t)\in \cal S$ is said to be ’regular’ if the set ${\bf T}_x\subset {\cal T}$ defined as \_x={sT/x(s)},\[CT\] has not limit points in $\cal T$.* ]{}

We suppose now that the ODE is analytic, i.e. the function $f$ is real-analytic on the domain $\cal A$. Thus, the sub-domain $\cal S\subset \cal A$ determines the [*analytic system*]{} $\cal S$.\

\[th:reg\] [*A solution, $(x,t)$, of an analytic system $\cal S$ is regular if and only if his curve passes through a point in S\_r=SS\^0.\[Dr\]* ]{}

[*Proof.*]{} If ${\bf T}_x$ has not limit points in $\cal T$, then, denoted ${\bf I}_{x,t}$, as usual, the domain of the solution $(x,t)$, we can take an interval $(\alpha,\beta)\subset {\bf I}_{x,t}\subset \cal T$, with $(\alpha,\beta)\cap {\bf T}_x=\emptyset$, thus for any $t\in(\alpha,\beta)$, $x(t)\not\in {\bf C}$ and, hence, $x(t)\in \cal S\setminus S^0$. Vice-versa, consider a solution passing through $x\in \cal S\setminus S^0$, and let $x(\cdot)$ be his curve . We have x\_i0,i=1,…,n.\[xne0\] Since the ODE is analytic the curve $x(\cdot)$ is analytic on his domain ${\bf I}_{x,t}$, and ${\bf T}_x$ is the set of zeroes of the maps $x_i(\cdot)$ (we can consider just the components $i\in {\cal I}^*$). All of these maps are not identically zero, by reason of their continuity and by (\[xne0\]). The thesis follows by well known properties of the analytic functions. $\bullet$\
Theorem \[th:reg\] implies that any solution of an analytic system either is regular – hence, it possibly crosses the coordinates hyperplanes only at discrete points forming a [*non dense*]{} sub-set in ${\bf C}$ – or it lies entirely in ${\bf C}$. This motivates the following definition.\

\[df:sin\] [*A solution $(x,t)$ of an analytic ODE, is said to be ’singular’ (resp: ’singular in $i$’) if $x_i\equiv 0$, for some $i\in \cal I^*$ (resp: for $i\in \cal I^*$). The set of all singular (resp: non singular) solutions is named ’singular part’ (resp: ’regular part’) of the $\sigma\pi$-ODE.*]{}

0.8cm

Let us denote by $\sigma$ (resp: $\sigma_i$) the set of all singular solutions (resp: in $i$). Then $\sigma$ is the singular part of the ODE. Accordingly, we also call $\sigma_i$ the [*singular part in $i$*]{}, or [*$i$th singular part*]{}. We have =\_[jI\^\*]{} \_j.\[sig\]

Note that, in (\[sig\]), there well may be $\sigma_j=\emptyset$ (there are not singular solutions in $j$). Now, recall that the sub-system projected on $\{i\}$ is the set of all solutions of the ODE obtained by removing the $i$-th equation and by setting $x_i=0$ in all the others, thus it is nothing else than the singular part in $i$. That said, and accounting of (\[barS\]), the following proposition ensues.\

\[pr:sin\] [ *The singular part in $i\in \cal I^*$ of an analytic system $\cal S$ is either empty, or is the sub-system projected on $\{i\}$. The domain of the $i$th singular part is \_i=S S\^0\_i.\[singi\]* ]{}

0.8cm

Let us denote by $\cal I_s\subset \cal I$ the set of all indeces $j\in \cal I^*$ such that $\sigma_j\ne \emptyset$. Denote by ${\bf S}$ the domain of $\sigma$, by (\[sig\]) and (\[singi\]) we have =S (\_[iI\_s]{} S\^0\_i).\[singdom\]

Since a solution is either regular or singular the following proposition holds.\

\[pr:sreg\] [*Let be given an analytic system $\cal S$. Define ${\bf R}=\cal S\setminus {\bf S}$. Then $\{{\bf S}, {\bf R}\}$ is a disjoint partition of the system $\cal S$, where ${\bf S}$ (resp: ${\bf R}$) is the singular (resp: regular) part.*]{}

0.8cm

Thus, by Proposition \[pr:sreg\] the sub-domain ${\bf R}\subset \cal S$ (the domain of the regular part, or the regular part [*tout court*]{}, by our usual identification point-solution) is given, accounting of (\[singi\]), by =S (\_[iI\_s]{} S\^0\_i) =S\_r (\_[iI\^\*I\_s]{} S\^0\_i), \[regdom\] where the last identity is obtained by exploiting (\[Dr\]). Note that the regular part has domain ${\bf R}$, but it can be described as well as a system [*generated by*]{} $\cal S_r$. As a matter of fact, $\cal S_r$ is not (in general) a domain, and ${\bf R}$ is the completion of $\cal S_r$. The generating set $\cal S_r$ is a domain if and only if $\cal I^*=\cal I_s$, i.e. all criticality indices are singularity indices, in which case we have $\cal S_r={\bf R}$. The latter case occurs as soon as all solutions $(x,t)\in \cal S_r$ have a curve that never crosses a coordinate hyperplane, whereas any solution $(x,t)\in S^0_i$, for some $i$, has a curve all lying in $S^0_i$. We denote by $\cal I_{ns}$ the set of all [*non-singular criticality indices*]{}: I\_[ns]{}=I\^\*I\_s.\[nsci\] If $\cal I_{ns}\ne \emptyset$, there exist regular solutions $(x,t)\in S^0_i$, for some $i\in \cal I_{ns}$, and we call them [*critical solutions*]{}.

Indices convention {#ssc:ic}
------------------

Hereafter we adopt the following convention: [*a vector*]{} ${\bf v}$ in $\re^\alpha$, $\alpha$ a finite integer, is always a [*column*]{} vector, and ${\bf v}'$ his transpose. Another convention we make, which will help a lot in reducing the length of the most used mathematical expressions, is the following [*indices convention*]{}: if a multi-indexed symbol, say $\xi_{i_1,\ldots,i_n}$, has been defined at some point of the paper as a [*scalar*]{} quantity, then thereinafter the same symbol with the last index removed: $\xi_{i_1,\ldots,i_{n-1}}$ shall denote the corresponding [*column vector*]{} gathering the scalars for all values of $i_n$. Also, if the convention has been applied once, then it can be applied to the aggregated vector again, so that $\xi_{i_1,\ldots,i_{n-2}}$ will be the stack of all column vectors for all values of $i_{n-1}$, and so on.

Writing $\sigma\pi$-ODEs {#ssc:wr}
------------------------

By using this indices convention in (\[SP-form1\]), we can rewrite a $\sigma\pi$-ODE in the following short form: x\_i=v\_i’X\_i,\[SP\] where $i=1,\ldots, n$. It should be noted that the $i$th equation of a $\sigma\pi$-ODE may be null [*in many ways*]{}, the extreme possibilities being: if $\nu_i=0$ (which means: there are not monomials), or the vector $v_i=0$, but it cannot be $X_i=0$ though, which is never verified. Indeed, $X_{i,l}$ is a [*function*]{} (of $x$), so with $X_{i,l}=0$ we shall mean that $X_{i,l}$ is identically zero (on its maximal domain), which is impossible by definition (there isn’t a choice of exponents $p^l_{i,j}$ such that $X_{i,l}\equiv 0$ on an open set of $\re^n$). So, which and how many monomials a $\sigma\pi$-ODE consists of, is a choice [*ad libitum*]{} rather than a material fact. Anyway, in the following, as we say that a $\sigma\pi$-ODE is [*given*]{}, we shall mean that all the monomial numbers $\nu_i$’s, coefficients $v_{i,l}$’s and monomials $X_{i,l}$’s have been fixed.

Multipartite systems {#ssc:MPS}
--------------------

We say that a system in $n$ indeterminates is [*multipartite*]{}, and, precisely: [*$m$-partite*]{}, with an integer $m$, if it is double-indexed as follows: x\_[i,l]{}=f\_[i,l]{}(x,t),\[mps\] where $i=1,\ldots,m$, $l=1,\ldots,\nu_i$ (the index $l$ depends in general of $i$), and $\nu_i$ are integers such that $\nu_1+\ldots,\nu_m=n$. Recall that, by the indices convention, $x'=[x_1',\ldots, x_m']$ where $x_i'=[x_{i,1},\ldots,x_{i,\nu_i}]$. In accordance with §\[ssc:ter\] we name [*parts*]{}, the $m$ sub-systems composing (\[mps\]). In particular, the [*$i$th part*]{} will denote the sub-system — written in [*vector*]{} form, by applying the indices convention: x\_i=f\_i(x,t).\[mpsv\] We say that a system in $n^2$ indeterminates is [*regularly partitioned*]{}, if it is $n$-partite with all parts $n$-dimensional.

The domain of a $\sigma\pi$-system
----------------------------------

For a $\sigma\pi$-system the system domain is ${\cal S}={\cal D}\times {\cal T}$, where $\cal T$ is the common domain of the coefficients $v_{i,l}$, and $\cal D$ is the maximal open set on which all the functions $X_{i,l}$, given in (\[SP-form1\]), are well defined and $C^\infty$. In the following, with some abuse of terminology, we sometimes refer to $\cal D$ (instead of $\cal S$) as the system domain.

Every $\sigma\pi$-system is analytic whereby is $C^\infty$ (and, in particular, locally Lipschitz therein). Consider the following (open) sets in $\re^n$, (where $j=1,\ldots,n$): \^+\_j={x\^n| x\_j&gt;0};\^-\_j={x\^n| x\_j&lt;0};\[dfSpm\] Let us see what is the structure of any possible $\cal D$, for any given set of exponents $p_{i,j}^l$. Let us denote by $ {\cal D}_{i,j}^{l}$ the maximal set where the function $\eta_{i,j}^{l}(x)=x_j^{p_{i,j}^{l}}$ is well defined. As directly entailed by well known properties of real exponents, there are only four cases:\
[*i)*]{} ${\cal D}_{i,j}^{l}=\re^n$: for any $p_{i,j}^{l}\ge 0$ rational, say $p_{i,j}^{l}=n_1/n_2$, irreducible, and $n_2$ odd;\
[*ii)*]{} ${\cal D}_{i,j}^{l}=\overline {\bf S}^+_j$: for any $p_{i,j}^{l}>0$ irrational, or rational as in [*i)*]{} but $n_2$ is even;\
[*iii)*]{} ${\cal D}_{i,j}^{l}={\bf S}^+_j$: as in [*ii)*]{} but $p_{i,j}^{l}<0$;\
[*iv)*]{} ${\cal D}_{i,j}^{l}={\bf S}^+_j\bigcup {\bf S}^-_j$: as in [*i)*]{} but $p_{i,j}^{l}<0$;\
In [@Ca15], it is shown that, by taking all intersections of the interiors of the above sets we obtain $\cal D$, and thus we can state the following proposition.\

\[pr:SPdom\] [*The domain ${\cal D}$ of a $\sigma\pi$-system of order $n$, is always an intersection of a choice of $m\le n^2+1$ sets among the following $n^2+1$ sets: $\re^n, {\bf S}^+_j, {\bf S}^+_{j'}\bigcup {\bf S}^-_{j'}$, $1\le j,j'\le n$.*]{}

0.8cm

From Proposition \[pr:SPdom\] we can see that the domain $\cal D$ of a $\sigma\pi$-system is always an [*open positive macro-orthant*]{}, namely ${\bf O}_{i_1,\ldots,i_L}$, of $\re^n$, which is defined, for a choice $L, i_1,\ldots, i_L\in\{1,\ldots,n\}$, as \_[i\_1,…,i\_L]{}={x\^n: x\_[i\_1]{} &gt; 0,…, x\_[i\_L]{} &gt; 0},\[MO\] to which [*possibly*]{} the points of some sub-family of the family of the [*coordinate hyperplanes*]{} $x_j=0$ for $j\not\in \{i_1,\ldots,i_L\}$ are drown away. Note that the [*positive orthant*]{} of $\re^n$ is the macro-orthant ${\bf O}_{1,\ldots,n}$, whereas $\re^n={\bf O}$ (i.e. no choice of the number $L$). Notice that a $\sigma\pi$-function has always a non-empty domain and all $\sigma\pi$-functions share at least the positive orthant in their domains. As a consequence we have the following proposition.\

\[pr:O\] [*Any $\sigma\pi$-ODE has a solution $(x,t)$ for all $x\in {\bf O}_{1,\ldots,n}$ and $t\in \cal T$*]{}

0.8cm

### Example {#ExdomSP}

Let us consider the following $\sigma\pi$-ODE with $n=3$, $\nu_1=2$, $\nu_2=1$, $\nu_3=4$ :

[l]{} x\_1=x\_2\^[13]{}x\_3+[2(2t)]{}x\_2x\_1\^[-[15]{}]{}\
x\_2=6x\_1x\_2\^[5]{}x\_3\
x\_3=3x\_1\^[-8]{}x\_2+4+x\_3\^[12]{}-[32]{}x\_1x\_3

\[eq-exdomSP\] We have $X_{11}=x_2^{1\over 3}x_3$, $X_{12}=x_2x_1^{-{1\over 5}}$, $X_{21}=x_1x_2^{5}x_3$, $X_{31}=x_1^{-8}x_2$, $X_{32}=1$, $X_{33}=x_3^{1\over 2}$, $X_{34}=x_1x_3$, and the coefficients $v_{i,j}$, $i=1,2,3$, $j=1,2,3,4$, are read out at once from (\[eq-exdomSP\]). The domain of the ODE (\[eq-exdomSP\]) is $\cal S=\cal D\times \cal T$, with $\cal D$ given, according to the four cases [*i)-iv)*]{} (cf [*supra*]{}) : $$\cal D=\bigcap_{i=1}^3\bigcap_{l=1}^{\nu_i} {\cal D}^l_{i,j} = \re^3\cap {\cal D}^1_{12} \cap {\cal D}^1_{31} \cap {\cal D}^3_{33}
=\{x\in\re^3\quad | x_1\ne 0, x_3>0\}$$ (recall that $\cal D^l_{i,j}$ is the domain of $x_j^{p^l_{i,j}}$, i.e. the $j$th term of the $l$th monomial in the $i$th equation, and the domains different from $\re^3$ are those of $x_1^{-{1\over 5}}$ in $X_{12}$, $x_1^{-8}$ in $X_{31}$, and $x_3^{1\over 2}$ in $X_{33}$). As for $\cal T$, it is the domain of $v_{12}=2/\cos(2\pi t)$ only: $$\cal T=\bigcup_{k=-\infty}^{+\infty} I_k,\qquad I_k= \{t\in\re \quad | \quad  -\pi+2k\pi < t < \pi + 2k\pi\}$$ The full system associated to ODE (\[eq-exdomSP\]), i.e. the set of all solutions, is the union of a countable family of sub-systems $\{(\cal D\times I_k, \{1,2,3\})\quad | k:{\rm integer} \}$ associated to the same ODE. In terms of the set defined in (\[dfSpm\]) and (\[S0i\]) we have $$\cal D={\bf S}^+_3\setminus S^0_1 \subset \re^3$$ where the sets (\[S0i\]) are considered in $\re^n$ (i.e. $\re^3$), neglecting the time coordinate. Thus, $\cal D$ is the macro-orthant of $\re^3$ $x_3>0$, with the points of the hypher-plane $x_1=0$ removed. Clearly the sub-systems $\cal D\times I_k$ and $\cal D\times I_s$, for $k\ne s$, just differ by translation along the time coordinate. The set ${\bf C}$ defined in (\[senind\]) (identified with its time-neglecting projection in $\re^3$) is the intersection of $\cal D$ with the coordinate hypherplanes $S^0$ thus: $${\bf C}= S^0_2\setminus \{x\in S^0_2\qquad | \quad x_1= 0, \quad \!\!{\rm and}\quad \!\!x_3\le 0\}$$ and the set of criticality indices is $\cal I^*=\{2\}$. We see that $x_2\equiv 0$ is a sub-solution in the second ODE of (\[eq-exdomSP\]), thus the singular part, ${\bf S}$ is nonempty and equal to the singular part in $\{2\}$, which in turn implies that ${\bf S}={\bf C}$, [^6] and ${\bf R}={\cal D}\setminus {\bf S}$.

Finally, and coming back to the notation in $\re^3\times I_k$, for any integer $k$ we have that the sub-system $(\cal D\times I_k, \{1,2,3\})$ can be decomposed in the two sub-systems $({\bf R}\times I_k, \{1,2,3\})$, and $({\bf S}\times I_k, \{1,2,3\})$ (the regular and singular parts). The singular part is the sub-domain given by the half-plane $S^0_2$ for $x_3>0$ without the line $x_1=0$, and all solutions of it are given by the projected (on $\{1,3\}$) sub-system: $({\bf S}\times I_k,\{1,3\})$, which is associated to the ODE obtained from (\[eq-exdomSP\]) by setting $x_2=0$ and removing the second equation:

[l]{} x\_1=0\
x\_3=4+x\_3\^[12]{}-[32]{}x\_1x\_3

\[ex-singpart\] The projected system associated to (\[ex-singpart\]) is regular (it has in fact an empty singular part) [^7] and, hence, ${\bf R}\times I_k$ and ${\bf S}\times I_k$ collect all solutions of the $k$th sub-system $\cal D\times I_k$.

In conclusion, all solutions $(x,t)$, for $t\in I_k$, and $x_2=0$, $x_1\ne 0$, and $x_3>0$ are singular: they all lie in the half-plane ${\bf S}$. All solutions $(x,t)$ with $x_1,x_2\ne 0$ and $x_3>0$ are regular, and lie entirely in ${\bf R}$. Notice that, the regularity domain is composed by the four open disjoint ortants of $\re^3$ with $x_3>0$: a solution passing through an $x\in {\bf R}$ remains in the same ortant and never crosses the plane $S^0_2$. This is because $S^0_2$ includes the singular part, and all solutions passing through a point of ${\bf S}$ have to lie entirely in ${\bf S}$.\

### Example {#example}

Let us modify ODE (\[eq-exdomSP\]) by adding the monomial $x_1$ to the second equation. Clearly the system domain does not change, and the set of criticality indices is still $\{2\}$. What changes is that the new system has no more a singular part, i.e. it is regular. In this case solutions $(x,t)$ with $x_2=0$ – i.e. solutions crossing the hyper-plane $S^0_2$ – exist for any $x\in \cal D\cap S^0_2$, but none of them lies entirely in $S^0_2$: for an enough small $\epsilon>0$ the solution $(x,t)\in \cal D\cap S^0_2 \times I_k$ satisfies $x(t\pm \epsilon)\not\in S^0_2$, i.e. [*escapes*]{} from $S^0_2$, and we can identify $(x,t)$ with $(x(t+\epsilon), t+\epsilon)$ (or $(x(t-\epsilon), t-\epsilon)$).\

Exact quadratization of $\sigma\pi$-ODE {#projODE}
---------------------------------------

The theorem we present now, is a version of the first statement of Theorem 2.5 in [@Ca15], in the formalism of the present paper, and in accordance with the definition of exact quadratization we have given in the Introduction. Let be given a $\sigma\pi$-ODE in $n$ indeterminates with $\nu_i$ monomials for $i=1,\ldots,n$, and consider the map $\Phi: {\rm Dom}\{\Phi\}\subset \re^{n+1}\to \re^d$ (recall $d=\nu_1+\ldots+\nu_n$): \_[i,l]{}(x,t)=x\_i\^[-1]{} X\_[i,l]{}.\[Ph\](actually, constant in the second argument). Recall that $l=1,\ldots,\nu_i$, and whether $\nu_i=0$ the corresponding equation is null ($\dot x_i=0$) and there is no $i$-labeled entry $\Phi_{i,l}$ in the aggregated vector $\Phi$. Also, recall that the $i$th equation may be zero with $\nu_i\ge 1$, $v_i=0$ and $X_{i,l}$ be [*any*]{} monomial for $l=1,\ldots,\nu_i$ (cf. §\[ssc:wr\]). The [*zero system*]{}, in particular, i.e. the ODE $\dot x_i=0$, for $i=1,\ldots,n$, is a $\sigma\pi$-ODE with an arbitrary domain. We conventionally assume $\re^n$ be the domain of the zero system in $n$ indeterminates.\

\[th:EQns\]

*(Exact Quadratization Theorem – non singular form) The regular part, $\cal S_r=\cal D_r\times \cal T$, of any $\sigma\pi$-ODE with domain $\cal S\subset \re^{n+1}$ is exactly quadratizable. Moreover, it undergoes the following quadratization under the map $\Phi^*: (x,t)\mapsto (x,\Phi(x,t))$, with $\Phi$ defined in (\[Ph\]): &&x\_i=(v\_i’Z\_i) x\_i \[fin\]\
&&Z\_[i,l]{}=\_[j=1]{}\^n\^l\_[i,j]{} (v’\_jZ\_j)Z\_[i,l]{} \[dri\] where \^l\_[i,j]{}=p\^l\_[i,j]{} -\_[i,j]{},\[pi\] $\delta_{i,j}$ being the Kronecker symbol: $=1$ (resp: $=0$) if $i=j$ (resp: $i\ne j$), and (\[fin\]) is (i.e.: remains) the zero equation, if $\nu_i=0$. Finally, the set of all solutions $(\Phi^*(x,t),t)=(x,\Phi(x,t),t)\in \re^{n+d+1}$, with $(x,t)\in {\rm Dom}\{\Phi\}$, of the ODE (\[fin\])–(\[dri\]) is the system (with domain) $\cal S^*\subset \re^{n+d+1}$:*

[rl]{} S\^\*&={(x,Z,t)D\_r\^d / Z=(x,t)},\
&=[D]{}\_r([S]{}\_r)T=\^\*([S]{}\_r)T.

\[sbdom\]

0.8cm

\[rm:EQ\] Of course the ODE (\[fin\])–(\[dri\]) is a quadratization according to the definition given in the Introduction, provided in the ODE (\[Quad1\]) we set $m=n+d$, replace in (\[fin\])–(\[dri\]) $Z_{i,l}$ with $x_{n+s}$, where s=\_i+l,\_i=\_[k=0]{}\^[i-1]{}\_[k]{} ,\_0=0[)]{}\[s\] and finally replace $v_{i,j}$ of (\[Quad1\]) with coefficients $v^*_{i,j}$ suitably defined as functions of the parameters $\pi, v$ in (\[fin\])–(\[dri\]).

0.8cm

[*Proof of Theorem \[th:EQns\].*]{} Let us consider a solution $(x,t)\in {\bf R}$, the domain of the regular part, in accordance with (\[regdom\]). The analyticity domain of the map $\Phi$ in (\[Ph\]) — considered as a map $\Phi:\re^{n+1}\to \re^d$ — is $\cal S_r$ defined in (\[Dr\]), and by Theorem \[th:reg\] we can always choose $(x,t)\in \cal S_r$. Therefore the corresponding curve $s\mapsto x(s)$ is defined and differentiable in an open neighborhood of $t$, and the function Z(t)=(x(t),t).\[Zt\] is defined and differentiable in the same neighborhood. That said, let us calculate the time derivatives of both sides of identity (\[Zt\]). By the chain rule, denoting by $\partial_{z}$ the partial derivative operator with respect to a real indeterminate $z$, and using (\[SP\]) we have Z\_[i,l]{}=\_[j=1]{}\^n (\_[x\_j]{}Z\_[i,l]{})x\_j+\_t Z\_[i,l]{}= \_[j=1]{}\^n(\_[x\_j]{} Z\_[i,l]{})v’\_jX\_j .\[Thchr\] Now, by using (\[pi\]), and the definition of $\Phi$ given in (\[Ph\]), we have that $Z$ satisfies, component-wise, the identity: $
Z_{i,l}=\prod_{j=1}^n x_j^{\pi^l_{i,j}},
$ therefore $$\partial_{x_j} Z_{i,l}=\pi_{i,j}^l \left(\prod_{\tau=1}^n x_\tau^{\pi^l_{i,\tau}}\right) x_j^{-1}=
\pi_{i,j}^l Z_{i,l}x_j^{-1},$$ which yields: (\_[x\_j]{} Z\_[i,l]{})v’\_jX\_j &&= \_[i,j]{}\^l Z\_[i,l]{}x\_j\^[-1]{} (v’\_jX\_j )=\_[s=1]{}\^[\_j]{} \_[i,j]{}\^l Z\_[i,l]{}x\_j\^[-1]{} v\_[j,s]{}X\_[j,s]{}\
&&=\_[s=1]{}\^[\_j]{} \_[i,j]{}\^l Z\_[i,l]{} v\_[j,s]{}Z\_[j,s]{} = \_[i,j]{}\^l (v’\_jZ\_j )Z\_[i,l]{}.Replacing the expression above in (\[Thchr\]), the result is (\[dri\]). The ODE (\[fin\]) directly comes from the $\sigma\pi$ short form (\[SP\]) on account of (\[Zt\]). Thus, not only $(x,t)$ is a local sub-solution, in the $x$ components, of the solution $\Phi^*(x,t),t)$ of (\[fin\])–(\[dri\]), but it holds (\[Zt\]) as well, at least in a neighborhood of $t$, which proves that $\Phi^*$ is locally compatible. Therefore all the solutions of (\[fin\])–(\[dri\]) of the form $(\Phi^*(x,t),t)=(x,\Phi(x,t),t)$ lie in the set $\cal S^*$ defined in (\[sbdom\]), which is then a sub-domain for the ODE (\[fin\])–(\[dri\]). This complete the proof of the theorem. $\bullet$\

\[df:GEQ\] [*We say that an analytic system $\cal S$ is ’globally (exactly) quadratizable’ if it can be partitioned by a finite collection, $\cal S_1,\ldots,\cal S_L$ of disjoint sub-systems and there exist maps $\Phi^1,\ldots,\Phi^L$, with ${\rm Dom}\{\Phi^i\}\cap \cal S_i\ne \emptyset$, such that the sub-system $\cal S_i$ is quadratizable under the map $\Phi^i$ for $i=1,\ldots,L$.* ]{}

0.8cm

\[th:GEQ\] [*(Exact Quadratization Theorem – Global form) Any analytic $\sigma\pi$-ODE is globally exactly quadratizable.* ]{}

0.8cm

[*Proof.*]{} As usual we consider a $\sigma\pi$-ODE in $n$ indeterminates, with a maximal domain $\cal S=\cal D\times \cal T$. Moreover we will write ${\bf R}(\cal S), {\bf S}(\cal S)$ for the regular and singular part respectively of a system $\cal S$. Let us build up the following sequence of systems: $\cal S_0={\bf S}(\cal S)$, $\cal S_{k+1}={\bf S}(\cal S_k)$ for $k=0,1,\ldots $, and let us suppose that $\cal S_k$ is singular (i.e. has a non empty singular part). By denoting $\cal I^{\bf s}_{k}$, the singularity indices for the system $\cal S_k$, we have $\cal I^{\bf s}_k\ne \emptyset$. Now, by definition $\cal S_{k+1}$ is $\cal S_k$ projected on $\cal I\setminus \cal I^{\bf s}_k$, thus $\cal S_{k+1}$ is a sub-set of solutions with at least one more zero component than $\cal S_k$, namely $x_i\equiv 0$ with $i\in \cal I^{\bf s}_k$. Therefore there is an $L<+\infty$, such that either $\cal S_L$ is regular (all of its solutions are regular) or it is constituted by the zero solution only $x\equiv 0$ (the ’zero system’). By denoting ${\bf R}_k={\bf R}(\cal S_k)$, for $k=0,\ldots, L$, we have that system $\cal S$ can be partitioned as $({\bf R}_0, \cal S_0)$, where $\cal S_0$ can be partitioned as $({\bf R}_1, \cal S_1)$ and so on, up to $\cal S_{L-1}$ which is partitioned as $({\bf R}_L, \cal S_L)$. Finally, $\cal S$ is partitioned by $({\bf R}_0, {\bf R}_1, \ldots, {\bf R}_{L-1}, \cal S_L)$ where either $\cal S_L=\emptyset$ or it is the zero-system, which is trivially quadratizable under the identity map[^8]. The theorem is proven by applying the non-singular form of the Exact Quadratization Theorem to the regular systems ${\bf R}_0, {\bf R}_1, \ldots, {\bf R}_{L-1}$. $\bullet$\
From the proof of Theorem \[th:GEQ\] we see that any solution non identically zero of an analytic $\sigma\pi$-ODE is a regular solution of some of its sub-systems, and in particular of the sub-system projected on the complement of a set of singularity indices. Thus, any non-trivial solution can be calculated as analytic prolongation of a sub-solution of a suitable quadratic system. For this reason, it will be without loss of generality the assumption we make hereinafter: as we talk about ’a solution’ of an analytic $\sigma\pi$-ODE, we always understand, unless explicitly specified, that it is a [*regular*]{} solution.\

Exact quadratization for implicitly defined $\sigma\pi$-systems {#sc:imp}
---------------------------------------------------------------

Let us consider the [*implicitly defined*]{} ODE: f(x,x,t)=0,\[IDO\] where $x\in \re^n$ and $f: U {\rm (open)}\subset \re^{2n+1}\to \re^n$ is a continuously differentiable function having full rank ($n$) in the $\dot x$ indeterminates, at a point $z=(x,\dot x,t)\in U$, i.e. $(\partial_{\dot x} f)(x,\dot x,t)$ has rank $n$ (with $\partial_{\dot x}$ denoting the Jacobian with respect to $\dot x$). Let $\cal U$ be the set in $U$ where the algebraic equation (\[IDO\]) is satisfied: U={(x,x, t)U| f(x,x,t)=0}.\[calU\] Then, by the Dini’s theorem (in many indeterminates), for any $(x,\dot x, t)\in \cal U$ there exist two open sets $U_1,U_2$, $(x,\dot x,t)\ni U_1\subset U$, $(x,t)\ni U_2\in \re^{n+1}$, and a continuously differentiable function $\phi:U_2 \to \re^n$, such that, denoting $y=(x,t)$: &&UU\_1={(x,(x,t),t) | (x,t)U\_2 },\[If1\]\
&&(\_y)(x,t) =-((\_[x]{} f)(x,(x,t) ,t))\^[-1]{}((\_[y]{} f)(x,(x,t) ,t)), yU\_2.\[If2\] Therefore, the [*explicit*]{} ODE x=(x,t),\[expO\] on the domain $U_2\subset \re^{n+1}$ has solutions that are local solutions of the implicit ODE (\[IDO\]). If $(\partial_{\dot x} f)(z)$ has full rank $\forall z\in U$, then all solutions of (\[IDO\]) are local solutions of some explicit ODE of the type (\[expO\]). Then we can naturally extend the definition given earlier of [*system*]{} (associate to an explicit ODE) to [*implicit systems*]{} (associated to an implicit ODE), by defining an implicit system as a [*domain*]{} $\cal S\subset \re^{n+1}$, and a solution $(x,t)\in \cal S$, as the maximal map $s\mapsto x(s)$, with $x(t)=x$, which solves the implicit ODE (\[IDO\]). From what said earlier it is $\cal S\ne \emptyset$, since for any $(x,\dot x,t)\in \cal U$ there is always an explicit ODE locally solved by $(x,t)$.

These facts carry over to the analytic case, and we have an analytic local inverse $\phi$ associated an analytic implicit system function $f$ at each point of its own domain. Now, the question we pose is the following: suppose $f$ is an analytic $\sigma\pi$-function in $2n$ indeterminates, which means it has the well known $\sigma\pi$ expression with respect to the indeterminates $x,\dot x$, for some set of analytic coefficients $v_{i,l}$, then is any associated local inverse $\phi$ an analytic $\sigma\pi$-function as well? The answer is in general negative, but the following Theorem shows that the explicit local system given by (\[expO\]) is always $\sigma\pi$-[*reducible*]{} (cf. §I) which is sufficient for our present purposes.\

\[th:imp\] [*Let be given an analytic $\sigma\pi$-ODE in the implicit form (\[IDO\]), with $f$ of full rank ($n$) with respect to the $\dot x$ indeterminates, at a point $z=(x,\dot x,t)\in \cal U$. Then, the associated explicit local representation (\[expO\]) is $\sigma\pi$-reducible.* ]{}

0.8cm

Before proving the Theorem, some comments are in order. Recall that with ’$\sigma\pi$-reducible’ we mean that the ODE (\[expO\]) has solutions that are all sub-solutions of a (larger-dimensional) $\sigma\pi$-ODE, i.e. there are additional indeterminates $x_{n+1},\ldots,x_{n+d}$, for some $d\in\na$, and analytic time-functions $v_{i,l}$ such that: x\_i= \_[l=1]{}\^[\_i]{}v\_[i,l]{}\_[j=1]{}\^[n+d]{} x\_j\^[p\^l\_[i,j]{}]{},i=1,…,n+d,\[SPes\] such that any solution $(x,t)\in \cal D\subset \re^{n+1}$ where $\cal D$ is the domain of (\[expO\]) is a local sub-solution of (\[SPes\]) in the first $n$ indeterminates. We can also say that the system $\cal D$ associated to ODE (\[expO\]) is $\sigma\pi$-reducible if it is the sub-system $(\cal D,\{1,\ldots,n\})$ associated to the ODE (\[SPes\]). Clearly, a $\sigma\pi$-reducible system undergoes a quadratization, and, hence, by Theorem 2 we can naturally build up the extension of the exact quadratization theory to implicit analytic full-rank $\sigma\pi$-systems of the type (\[IDO\]).\
Some classes of $\sigma\pi$-reducible (explicit) systems have been studied in [@Ca15], these classes of systems are defined by the class of functions the system-function belongs to. We here recall one of them, which is very basic, referring reader to [@Ca15] for more detail: (system function given by) a [*finite composition*]{} of $\sigma\pi$-functions, which in [@Ca15] is named, with some abuse of terminology, an [*algebraic*]{} function. For instance, consider the following ODE:

[l]{} x\_1= [x\_1\^2x\_2\^[-[14]{}]{} x\_1+x\_2\^2]{}\
x\_2=x\_1\^2x\_2

\[exNSP\] the first equation is not $\sigma\pi$, but the system function can be written as $\phi_1\circ \phi_2$, where $\phi_1(x_1,x_2,x_3)=x_1^2x_2^{-{1\over 4}}x_3^{-1}$ and $\phi_2(x_1,x_2)=(x_1,x_2, x_1+x_2)$, and thus, it is the composition of two $\sigma\pi$-functions. By taking the time derivative of $x_3=x_1+x_2^2$ we have $\dot x_3=\dot x_1+2x_2\dot x_2$, which is $\sigma\pi$, since $\dot x_1=\phi_1(x_1,x_2,x_3)$ is $\sigma\pi$, $\dot x_2$ is $\sigma\pi$ by the second of (\[exNSP\]), and sums and products of $\sigma\pi$-functions results in a $\sigma\pi$-function. By adjoining the equation for $x_3$ to (\[exNSP\]), and replacing $x_3$ for $x_1+x_2^2$, we get a $\sigma\pi$-ODE in the three indeterminates $x_1,x_2,x_3$. We recall also this (quite straightforward) result: let $\phi$ a scalar $\sigma\pi$-function in $x\in\re^n$, then $\partial _{x_1}f$, the partial derivative of $f$ in $x_i$, is a $\sigma\pi$-function (in $x$).\
[*Proof of Theorem \[th:imp\].*]{} By the Dini’s theorem we get the explicit ODE (\[expO\]) whose domain $U_2$ is such that every $(x,t)\in U_2$ is a local solution of (\[IDO\]) as well. Moreover $\phi$ is continuously differentiable, thus we can set $z=\phi(x,t)$ and take the time-derivative, the result is

[l]{} x=z,\
z=(\_x) x+\_t =(\_x) z+\_t

\[tran\] In (\[If2\]), note that, for $y=(x,t)$, $\partial_y\gamma = (\partial_x \gamma,  \partial_t\gamma)$ (where the symbol $\gamma$ stands for $\phi$ or $f$) thus (\[If2\]) can be rewritten as

[l]{} (\_x)(x,t) =-((\_[x]{} f)(x,(x,t) ,t))\^[-1]{}((\_[x]{} f)(x,(x,t) ,t)),\
(\_t)(x,t) =-((\_[x]{} f)(x,(x,t) ,t))\^[-1]{}((\_[t]{} f)(x,(x,t) ,t)).

\[If2\] Now, since $f(x,\dot x,t)$ is $\sigma\pi$ with respect to $x$ and $\dot x$, any component of the matrix $\partial_{\dot x}f$ or $\partial_{x}f$ is a $\sigma\pi$-function in $x,\dot x$ as well. Moreover, any component of $(\partial_xf)^{-1}$ is the ratio of the determinant of some minor and the determinant of the whole matrix, thus a ratio of polynomials – hence a $\sigma\pi$-function – of the components of $\partial_xf$. Therefore, we can set -((\_[x]{} f)(x,(x,t) ,t))\^[-1]{}((\_[x]{} f)(x,(x,t) ,t))=h(x,(x,t),t), \[conc1\] where the matrix function $h$ has finite compositions of $\sigma\pi$-functions in all entries. The function $\partial_tf$ has the usual $\sigma\pi$ expression (\[SP-form1\]) but the coefficients $\dot v_{i,j}$ replaces $v_{i,j}$, therefore is an analytic $\sigma\pi$-function, and we can set -((\_[x]{} f)(x,(x,t) ,t))\^[-1]{}((\_[t]{} f)(x,(x,t) ,t))=g(x,(x,t),t)., \[conc2\] and $g$ is still a finite composition of $\sigma\pi$-functions. By (\[If2\])–(\[conc2\]) eq.s (\[tran\]) rewrite

[l]{} x=z,\
z=h(x,z,t) z+g(x,z,t)=(x,z,t)

\[tran2\] where $\psi$ is clearly a finite composition of $\sigma\pi$-functions. In conlusion, the ODE (\[tran2\]) is $\sigma\pi$-reducible, and since the solution $(x,t)$ of the explicit ODE (\[expO\]) is a sub-solution of (\[tran2\]), we conclude that (\[expO\]) is $\sigma\pi$-reducible as well. $\bullet$\
Theorem \[th:imp\] implies in turn the following corollary.\

\[co:imp\] [*Any implicit analytic $\sigma\pi$-ODE of the type (\[IDO\]), $f$ having rank $n$ in the $\dot x$ indeterminates, for any $z=(x,\dot x,t)\in \cal U$, is globally quadratizable.* ]{}

0.8cm

It should be stressed that, unlike the case of an explicit $\sigma\pi$-ODE, any solution $(x,t)$ in the domain of an implicit analytic $\sigma\pi$-ODE is, by Corollary \[co:imp\], a local sub-solution of a quadratic system of the type (\[Quad1\]), but the latter is not in general unique, and in fact depends of the solution $(x,t)$ itself.\

\[rm:impl\] A more important remark is that, even if the function $\phi$ of the locally defined explicit ODE (\[expO\]) could be difficult to calculate, we have not to perform this calculation though in order to calculate a solution of the original implicit ODE. For, by the proof of Theorem \[th:imp\] we realize that any ODE solution is a sub-solution (in the first $n$ components) of the explicit ODE (\[tran2\]), which is defined by a vector field that can be calculated through (\[conc1\]), and (\[conc2\]), from the initial implicit vector field $f$: we can then quadratize system (\[conc1\]), (\[conc2\]) and the solution of the original implicit system is [*a fortiori*]{} a local sub-solution of the quadratized ODE, thus still recoverable by analytic extension.

0.8cm

The structure of the canonical quadratization {#sc:SCO}
---------------------------------------------

As for the ODEs (\[fin\]) and (\[dri\]), we here follow the terminology introduced in [@Ca15], but accounting as well of the new definitions given in the previous sections. Hereinafter we name the ODE (\[SP\]), with domain $\cal S$, the [*original $\sigma\pi$-system*]{}.\

\[df:quahu\][*We name the system (on $\re^{n+d}\times {\cal T}$) defined by the ODEs (\[fin\]) and (\[dri\]), the quadratic hull (of the original $\sigma\pi$-system)*]{}

0.8cm

\[df:quaim\] [*We name ’quadratic image’ or ’image system’ or ’canonical quadratization’ (of a given $\sigma\pi$-system, that is referred to as ’original system’) – the system defined by the ODEs (\[fin\]) and (\[dri\]), on the sub-domain ${\cal S}^*$ defined by (\[sbdom\]). Moreover, the two sub-systems of the quadratic image, given by the ODEs (\[fin\]) and (\[dri\]), are said the ’final stage’ (of the quadratization) and the ’Driver’ (of the original system) respectively.* ]{}[^9]

0.8cm

Note that, in accordance with §\[ssc:MPS\] the Driver is a [*$n$-partite*]{} system, whose $i$th part consists of $\nu_i$ (scalar) equations. Moreover is an [*autonomous*]{} system, i.e. the sub-system $({\cal S}^*, {\cal I}^{\bf d})$ of the image system, where $\cal I^{\rm d}\subset \cal I$ the sub-set of indices addressing the last $d$ equations of the image system, thus following §\[ssc:aut\], and accounting as well of the structure of $\cal S^*$ given in (\[sbdom\]), it can be identified with the system (in $\re^d\times \cal T$): $(\Phi(\cal D_r\times \cal T)\times \cal T, \cal I^{\rm d})$. We can omit the (obvious) symbol $\cal I^{\rm d}$, and consider the Driver with the full system: $\Phi(\cal D_r\times \cal T)\times \cal T$.

In the following, as we talk about a Driver, it shall be always entailed [*of which*]{} original ($\sigma\pi$) system, that Driver is. By reason of that, we also name [*pre-Driver*]{} the original $\sigma\pi$-system, insofar we move the focus from the original system to the associated Driver. In order to refer ODE (\[dri\]) as a sub-system of the quadratic hull we use the different term of [*Driver ODE*]{}, instead of the single word ’Driver’.

The canonic quadratization holds under the map $\Phi^*$, which we call the [*canonic quadratizing map*]{} or the [*quadratic immersion*]{}.

Non-canonical quadratizations, and some examples. {#sc:SCO}
-------------------------------------------------

According to the definition given in the previous paragraph, the Driver is a part of the canonical quadratization, as well as the final stage. Nevertheless, there is a kind of quadratization, that can be always applied to the same original system, in which the final system is redundant, and can be considered consisting of the sole Driver.\

\[df:inc-qua\] [*We say that the canonical quadratization (of a given $\sigma\pi$-ODE) is ’inclusive’, if all solutions of the original ODE are local sub-solutions of the Driver.*]{}

0.8cm

By the definition given above, a canonical quadratization is inclusive if and only if for $i=1,\ldots,n$, $\Phi_{i,l}=x_i$ for some $l=1,\ldots, \nu_i$, which is verified if for all $i$ there is an $l$ such that $X_{i,l}=x_i^2$. If for some $i$ the original equation has not the monomial $x_i^2$ one can always add the fictitious term $0\cdot x_i^2$ in the $i$th equation. The addition of a fictitious term obviously doesn’t change the solution of the ODE (\[SP-form1\]), but in the canonical quadratization it does change the dimension ($d$) of the target space for the aggregate map $\Phi$ given by (\[ChV\]), and in fact changes the map $\Phi$ itself, making redundant the quadratizing map $\Phi^*$. Indeed, we can easily verify, by applying formula (\[dri\]) that the driver equation for $Z_{i,l}=\Phi_{i,l}=x_i$ (thus, associated to the monomial $X_{i,l}=x_i^2$) is equal to the $i$th equation of the final sub-system (\[fin\]). For, let us add $0\cdot x_i^2$ at the $i$th equation as the last term, then the new $\sigma\pi$-ODE has $\nu_i^*=\nu_i+1$ monomials, and the Driver has the further variable $Z_{i,\nu_i^*}=x_i^{-1}x_i^2=x_i$. By (\[dri\]) we have $$\dot Z_{i,\nu_i^*}=\sum_{j=1}^n \pi^{\nu^*_i}_{i,j} (v_j'Z_j) Z_{i,\nu^*_i},$$ but from (\[pi\]) $\pi^{\nu^*_i}_{i,j}=\delta_{i,j}$, ($X_{i,\nu^*_i}=x_i^2=x_1^0\cdots x_{i-1}^0x_i^2x_{i+1}^0\cdots x_n^0$ from which $p_{i,j}=0$ for $i\ne j$ and $p_{ii}=2$) and thus $\dot Z_{i,\nu_i^*}=(v_i'Z_i)Z_{i,\nu_i^*}$ which is the $i$th final equation (\[fin\]), since $Z_{i,\nu_i^*}=x_i$ and accounting as well that, with respect to the canonical quadratization, $v_i$ has the further last component $0$ and $Z_i$ has the further component $Z_{i,\nu_i^*}$, which leaves the product $v_i'Z_i$ unchanged.

In conclusion, given any $\sigma\pi$-ODE, we naturally obtain an inclusive quadratization by applying the canonical quadratization to the ODE modified by adding, for $i=1,\ldots,n$, the fictitious term $0\cdot x_i^2$ in the $i$th equation, provided the monomial $x_i^2$ is not therein yet present. Clearly, the new quadratization has as many new Driver variables, say $Z_{i,\nu_i+1}=x_i$, as the original ODE has equations $i$-labeled lacking the monomial $x_i^2$. Correspondingly, a subset of the final equations in the canonical quadratization turns out to be replied among the new Driver equations. Since the natural inclusive quadratization is [*redundant*]{}, we usually omit the final stage, thus for an inclusive quadratization we assume that the quadratic hull of the original $\sigma\pi$-ODE is the Driver ODE, and identify the quadratization with the sole Driver.\

### Example {#Exnoncan}

Let us consider the following $\sigma\pi$-ODE with $n=3$, $\nu_1=2$, $\nu_2=1$, $\nu_3=2$: &&x\_1=x\_2x\_3+2x\_1\^[-1/3]{}\
&&x\_2=4x\_1x\_2\^4x\_3,\
&&x\_3=5x\_1\^[-3]{} x\_2+3.The canonical quadratization is composed by a Driver in 5 indeterminates (as many as is the number of monomials) defined as &&Z\_[11]{}=x\_2x\_3x\^[-1]{}\_1; Z\_[12]{}=x\_1\^[-4/3]{},\
&&Z\_[21]{}=x\_2\^5x\_3,\
&&Z\_[31]{}=x\_1\^[-3]{}x\_2x\_3\^[-1]{}; Z\_[32]{}=x\_3\^[-1]{}.By (\[fin\] ), the final ODE is:

[ll]{} x\_1=(Z\_[11]{}+2Z\_[12]{})x\_1 &[(]{} =(v\_1’Z\_1)x\_1 [)]{}\
x\_2=4Z\_[21]{}x\_2, &[(]{} =(v\_2’Z\_2)x\_2 [)]{}\
x\_3=(5Z\_[31]{}+3Z\_[32]{})x\_3, &[(]{} =(v\_3’Z\_3)x\_3 [)]{}.\

\[fex\] By (\[dri\]) we calculate the Driver ODE. A trick for calculating the equation of $Z_{i,l}$ is looking at the monomial which defines it: just write the linear composition of $v_j'Z_j$ using the corresponding powers of $x_j$ in the monomial as coefficients, and finally multiply by $Z_{i,l}$: &&Z\_[11]{}=(v’\_2Z\_2+v’\_3Z\_3-v\_1’Z\_1)Z\_[11]{}; Z\_[12]{}-(4/3)(v\_1’Z\_1)Z\_[12]{},\
&&Z\_[21]{}=(5v\_2’Z\_2+v\_3’Z\_3) Z\_[21]{},\
&&Z\_[31]{}=(-3v\_1’Z\_1+v’\_2Z\_2-v\_3Z\_3)Z\_[31]{}; Z\_[32]{}=-(v\_3’Z\_3)Z\_[32]{},and the $v_j'Z_j$, $j=1,2,3$, are given from (\[fex\]). The natural inclusive quadratization is obtained by simply introducing the new Driver indeterminates: $Z_{13}=x_1$, $Z_{22}=x_2$, $Z_{33}=x_3$, whose equations are directly given by (\[fex\]) upon substituting $x_1,x_2,x_3$.\
Adding a fictitious monomial is a method that was not investigated in [@Ca15]: we now show that certain side results of [@Ca15] can be derived in a straightforward way, and even generalized, by using this method. Let us consider a general monomial function of the $n$ indeterminates of a given $\sigma\pi$-ODE, namely $g(x)\in \re$: g(x)=x\^[q\_1]{}\_1x\^[q\_2]{}\_2x\^[q\_n]{}\_n;q\_1,…,q\_n,\[funcg\] By using the exact quadratization we can prove that $g$ satisfies a scalar quadratic ODE: let us add the monomial $0\cdot g$ in the $i$th equation, and apply the canonical quadratization. With the same symbols as above, we have $Z_{i,\nu_i^*}=x_i^{-1}g$, that is $g=Z_{i,\nu_i^*}x_i$, and by taking the time derivative: $$\dot g=\dot Z_{i,\nu_i^*}x_i+Z_{i,\nu_i^*}\dot x_i$$ and by the Driver and final equations we have

[rl]{} g&=(\_[ij=1]{}\^n q\_j (v\_j’Z\_j) +(q\_i-1) v\_i’Z\_i) Z\_[i,\^\*\_i]{}x\_i+ (v\_i’Z\_i)x\_iZ\_[i,\^\*\_i]{}\
&=(\_[j=1]{}\^n q\_j (v\_j’Z\_j)) g=hg

\[eqg\] Let $(x,t_0)\in \cal D_r\times \cal T$ be a solution of the $\sigma\pi$-ODE, and $I_{x,t_0}\subset \cal T$ be its maximal open interval of definition. Then, there is an interval $t_o\ni I$(open) $\subset I_{x,t_0}$ where the time-functions $$h_t= \sum_{j=1}^n q_j (v_j'(t)Z_j(t)); \qquad g_t=g(x(t))$$ are well defined (and, hence, analytic). Using (\[eqg\]), $h_t$ in $I$ is given by: g\_t=e\^[\_[t\_0]{}\^t h\_d]{} g\_[t\_0]{},\[gsolexp\] which allows the calculation of $g_t$ from a Driver solution $Z(t)$ (through $h_t$). If $g(x)=x_i$, and thus $q_j=0$ for $j\ne i$ and $q_i=1$, (\[eqg\]) turns out into the well known $i$th final equation (we come back in this way to the canonical quadratization). In the latter case, $g_t=x_i(t)$, and it is defined by definition on the whole of $I_{x,t_0}$, which is in general a larger interval than $I$. Thus, the exponential solution formula (\[gsolexp\]) gives $x_i$ only on the interval $I$ (where $Z(t)$ is defined). Outside of $I$, $x(t)$ is determined though, by analytic prolongation.

We have found out in a very simpler way the same quadratization shown in [@Ca15], where a further system: the [*medial system*]{} given by the ODE (\[eqg\]) for $g=X_{i,l}$, is ’driven’ by $Z(t)$ (by the Driver) and gives all the monomial functions $X_{i,l}(t)$, which in turn give the original solution $x(t)$ through integral action: $$x_i(t)=\int_{t_0}^t v_i'(\tau)X_i(\tau) d\tau$$ We name this kind of quadratization: a [*medial quadratization*]{}.\

### Example {#Ex2}

Let $a,b$ be two fixed positive real numbers, $t_0\in\re$, and let us consider the affine ODE: x=-a x-b.\[ex-lsys\] In our terminology, we are considering the system $(\re^2,\{1\})$ associated to ODE (\[ex-lsys\]). Since $b\ne 0$ there aren’t singular solutions, and the regularity domain is $\cal S_r=\cal D_r\times \re$ with $\cal D_r=\re^+\cup \re^-$. As a $\sigma\pi$-system, (\[ex-lsys\]) is characterized by $n=1$, $\nu_i\equiv \nu_1=2$, $v_{1,1}=-a$, $v_{1,2}=-b$, $X_{1,1}=x$, $X_{1,2}=1$, $p^1_{1,1}=1$, $p^2_{1,1}=0$. The components of the canonical Driver are $Z_{1,1}=1$, $Z_{1,2}=x^{-1}$, and by (\[dri\]) we obtain the canonical Driver equations: &&Z\_[1,1]{}=0,\[ex-dri1\]\
&&Z\_[1,2]{}=aZ\_[1,1]{}Z\_[1,2]{}+bZ\^2\_[1,2]{}.\[ex-dri2\] By (\[fin\]) we get the final equation: x=(-aZ\_[1,1]{}-bZ\_[1,2]{})x.\[ex-fin\] The Driver map is $\Phi(x,t)=(1,x^{-1})$, which has domain $\cal S_r$. The quadratic hull is the system, with domain $\re^4$, associated to the three equations (\[ex-dri1\])-(\[ex-fin\]), whereas the image system is given by the same ODEs on the domain $\re\times \Phi(\cal S_r) \times \re$.

Let us consider the solution $(x,t_0)$ of (\[ex-lsys\]), and denote $t\mapsto x(t)$ the corresponding solution curve. We can easily calculate $x(t)$: x(t)=x e\^[-a(t-t\_0)]{}-[ba]{} (1-e\^[-a(t-t\_0)]{}),\[solex\] and we have ${\bf I}_{x,t_0}=\re$. Let us define ${\bf T}_{x,t_0}=\{t\in\re: x(t)=0\} $ which is the set (\[CT\]) for the case at issue. ${\bf T}_{x,t_0}$ is a countable set, by Theorem \[th:reg\], and we can calculate it by distinguishing three cases:\
[*(i)*]{} $x\ge 0$. In this case ${\bf T}_{t_0,\bar x}=\{\bar t\}$, with |t=t\_0+[ca]{};c=-([ba|x+b]{}) &gt;0.\[ex-tmp\] and $t_0\le \bar t$.\
[*(ii)*]{} $-(b/a)<x<0$. In this case ${\bf T}_{t_0,\bar x}=\{\bar t\}$, with $\bar t$ as in (\[ex-tmp\]), but $t_0>\bar t$.\
[*(iii)*]{} $x\le-(b/a)$. In this case ${\bf T}_{0,\bar x}=\emptyset$.\
Now, consider the solution $(\Phi(x,t),t)=(1,x^{-1},t)$ of the Driver. By (\[ex-dri1\]) we have $Z_{1,1}\equiv 1$ for any $t$. The component $Z_{2,2}$ is given by the ODE (\[ex-dri2\]), which is a [*Bernoulli*]{} quadratic differential equation, and thus the general solution can be calculated by a well known method: define $W_{1,2}=Z^{-1}_{1,2}$, then (as it is easy to verify) $W_{1,2}$ satisfies the [*affine*]{} equation: W\_[1,2]{}=-aZ\_[1,1]{} W\_[1,2]{}-b.\[ex-W12\] Notice that, since $Z_{1,1}=1$, (\[ex-W12\]) is the same equation as the original system (\[ex-lsys\]). Moreover $W_{1,2}(t_0)=Z^{-1}_{1,2}(t_0)=x$. The solution is then equal to (\[solex\]), and by $Z_{1,2}=W^{-1}_{1,2}$, we get the general expression for $Z_{1,2}$: Z\_[1,2]{}(t)=[a (ax +b)e\^[-a(t-t\_0)]{}-b]{}.\[ex-SG\] The above function is (as expected) not defined for $t\in {\bf T}_{x,t_0}$. We can use formula (\[gsolexp\]) with $g(x)=x$ in order to express the solution given by the final system, i.e. ODE (\[ex-fin\]) [*driven*]{} by the Driver: x(t)=e\^[\_[t\_0]{}\^t Z\_[1]{}\^Tv\_[1]{}]{} x=e\^[-\_[t\_0]{}\^t (aZ\_[1,1]{}+bZ\_[1,2]{})]{} x, \[ex-solfinal1\] If $x\le -(b/a)$ – case [*(iii)*]{} above – (\[ex-SG\]) is always defined and non-zero on $\re$ and (\[ex-solfinal1\]) holds $\forall t\in \re$, and gives back the time-function (\[solex\]) $\forall t$. Thus. according to the Exact Quadratization Theorem, $x(t)$ is a sub-solution of the quadratic image, and for $x\le -(b/a)$ in particular is a [*global*]{} sub-solution (cf §I) defined on ${\bf I}_{x,t_0}=\re$. If $x\ge 0$ (resp: $-(b/a)<x<0$) – case [*(i)*]{} (resp: [*(ii)*]{}) – the solution $(x,t_0)$ of the original ODE is a local sub-solution defined on the sub-interval ${\bf I}=(-\infty, \bar t )$, (resp: $(\bar t, +\infty)$) with $\bar t$ given by (\[ex-tmp\]), and (\[ex-solfinal1\]) – calculated for $t\in {\bf I}$ – gives the restriction on ${\bf I}$ of the function (\[solex\]). It should be stressed that, a strict-sense local sub-solution of the original solution $(x,t_0)$ – which occurs as $x>-(b/a)$ – is given by (\[ex-solfinal1\]) only in ${\bf I}$, and it couldn’t be otherwise since the time function defined by (\[ex-solfinal1\]) is either positive or negative (which depends of the sign of $x$) everywhere in ${\bf I}$. Outside ${\bf I}$ though, it shall be still equal to $(x,t_0)$ by analytic prolongation.\

Inverse Driver {#Inv-driv}
--------------

By the Exact Quadratization Theorem (any of the two forms) a Driver solution is a [*$Z$-sub-solution*]{} of the quadratization i.e. the components from the $(n+1)$th to the $(n+d)$th of a solution lying in the domain $\cal S^*$ defined in (\[sbdom\]). Thus, we have pointwise $Z(t)=\Phi(x(t))$ for any $t\in I_{x,t}$ such that $x(t)\in \cal D_r$, and hence — by (\[Ph\]) — $Z^{-1}(t)$ is defined for the same set of $t$ as $Z(t)$ i.e. $\forall t\in I_{x,t}$ with the possible exception of a subset that hasn’t limit points. in $I_{x,t}$. Therefore, we can define W\_[i,l]{}=Z\_[i,l]{}\^[-1]{},\[dfW\] so everywhere $Z$ is defined and differentiable so does $W$ as well. By taking the time derivative, we have that $W_{i,l}$ satisfies the equation: W\_[i,l]{}=-Z\^[-2]{}\_[i,l]{}Z\_[i,l]{}=-(\_[j=1]{}\^n \_[i,j]{}\^[l]{}v’\_jZ\_j)W\_[i,l]{}\[inverse-driver\] Thus, the aggregate $W$ satisfies another $\sigma\pi$-ODE that we call [*the inverse Driver*]{} (associated to a given original system) [^10] In accordance with the terminology introduced before, the inverse Driver is the sub-system $(\cal W,\cal I_{W})$ of the ODE (\[fin\]), (\[inverse-driver\]) where $\cal I_{W}$ is the set of indices addressing (\[inverse-driver\]) and $\cal W$ is defined as

[rl]{} [W]{}&={(x,W,t)D\_r\^d/W\_[i,l]{}=Z\^[-1]{}\_[i,l]{}, Z=(x,t)},\
&=[D]{}\_r\^W([S]{}\_r)T,

\[SWst\] where $\Phi^W:{\rm Dom}\{\Phi^W\}={\rm Dom}\{\Phi\}=\cal S_r\subset \re^{n+1}\to \re^d$ has components $\Phi^W_{i,l}$ given by \^W\_[i,l]{}(x)=(\_[i,l]{}(x,t))\^[-1]{},xD\_r.\[dfPhist\] We see that the inverse Driver, as well as its [*direct*]{} counterpart, entails the calculation of the original solution, provided it feeds an [*inverse final*]{} system. Indeed, any (regular) solution $(x,t_0)\in \cal D_r\times \cal T$ is given by x\_i=(\_[l=1]{}\^[\_i]{} W\_[i,l]{}\^[-1]{} v\_[i,l]{})x\_i, x\_i(t)=e\^[\_[l=1]{}\^[\_i]{}\_[t\_0]{}\^t W\_[i,l]{}\^[-1]{}v\_[i,l]{}]{} x\_i.\[inverse-final\] We name the system $\cal W$ associated to the ODEs (\[inverse-driver\]), (\[inverse-final\]), an [*inverse quadratization*]{} of the original $\sigma\pi$-system.\

### Bernoulli differential equations.

In the second example of §\[sc:SCO\] we have met a scalar Driver ODE belonging to a known class of scalar differential equations, the [*Bernoulli differential equations*]{}, that are all solvable through a change of variables that turns them into linear equations. Recall that a Bernoulli ODE, in our notation, is a $\sigma\pi$-ODE with $n=1$ (it is scalar), $\nu=2$, $X_{11}=x$, $X_{12}=x^\alpha$ with $1\ne \alpha\in \re$: x=v\_1 x+v\_2 x\^.\[BEs\] Following the classic method, one divides (\[BEs\]) by $x^{\alpha}$ and obtains $
x^{-\alpha} \dot x= v_1 x^{1-\alpha} +v_2
$. Then, by defining $w=x^{1-\alpha}$, one has $\dot w=(1-\alpha) x^{-\alpha} \dot x$, and finally w=(1-)(v\_1w+v\_2).\[BEw\] Therefore, $w$ can be found as solution of the linear ODE (\[BEw\]) and $x=w^{1/(1-\alpha)}$.

We get the same result (even relaxing the hypothesis $\alpha\ne 1$) by applying to (\[BEs\]) the exact quadratization in [*inverse form*]{} given by the inverse Driver (\[inverse-driver\]). For, the Driver components associated to the $\sigma\pi$-ODE (\[BEs\]) are $Z_{11}=1$, and $Z_{12}=x^{\alpha-1}$, thus the inverse Driver components are $W_{11}=1$, $W_{12}=x^{1-\alpha}$, and we see that $W_{12}=w$. We get the ODE for $W_{12}$ directly from its monomial, $x^{1-\alpha}$, by applying the inverse-Driver formula (\[inverse-driver\]) : $\dot W_{12}= (1-\alpha)(v'Z)W_{12}=(1-\alpha)(v_1Z_{11}+v_2Z_{12})W_{12}$, then, since $Z_{11}=1$, and $Z_{12}W_{12}=1$ we have W\_[12]{}= (1-)(v\_1W\_[12]{}+v\_2),\[IDex\] which is the linear equation (\[BEw\]). For $\alpha=1$ the original ODE is linear, thus solvable with the known classic methods, and unlike Bernoulli method, that fails for $\alpha=1$, such solution is still a local sub-solution of the inverse quadratization.

We can say that a $\sigma\pi$-ODE is a generalized version, in multiple indeterminates, of the Bernoulli equation. Unfortunately though, besides the Bernoulli equation, there are very few particular cases, and not really meaningful, for which the inverse quadratization leads to the solution calculation via linearization.\

The Driver of a linear ODE
--------------------------

A linear ODE is a $\sigma\pi$-ODE where all $p_{i,j}^l$ in (\[pi\]) may be one or zero, and for any $i,l$ there is one and only one $j$ such that $p_{i,j}^l=1$. For a linear ODE we can always assume that, for any $i$, $\nu_i=n$, which amounts to saying that all the $n$ monomials: $X_{i,1}=x_1,\ldots, X_{i,n}=x_n$ of the $i$-th equation are actually written (otherwise, write it multiplied by a zero coefficient). This in turn entails that the monomials do not depend of the equation index $i$. In such a framework we have coefficients $v_{i,j}$, for $i,j=1,\ldots,n$ and $p^l_{i,j}=\delta_{j,l}$ for $i,j,l=1,\ldots,n$, (not depending of $i$). That said, consider a linear system: x\_i=\_[j=1]{}\^n v\_[i,j]{}x\_j, x\_i=\_[j=1]{}\^[n]{} v\_[i,j]{}X\_[i,j]{} X\_[i,j]{}=x\_j [)]{}.\[aff\] Since (\[aff\]) is a $\sigma\pi$-ODE, it undergoes a quadratization, and then it has an associated Driver, that we readily compute by the general Driver expression (\[dri\]). By (\[pi\]): $\pi^l_{i,j}=\delta_{j,l}-\delta_{i,j}$, which replaced in (\[dri\]) yields Z\_[i,l]{}=(v’\_lZ\_l-v\_i’Z\_i)Z\_[i,l]{} \[ldri\] The Driver lives on the domain ${\cal S}^*\subset \re^{n+d+1}$ given by (\[sbdom\]) with $\Phi$ given by Z\_[i,l]{}=\_[i,l]{}(x,t)=[x\_l x\_i]{}, l=1,…,n. \[Zaff\] As explained earlier, we get an inclusive quadratization by enlarging the Driver (\[ldri\]) with the further $n$ ODEs: Z\_[i,n+1]{}=(v\_i’Z\_i)Z\_[i,n+1]{},\[IDL\] with $Z_{i,n+1}=x_i$, which replaces the final equation.

Quadratic ODEs, Drivers and self-Drivers {#SD:sec}
----------------------------------------

We name [*Driver-type ODE*]{} (in $n$ indeterminates) an homogeneous quadratic ODE of the type x\_i=\_[j=1]{}\^[n]{} v\_[i,j]{}x\_ix\_[j]{}(=(v’\_i x) x\_i) \[SD\] where the $v_{i,j}$ are real analytic functions all defined on some open interval $\cal T\subset \re$. Therefore, a [*Driver-type system*]{} shall be the system $\re^n\times \cal T$ associated to the ODE (\[SD\]).

This name is motivated by the following theorem.\

\[th:Drty\] [*Any Driver-type ODE has always one sub-system that is the Driver of some $\sigma\pi$-system* ]{}

[*Proof.*]{} It is sufficient to prove that (\[SD\]) is the Driver ODE of some $\sigma\pi$-ODE. For, let us choose an integer $m\le n$, and $m$ integers $\nu_1,\ldots,\nu_m$ such that $\nu_1+\ldots,\nu_m=n$. Then we can rewrite (\[SD\]) as the following $m$-partite system with $k$th part $\nu_k$ dimensional, for $k=1,\ldots,m$: x\_[k,s]{}=\_[r,l=1]{}\^[m,\_r]{} v\_[k,s,r,l]{} x\_[k,s]{}x\_[r,l]{},\[mpart\] where in the summation $r=1,\ldots,m$ runs first, and $s=1,\ldots,\nu_r$, and we have simply replaced $i$ with $(k,s)$, and $j$ with $(r,l)$. Now, let us decompose $v_{k,s,r,l}$ as a product of real numbers $\pi^s_{k,r}$ and $v^*_{r,l}$: v\_[k,s,r,l]{}=\^s\_[k,r]{}v\^\*\_[r,l]{},\[mpartT1\] and define the real number $p^s_{k,r}$: p\^s\_[k,r]{}=\^s\_[k,r]{}+\_[i,j]{}. By using (\[dri\]), and (\[pi\]) one readily verify that (\[mpart\]) — with $x$ replaced by the symbol $Z$ — is the Driver ODE associated to the $\sigma\pi$-ODE with system function given by (\[SP-form1\]), where $m$ replaces $n$ and $v^*$ replaces $v$. $\bullet$\
Notice, by the proof of the theorem, that there are indeed many different $\sigma\pi$-systems generating the same Driver ODE (\[SD\]). Nonetheless, they generate each a different Driver as a different sub-system of the (same) Driver ODE, in other words, as a set of sub-solution (in the $Z$ indeterminates) lying in a different domain $\cal S^*$ — which is in fact realized by (\[sbdom\]) through a different map $\Phi^*$.

We name [*self-Driver*]{} (in $n$ indeterminates) any sub-system of the Driver-type ODE (\[SD\]). Thus, a sub-set $\cal F\subset \re^n\times \cal T$ invariant with respect to the ODE (\[SD\]), determines a self-Driver, and the ODE (\[SD\]) itself is a self-Driver (with domain $\re^n\times \cal T$). The set of all self-Drivers associated to the ODE (\[SD\]) that are Drivers (of some pre-Driver) is a sub-set (in general proper) of the set of all self-Drivers. Clearly, a self-Driver $\cal F\subset \re^n\times \cal T$ is a Driver if and only if there exists an integer $m\le n$, $m$ integers $\nu_1,\ldots,\nu_m$ such that $\nu_1+\ldots+\nu_m=n$, and $n$ monomial functions $X_{i,l}$, for $i=1,\ldots,m$ and $l=1,\ldots,\nu_i$, such that $\cal F=\Phi(\cal D_r\times \cal T)$ with $\Phi_{i,l}=x_i^{-1}X_{i,l}$.

Since (\[SD\]) is quadratic, it is a $\sigma\pi$-ODE as well, of order $n$ and with $n^2$ monomials. Thus, it undergoes a canonical quadratization and has in turn an associated Driver. The name ’self-Driver’ comes from the fact that any system associated with the ODE (\[SD\]) has a Driver whose solutions include, as sub-solutions, all the solutions of the system itself, and only these. Indeed, we can prove the following theorem.\

\[th:dristrut\] [*Let us consider a self-Driver $\cal F\subset \re^n\times \cal T$, associated to the ODE (\[SD\]). Then the Driver of $\cal F$ is defined by the following ODE $n$-partite in the $n^2$ indeterminates $Z_{i,j}$, $i,j=1,\ldots,n$: Z\_[i,j]{}=(v’\_jZ\_j)Z\_[i,j]{},\[SSD\] on the domain $\cal F^{\bf d}\subset \re^{n^2}\times \cal T$: F\^[**d**]{}=\_[(x,t)]{}{(Z,t)\^[n\^2]{}T/Z\_[i,j]{}=x\_i, (x,t)F}.\[DomSD\]* ]{}

[*Proof.*]{} Let us apply (\[dri\]) to (\[SD\]). We have $\pi^l_{i,j}=p^l_{i,j}-\delta_{i,j}$, where $p^l_{i,j}$ is the power of $x_j$ in the monomial $x_ix_l$, and hence is $1$ for $j=l$ or $j=i\ne l$. It follows that $\pi^l_{i,j}=1$ if and only if $j=l$. Therefore, the Driver ODE is (\[SSD\]). The Driver map $\Phi$, whose general formula is (\[Ph\]), applied to the $\sigma\pi$-ODE (\[SD\]) gives: $Z_{i,l}=\Phi(x)=x_{l}$, and hence (\[DomSD\]) is an invariant set for the ODE (\[SSD\]), thus the Driver domain. $\bullet$\
Although the Driver of the self-Driver $\cal F$ has $n^2$ indeterminates, any solution has many equal components, so that just suitably picking up $n$ among the $n^2$ solution components, one get the original self-Driver solution $(x,t)\in \cal F$. This motivates the following identification, stated by the definition below, of a self-Driver with its Driver, which amounts to consider the latter just a different (augmented) [*representation*]{} of the former.\

\[df:augR\] [*Given a self-Driver $\cal F\subset \re^n\times \cal T$ described by the ODE (\[SD\]), we call the system (in $\re^{n^2}\times \cal T$) defined by the ODE (\[SSD\]) on the domain $\cal F^{\bf d}$ given by (\[DomSD\]) the ’augmented representation’ of the self-driver $\cal F$.*]{}

0.8cm

We finally give the following definition, acconding to which, s Driver-type ODE is characterized by its own matrix of coefficients $v_{i,j}$-\

\[df:frame\] [*We name ’the frame’ of a Driver-type ODE in $n$ indeterminates, the matrix $V\in\re^{n\times n}$, $V=(v_{i,j})$, collecting all of its own coefficients.* ]{}

0.8cm

The solution formula for Driver-type ODEs {#sc:SDTheo}
=========================================

Statement of the main theorem {#ssc:mai}
-----------------------------

The importance of Driver type quadratic ODEs, and, hence, of the exact quadratization, relies on the result stated in the next theorem, showing that all solutions of a Driver type ODE can be calculated through power series.

Let us consider an ODE of the type (\[SD\]) with analytic coefficients defined on the interval $\cal T\subset \re$. The associated [*full*]{} system is then $(\re^n\times \cal T, \{1,\ldots,n\})$, but on account of the properties of analytic systems, described in §\[ssc:sinreg\], without loss of generality we can consider the sub-system $(\cal D\times \cal T, \{1,\ldots,n\})$, where $\cal D=\re^n\setminus S^0$, and $S^0$ is the set of all coordinate hyperplanes in $\re^n$. Indeed, any solution $(x,t)\in \re^n\times \cal T$ with $x_i=0$ for some $i$, is either a regular solution, and, hence, it can be rewritten as $(x^*,t^*)\in \cal D\times \cal T$, or is singular, and then it can be rewritten as $(x^*,t)\in \cal D^*\times \cal T$, where $\cal D^*\times \cal T$ is the domain of system projected on $\{1,\ldots, n\}\setminus \{i\}$. Besides, if the system is a self-Driver, its domain, say $\cal D\times \cal T\subset \re^n\times \cal T$, is the image: $\cal D=\Phi(\cal D^*)$ where $\cal D^*\in \re^m$, for some $m\le n$, $\cal D^*\cap D^0=\emptyset$ ($D^0$ denoting now the set of points of all hypherplanes in $\re^m$), and $\Phi$ is an $\re^n$-valued map whose components are all [*monomials*]{} of the type $z^{p_1}_1\dots z^{p_m}_m$ with $z_1,\ldots, z_m$ real indeterminates, and $p_1,\ldots p_m$ real exponents. Thus, any component of $\Phi$ can be zero only on $D^0$, which implies that $\cal D\cap S^0=\emptyset$. The hypothesis of a domain not including any point of the coordinate hypherplanes, is then directly verified for the more expected situation of our study: i.e. when the ODE is a quadratization of some original $\sigma\pi$-reducible system.\
With that being said, for the solution $(x,t_0)$ of the ODE (\[SD\]) with coefficients $v_{i,j}$, where $(x,t_0)\in \cal D\times \cal T$, we can consider the Taylor expansion of the corresponding curve-solution, $x(t)$, with $t\in {\bf I}_{x,t_0}\subset \cal T$: x\_i(t)=\_[k=0]{}\^c\_k(i) [(t-t\_0)\^kk!]{},c\_k(i)=x\^[(k)]{}(t\_0), \[SPSD\] where $x^{(k)}(\cdot)$, defined for all $t\in {\bf I}_{x,t_0}$, denotes the curve of the $k$th time-derivative of $x(\cdot)$. Identity (\[SPSD\]) holds for any $t\in {\bf I}_{x,t_0}$ such that the series in the right hand side converges, which happens in general on a sub-interval: (t\_0,r)\_[x,t\_0]{},\[intconv\] where ${\bf I}(t_0,r)$ is the open interval centered in $t_0$ with ray $r>0$.

The following theorem gives the general formula for the coefficients $c_k(i)$. In the following the partial time-derivative operator $\partial_t$ will be used in place of the ’dot’ notation for symbols including superscripts, thus $\partial_t v^k$ in place of $\dot {\bf v}$ for ${\bf v}=v^k$.\

\[th:serie\] [*The coefficients $c_k(i)$ in the Taylor expansion (\[SPSD\]) are given by &&c\_0(i)=x\_i,\[GE0\]\
&&c\_k(i)=\_[s=2]{}\^[k+1]{}\_[i\_1,…,i\_[s-1]{}I]{} v\^[k+1,s]{}\_[i,i\_1,…,i\_[s-1]{}]{}(t\_0) x\_ix\_[i\_1]{}x\_[i\_[s-1]{}]{}, I={1,…,n},k1,\[GEs\] where $i_l\in \cal I$, for $l=0,1,\ldots, k$ (we set $i_0=i$), and the coefficients $v^{k+1,s}_{i,i_1,\ldots,i_{s-1}}$, ($=v^{k+1,s}_{i,i_1,\ldots,i_{s-1}}(t)$, $t\in {\bf I}_{x,t_0}$) are given by the following recursive equations ($k\ge 2$): &&v\^[k+1,k+1]{}\_[i,i\_1,…,i\_k]{}=v\^[k,k]{}\_[i,i\_1,…,i\_[k-1]{}]{}(\_[j=0]{}\^[k-1]{} v\_[i\_j, i\_k]{}), \[Vkk\]\
&&v\^[k+1,s]{}\_[i,i\_1,…,i\_[s-1]{}]{}=v\^[k,s-1]{}\_[i,i\_1,…,i\_[s-2]{}]{}(\_[j=0]{}\^[s-2]{} v\_[i\_j, i\_[s-1]{}]{})+\_t v\^[k,s]{}\_[i,i\_1,…,i\_[s-1]{}]{},s=3,…,k,\[Vks\]\
&&v\^[k+1,2]{}\_[i,i\_1]{}=\_t v\^[k,2]{}\_[i,i\_1]{}, ( = \_t\^[k+1]{} v\_[i,i\_1]{}),\[Vk2\] which are initialized by v\^[2,2]{}\_[i,i\_1]{}=v\_[i,i\_1]{}.\[initM\]* ]{}

0.8cm

### Comments to the statement of Theorem \[th:serie\] {#ssc:com}

for any $k,s$, the number of coefficients $v^{k,s}_{i,i_1,\ldots,i_{s-1}}$ is $n^s$. in order to calculate $c_k(i)$, $i$ is fixed, and one needs $n+\ldots +n^k$ coefficients, one for any fixed sequence $i_1,\ldots,i_{s-1}\in \cal I$, $s=2,\ldots, k+1$. The coefficients have to be multiplied each for the corresponding monomial $x_ix_{i_1}\cdots x_{i_{s-1}}$ obtained by the components of the initial value $x=x(t_0)$. By using the index convention we have $v^{k,s}_i(t)\in \re^{n^s-1}$, and we can write (\[GEs\]) in vector notation as well: $$c_k(i)=\sum_{s=2}^{k+1}\left(v_i^{k,s} (t_0)\right)' x^{\otimes s-1} x_i$$ where $x^{\otimes l}$ denotes the $l$-times tensor product $x\otimes \ldots\otimes x$. Note that the scheme (\[Vkk\])–(\[initM\]) is recursive in $k$: from the set of all coefficients involved in $c_{k-1}(i)$, namely: $$\cal V_i(k-1)=\{ v^{r,s}_i\quad |\quad r=2,\ldots,k,\quad s=2,\ldots, r\}$$ one directly calculates from (\[Vkk\])–(\[Vk2\]) the set $\{v_i^{k+1,s}, \quad | s=2,\ldots, k+1\}$, which, joined to $\cal V_i(k-1)$, yields $\cal V_i(k)$. Note that the coefficients $v^{k,s}_i$ can be calculated independently of the initial point $x$. If formula (\[SPSD\]) is used numerically, then for any fixed approximation degree, and initial point $x\in \cal D$, it allows to calculate separately each component $x_i(t)$ of the solution, within the given approximation, and within the convergence interval ${\bf I}(t_0,r)$. One is able as well to calculate a possible extension of the solution outside ${\bf I}(t_0,r)$, by moving the initial time to another $t_0'\in {\bf I}(t_0,r)$ – provided that $x(t'_0)\in \cal D$, and ${\bf I}(t_0',r')\cup {\bf I}(t_0,r)\setminus {\bf I}(t_0,r)\ne \emptyset$ – and calculating the value at $t_0'$ of always the same unique pre-calculated set of coefficients $v_i^{k,s}$.

A preliminary result {#ssc:prel}
--------------------

In order to prove Theorem \[th:serie\] we need some additional definitions. Let us fix $i\in \cal I$, where $\cal I$ is the set of indices in (\[GEs\]). For any integer $k\ge 2$, we define $\cal P_i^k$ to be the set of all homogeneous polynomials of degree $k$ of the $n$ indeterminates $x_1,\ldots,x_n$, with $x_i$ always in, and where the coefficients are time-functions, defined on some open real interval, and having time-derivative of any order. Therefore, $p_i^k\in \cal P_i^k$ has the form: p\_i\^k(x)=\_[i\_1,…,i\_[k-1]{}I]{} v\^k\_[i,i\_1,…,i\_[k-1]{}]{} x\_ix\_[i\_1]{}…, x\_[i\_[k-1]{}]{}.\[pk\] Now, let us consider a Driver-type ODE with frame $V$, and, related to $V$, define a map $A_i:\cal P_i^k\to \cal P_i^{k+1}$, for any $k\ge 2$ as follows: A\_ip\_i\^k=p\_i\^[k+1]{}=\_[i\_1,…,i\_kI]{} v\^[k+1]{}\_[i,i\_1,…,i\_[k]{}]{} x\_ix\_[i\_1]{}…, x\_[i\_[k]{}]{}, v\^[k+1]{}\_[i,i\_1,…,i\_[k]{}]{}=v\^[k]{}\_[i,i\_1,…,i\_[k-1]{}]{}(\_[j=0]{}\^[k-1]{}v\_[i\_j,i\_k]{}),\[dA\] where $v_{i_j,k}$ are coefficients of the frame $V$, and conventionally we set $i_0=i$. Note that any $\cal P_i^k$ is a vector space, with naturally defined operations and null vector, thus we can define $\cal P_i$ as the (infinite) direct sum of all $\cal P_i^k$’s: $
\cal P_i=\oplus_{s=0}^{+\infty} \cal P_i^s$. Therefore, for any $q\in\cal P$ there is a $k\ge 2$ and $p_i^s\in \cal P_i^s$ for $s=2,\ldots,k$ such that $q_i=\sum_{s=2}^kp^s$, and we say that $k$ is the degree of $q_i$. We can naturally extend as follows the definition of $A_i$ given earlier, so that $A_i:\cal P\to \cal P$: let $q_i^k\in \cal P_i$ a $k$-degree polynomial, then $q_i^k=p_i^2+\ldots+p_i^k$, for some $p_i^s\in \cal P_i^s$, then we define $
A_iq_i^k=\sum_{s=2}^k A_ip_i^s
$. The exponent $A_i^s$ is defined at once as a composition of maps, and we have $A_i^s: \cal P_i^k\to \cal P_i^{k+s}$ for any integer $s\ge 0$, provided we set $A_i^0=I$ (the identity map). Besides $A_i$ we make use of the partial time-derivative operator $\partial_t$, for which a similar convention is applied: $\partial_t^s$, for $s\ge 0$ shall denote composition, and $\partial^0_t$ is assumed to be the identity. Note that $\partial^s_t: \cal P_i^k\to \cal P_i^k$ for any $k\ge 2$, $i\in \cal I$, and we have: \^s\_tp\_i\^k=\_[i\_1,…,i\_[k-1]{}]{} \_t v\^[k]{}\_[i,i\_1,…,i\_[k-1]{}]{} x\_ix\_[i\_1]{}…, x\_[i\_[k-1]{}]{}.\[papk\] In the same way all compositions of the maps $A$ and $\partial_t$ are well defined, as well as the map $(A+\partial_t)^s$ for any $s\ge 0$, and we have $(A_i+\partial_t)^s \cal P_i\subset \cal P_i$, and $(A_i+\partial_t)^s \cal P_i^k\subset \cal P_i^{k+s}$.

Now, hereinafter we fix the symbol $p^2_i\in \cal P^2$ to be the polynomial whose coefficients are just the components of the frame $V$: p\^2\_i(x)=\_[i\_1]{} v\_[i,i\_1]{} x\_i(t)x\_[i\_1]{}(t), \[p2\] which, in the formalism defined above means that we are taking $v^{2}=v$. It follows at once that we can rewrite a Driver-type ODE with frame $V$ as: x\_i\^[(1)]{}=p\^2\_i(x). \[DTOrwt\] The proof of Theorem \[th:serie\] requires the following result, that we single out in a lemma because it has a significance in its own right.\

\[lm:serie\] [*There exists an open interval \^\*\_[x,t\_0]{}\_[x,t\_0]{},\[opin\] with $t_0\in {\bf I}^*_{x,t_0}$, such that the $k$th time-derivative of the curve $x(\cdot)$ of the solution $(x,t_0)\in \cal D\times \cal T$ of (\[SD\]), is given component-wise by x\_i\^[(k)]{}(t)=q\_i\^[k+1]{}(x(t)) k1,t\^\*\_[x,t\_0]{},\[kder\] where $q_i^{k}\in \cal P_i$, is a $k$th degree polynomial for any $k\ge 2$ satisfying the recursive equation: q\_i\^[k+1]{}=(A\_i+\_t) q\_i\^k, q\_i\^2=p\_i\^2.\[qkrec\]* ]{}

0.8cm

[*Proof.*]{} Since $t\mapsto x(t)$ is continuous in ${\bf I}_{x,t_0}$, and $x_i\ne 0$ $\forall i$, there exists an interval ${\bf I}^*_{x,t_0}\subset {\bf I}_{x,t_0}$, with $t_0\in {\bf I}^*_{x,t_0}$, such that x\_i(t)0, t\^\*\_[x,t\_0]{}.\[xine0\] Note that the recursive equation (\[qkrec\]) is equivalent to q\_i\^[k+1]{}=(A\_i+\_t)\^[k-1]{}p\_i\^2,\[qnonrec\] therefore, (\[kder\]) is equivalent to x\_i\^[(k)]{}(t)=(A\_i+\_t)\^[k-1]{}p\_i\^2(x(t)), k1,t\^\*\_[x,t\_0]{}.\[kderts\] We prove the lemma by induction. For $k=1$ the lemma holds, since (\[kderts\]) turns into (\[DTOrwt\]). Suppose (\[kderts\]) is true, and let us prove (\[kderts\]) for $k+1$ replacing $k$. We have ($D_t$ denoting the total time-derivative, replacing the [*dot*]{} notation upon symbols with other superscripts): x\_i\^[(k+1)]{}(t)=D\_[t]{} x\_i\^[(k)]{} (t)=\_[l=1]{}\^n (\_[x\_l]{} x\_i\^[(k)]{})|\_[x=x(t)]{} x\_l(t)+ \_t x\^[(k)]{}\_i |\_[x=x(t)]{}\[lmM1\] As for the last term, by using the inductive hypothesis, it results in \_t x\^[(k)]{}|\_[x=x(t)]{}=\_t(A\_i+\_t)\^[k-1]{} p\^2\_i(x(t)).\[lmM2\] By (\[qnonrec\]), $q_i^{k+1}$ is a polynomial of degree $k+1$ in $\cal P_i$, therefore it has the form: q\_i\^[k+1]{}=\_[s=2]{}\^[k+1]{}p\^s\_i, p\^s\_iP\^s,\[dq\] where the $p_i^s$’s are homogeneous polynomials of the type (\[pk\]). Thus, the first term in the left hand side of (\[lmM1\]) becomes, by replacing the polynomial expression (\[dq\]) for $q^{k+1}_i$: \_[l=1]{}\^n (\_[x\_l]{} x\_i\^[(k)]{})|\_[x=x(t)]{} x\_l(t) = \_[s=2]{}\^[k+1]{}\_[i\_1,…,i\_[s-1]{}]{} (v\^[s]{}\_[i,i\_1,…,i\_[s-1]{}]{} \_[l=1]{}\^n \_[x\_[l]{}]{} (x\_ix\_[i\_1]{}x\_[i\_[s-1]{}]{})|\_[x=x(t)]{}x\_i(t)) \[dxtmp\] Now, on account of (\[xine0\]), we have \_[l=1]{}\^n \_[x\_[l]{}]{} (x\_ix\_[i\_1]{}x\_[i\_[s-1]{}]{})|\_[x=x(t)]{} x\_l(t) &&= \_[j=0]{}\^[s-1]{} \_[x\_[i\_j]{}]{} (x\_ix\_[i\_1]{}x\_[i\_[s-1]{}]{}) |\_[x=x(t)]{}x\_[i\_j]{}(t),\
&&= \_[j=0]{}\^[s-1]{} x\_ix\_[i\_1]{}x\_[i\_k]{} x\^[-1]{}\_[i\_j]{} |\_[x=x(t)]{} x\_[i\_j]{}(t),thus, by replacing $\dot x_{i_j}$ with the right hand side of (\[p2\]) where we also rename $i_1$ after the new index $i_{s}$, we can write &&-1.5cm \_[i\_1,…,i\_[s-1]{}]{} (v\^[s]{}\_[i,i\_1,…,i\_[s-1]{}]{} \_[j=0]{}\^[s-1]{} x\_ix\_[i\_1]{}x\_[i\_[s-1]{}]{} x\^[-1]{}\_[i\_j]{} \_[i\_[s]{}]{} v\_[i\_j,i\_[s]{}]{} x\_[i\_j]{}x\_[i\_[s]{}]{},)|\_[x=x(t)]{}\
&&=\_[i\_1,…,i\_[s]{}]{} (v\^[s]{}\_[i,i\_1,…,i\_[s-1]{}]{} x\_ix\_[i\_1]{}x\_[i\_[s-1]{}]{} \_[j=0]{}\^[s-1]{} x\^[-1]{}\_[i\_j]{} v\_[i\_j,i\_[s]{}]{} x\_[i\_j]{}x\_[i\_[s]{}]{},)|\_[x=x(t)]{}\
&&=\_[i\_1,…,i\_[s]{}]{} (v\^[s]{}\_[i,i\_1,…,i\_[s-1]{}]{} x\_ix\_[i\_1]{}x\_[i\_[s]{}]{}) |\_[x=x(t)]{} (\_[j=0]{}\^[s-1]{} v\_[i\_j,i\_[s]{}]{}).\
&&=A\_ip\^[s]{}(x(t)), where, in the last step, the definition of the map $A_i$, given in (\[dA\]) has been used. By replacing the above identity in (\[dxtmp\]), and using again the inductive hypothesis – in the form (\[kder\]) – we obtain $$\sum_{l=1}^n \left(\partial_{x_l} x_i^{(k)}\right)\Big|_{x=x(t)} \dot x_l(t)= \sum_{s=2}^{k+1}A_ip^s_i=A_iq^{k+1}_i=
A_i(A_i+\partial_t)^{k-1}p^2_i(x(t))$$ which used in (\[lmM1\]) at once with (\[lmM2\]) gives $$x_i^{(k+1)}(t)=(A_i+\partial_t)^k p^2_i(x(t)),$$ Q.E.D.$\bullet$\

Proof of Theorem \[th:serie\]. {#scc:pro}
------------------------------

By Lemma \[lm:serie\] the coefficients $c_k(i)$ are given by c\_k(i)=x\^[(k)]{}(t\_0)=q\^[k+1]{}\_i(x), k0\[ck\] and since $q^{k+1}_i$ is a $k+1$ degree polynomial in $\cal P_i$, there exist $k-1$ polynomials $p^{k+1,s}_i\in \cal P_i^{s}$, $s=2,\ldots,k+1$ such that q\^[k+1]{}\_i=\_[s=2]{}\^[k+1]{}p\^[k+1,s]{}\_i=\_[s=2]{}\^[k+1]{} \_[i\_1,…,i\_[s-1]{}]{} v\^[k+1,s]{}\_[i,i\_1,…,i\_[s-1]{}]{} x\_ix\_[i\_1]{}…x\_[i\_[s-1]{}]{}\[genfpk\] which gives the general expression (\[GEs\]). In order to proof formulas (\[Vkk\])–(\[Vk2\]), consider the recursive formula (\[qkrec\]). where $q^k_i$ has a similar expression like (\[genfpk\]) but the degree is $k$: q\^[k]{}\_i=\_[s=2]{}\^[k]{}p\^[k,s]{}\_i=\_[s=2]{}\^[k]{} \_[i\_1,…,i\_[s-1]{}]{} v\^[k,s]{}\_[i,i\_1,…,i\_[s-1]{}]{} x\_ix\_[i\_1]{}…x\_[i\_[s-1]{}]{}\[genfpk2\] thus we have (omit all subscripts)

[rl]{} q\^[k+1]{}\_i=p\^[k+1,k+1]{}+p\^[k+1,k]{}+…+p\^[k+1,2]{}&=Ap\^[k,k]{}+Ap\^[k,k-1]{}+…,Ap\^[k,2]{}\
&+p\^[k,k]{}+p\^[k,k-1]{}+…,p\^[k,2]{}

\[q=p\] where in the right hand side, by definition given in (\[dA\]), we have that $Ap^{k,s}\in \cal P^{s+1}$, and $\partial p^{k,s}\in \cal P^s$, for $s=2,\ldots,k$. By equating the homogeneous polynomials with the same degree, (\[q=p\]) is equivalent to the following set of equalities: p\^[k+1,k+1]{}&=&Ap\^[k,k]{},\
p\^[k+1,s]{}&=&Ap\^[k,s-1]{}+p\^[k,s]{},s=3,…, k\
p\^[k+1,2]{}&=&p\^[k,2]{}which directly gives equalities (\[Vkk\])–(\[Vk2\]), as soon as the definition of A, given in (\[dA\]) is applied to the polynomial expression of $p^{k,s}$ and $p^{k+1,s}$, for $s=2,\ldots k+1$, which can be read out from (\[genfpk2\]) and (\[genfpk\]). $\bullet$\
Let $\cal S\subset \{1,\ldots, n\}$ be the set of indices: S ={ j{1,…,n}| i, v\_[i,j]{}0},\[defSc\] In other words, the set $\cal S$ is the set of indices addressing the non-zero columns of the frame $V$. The following theorem can ease the calculation of formula (\[GEs\]).\

\[th:S\] [*In formula (\[GEs\]) the indices $i,i_1,\ldots,i_{s-1}$, don’t have to run across all values $\{1,\ldots, n\}$: it’s enough they span just the subset $\cal S$ defined in (\[defSc\]).* ]{}

0.8cm

[*Proof.*]{} The thesis is equivalent to say that for any $k\ge 2$ and $2\le s\le k$, the following condition is satisfied:\
[*Condition (C)*]{}: $v^{k,s}_{i,i_1,\ldots,i_{s-1}}=0$ if there is an $m$, $1\le m\le s-1$ such that $i_{m}\not \in \cal S$.\
Let is denote by ${\bf v}(l)$, with $l\in \na$, the set of all pair $(k,s)$ with $k\ge 2$, $2\le s\le k$, and $k-s=l$, then the theorem is proven as soon as we prove, by induction over $l\in \na$, that condition (C) is verified by all pairs in $\cup_{l\in \na}{\bf v}(l) $. For, let us consider first the pairs of ${\bf v(0)}$, which address the coefficients $v^{k,k}$ generated by (\[Vkk\]). Suppose that there is an index, say $i_m$, $1\le m\le s-1< k$, such that $i_m\not\in \cal S$. Then $v_{j,i_m}=0$ (for any $j$) and by formula (\[Vkk\]) for $k=m$, we have $$v^{m+1,m+1}_{i,i_1,\ldots,i_m}=v^{m,m}_{i,i_1,\ldots,i_{m-1}} \left(\sum_{j=0}^{m-1} v_{i_j, i_m}\right)=0,$$ and, hence, again by the recursive formula (\[Vkk\]): $$v^{k+1,k+1}_{i,i_1,\ldots,i_k}=v^{m+1,m+1}_{i,i_1,\ldots,i_{m}} \left(\sum_{j=0}^{m} v_{i_j, i_{m+1}}\right)\cdots 
 \left(\sum_{j=0}^{k} v_{i_j, i_{k+1}}\right)=0.$$ Now, suppose that condition (C) is verified by all pairs in ${\bf v}(l)$, and let us prove that it is verified by all pairs of ${\bf v}(l+1)$ as well. For, suppose that $k-s=l$ and there is an $i_m$, $1\le m\le s-1< k$, such that $i_m\not\in \cal S$, then, by the inductive hypothesis, we have $v^{k,s}_{i,i_1,\ldots,i_{s-1}}=0$, which simplifies the recursive formula (\[Vks\]) as follows: $$v^{k+1,s}_{i,i_1,\ldots,i_{s-1}}=v^{k,s-1}_{i,i_1,\ldots,i_{s-2}}\left(\sum_{j=0}^{s-2} v_{i_j, i_{s-1}}\right),$$ from which, by successive substitutions up to $m+1$ we obtain $$v^{k+1,s}_{i,i_1,\ldots,i_{s-1}}=v^{m+l+1,m}_{i,i_1,\ldots,i_{m-1}} \left(\sum_{j=0}^{s-2} v_{i_j, i_{s-1}}\right)\cdots 
 \left(\sum_{j=0}^{m-1} v_{i_j, i_{m}}\right)=0.$$ since, again by the inductive hypothesis, $v_{j,i_m}=0$ for all $j$. We have proven that $v^{k+1,s}=0$ for $k-s=l$, which means $v^{k,s}=0$ for $k-s=l+1$, i.e. for all pairs in ${\bf v}(l+1)$, Q.E.D. $\bullet$\

The stationary case {#ssc:statcase}
-------------------

Worth of mention is the [*stationary*]{} case, i.e. when the coefficients $v_{i,l}$ are all constants, because in this case the calculation of the coefficients $c_k(i)$ greatly simplifies, as shown in the following theorem\

\[th:seriest\] [*In the stationary case, the coefficients $c_k(i)$ in the Taylor expansion (\[SPSD\]) are constant and are given by c\_k(i)=\_[i\_1,…,i\_kS]{} v\^[k+1]{}\_[i,i\_1,…,i\_k]{}x\_ix\_[i\_1]{}x\_[i\_k]{},\[cstat\] where $i_s=1,\ldots, n$, for $s=0,1,\ldots, k$ (we set $i_0=i$), and the (constant) coefficients $v^{k+1}_{i,i_1,\ldots,i_k}$, are given by the following recursive equation v\^[k+1]{}\_[i,i\_1,…,i\_k]{}=v\^[k]{}\_[i,i\_1,…,i\_[k-1]{}]{}(\_[j=0]{}\^[k-1]{} v\_[i\_j, i\_k]{}), v\^1\_i=1.\[Vkstat\]* ]{}

0.8cm

[*Proof.*]{} Let us prove that $v^{k,s}=0$, for $k-s>0$ (note that $k-s<0$ doesn’t apply). The proof is by induction over $k-s$, thus consider first the case $k-s=1$, which amounts to prove that $v^{k+1,k}=0$, $\forall k>2$. For, let us proceed by induction over $k$: for $k=2$ by (\[Vk2\]) we have $v^{3,2}=\partial_t v^{2,2}=0$ by the stationarity hypothesis. Suppose $v^{k,k-1}=0$, then (\[Vks\]) for $s=k$ yields $v^{k+1,k}=\partial _t v^{k,k}=0$, since by (\[Vkk\]) initialized by (\[Vk2\]), and by stationarity, we have $\partial_t v^{k,k}=0$. Coming back to the main induction, let us suppose that for $k-s=l$, $v^{k,s}=0$. Then by (\[Vks\]): v\^[k+1,s]{}\_[i,i\_1,…,i\_[s-1]{}]{}=v\^[k,s-1]{}\_[i,i\_1,…,i\_[s-1]{}]{}(\_[j=0]{}\^[s-2]{}v\_[i\_j,i\_[s-1]{}]{}),\[k-s=l\] using which, we can prove that $v^{k,s}=0$ for $k-s=l+1$, by induction over $s$: for $s=2$ we have $k=l+3$, and $v^{l+3,2}=0$ by (\[Vk2\]) and by stationarity. Then, suppose $v^{k,s-1}=0$, and $k-s+1=l+1$, then $k-s=l$ and by (\[k-s=l\]) we have $v^{k+1,s}=0$ and $k+1-s=l+1$, which completes the induction for $k-s>0$. In conclusion, there remains only the coefficients $v^{k,k}$ given by (\[Vkk\]), initialized by (\[initM\]): the thesis follows, by renaming $v^{k,k}\to v^k$, and extending $v$ with $v^1_i=1$ – which makes (\[Vkstat\]) returning (\[initM\]) for $k=1$. $\bullet$\
Let $\alpha^l_{i,i_1,\ldots,i_{m}}$ for $m\in \na$ be the number of times that $l=i_j$ for $j=0,1,\ldots, m$. Thus, given the sequence of indices $i,i_1,\ldots,i_{m}$, we have recursively: \^l\_[i,i\_1,…,i\_m]{} = \^l\_[i,i\_1,…,i\_[m-1]{}]{}+\_[l,i\_m]{}\^l\_i=\_[l,i]{},\[alpharec\] and \^l\_[i,i\_1,…,i\_m]{}=\_[j=0]{}\^m \_[l,i\_j]{}\[alpha\] It is possible to reduce the computational burden of formulas (\[Vkstat\]) and (\[cstat\]) by avoiding the calculation of many zero terms, as shown in the following corollary.\

\[co:seriest\] [*Formula (\[Vkstat\]) can be rewritten as follows: &&v\^[k+1]{}\_[i,i\_1,…,i\_k]{}=v\^[k]{}\_[i,i\_1,…,i\_[k-1]{}]{} \_[i,i\_1,…,i\_k]{} v\^1\_i=1,\[Vkd\]\
&&\_[i,i\_1,…,i\_k]{}=\_[l(i\_k){i}]{} \^[l]{}\_[i,i\_1,…,i\_[k-1]{}]{} v\_[l, i\_k]{}.\[Vkb\] where, for any $j\in \cal S$: (j) ={ lS | v\_[l,j]{}0}\[rho\]* ]{}

0.8cm

[*Proof.*]{} Define (j) ={l{1,…,n}| v\_[l,j]{}0}.We see in (\[Vkstat\]) that the index $i_j$ for $j\ne 0$ (thus, with the exception of $i$) doesn’t have to take values that are not in $\tilde \rho(i_k)$. Since $i_j$, for $j\ne 0$, doesn’t have as well to take values that are not in $\cal S$, the thesis follows. $\bullet$\
Denote by $\sigma$ the cardinality of $\cal S$, $v_M=\max_{i,j\in {\cal I}} |v_{i,j}|$, and for $x\in \cal D$, (thus, $\forall i$, $x_i\ne 0$) $x_M=\max_i |x_i|$. In the following Theorem, we find a lower bound, $\bar r$, for $r$, the convergence ray of the series for all the functions $t\mapsto |x_i(t)|$, components modules of the solution $(x, t_0)$, as well as an upper bound for all of them, in the interval $(t_0-\bar r, t_0+\bar r)\subset {\bf I}(t_0,r)$.\

\[th:seriest\] [ *For all $t\in (t_0-\bar r, t_0+\bar r)\subset {\bf I}(t_0,r)$, with |r=[1v\_Mx\_M]{},\[LB\] and for any $i=1,\ldots,n$ we have: |x(t)|, \[UB\]* ]{}

0.8cm

[*Proof.*]{} By (\[Vkstat\]) we have $|v^{k+1}|\le |v^k| k v_{M}$ which, iterated, yields (on account that $v^1=1$): $$|x(t)| \le k! v_M^k.$$ Thus, by (\[cstat\]), we have $$|c_k(i)| \le \sigma^k k! v_M^kx_M^{k+1}$$ which gives a geometric series whose generator is $\sigma v_Mx_M t$, and we get the upper-bound (\[UB\]). The lower-bound $\bar r$ in (\[LB\]), is readily derived by an application of the ratio criterion. $\bullet$\

Examples {#ssc:exfin}
--------

### The basic linear ODE {#sss:blo}

$\dot x=ax$, with $a\in \re$. The solution $(x,0)\in \re^2$ is $x(t)=x e^{at}$. Then, let us verify that the Taylor expansion of the solution, calculated by using Theorem \[th:seriest\] is: x(t)=x\_[k=0]{}\^.\[ex-eat\] For, we have first to quadratize the ODE, which is readily done in this simple case by defining the two indeterminates $x_1, x_2$ and setting

[l]{} x\_1=ax\_1x\_2,\
x\_2=0,

\[ex-newe\] from which we can recover the original solution as the sub-solution (in $x_1$) of the solution $( (x_1,1), 0)$ of (\[ex-newe\]). Now, the frame of \[ex-newe\] is V=(

[cc]{} 0& a\
0&0

)\[Vex\] thus $\cal S=\{2\}$, and by setting $i=1$, the indices $i,i_1,\ldots,i_k$ in (\[cstat\]) can take just the single value $1,2_{k}$, where we denote $2_k=2,\ldots,2$ ($k$-times). Therefore, formula (\[Vkstat\]) yields: $$v^{k+1}_{1,2_k}=v^k_{1,2_{k-1}} (v_{1,2}+kv_{2,2})=v^k_{1,2_{k-1}}a=a^k$$ and by (\[cstat\]), recalling that $x_2=1$: $$x_1(t)=\sum_{k=0}^\infty v^{k+1}_{1,2_k} x_1 x_2\cdots x_2 {t^k\over k!}=x_1\sum_{k=0}^\infty {a^kt^k\over k!}.$$ which is (\[ex-eat\]).\

### A simple linear non-stationary ODE. {#sss:blons}

Consider the ODE: $\dot x=2t x$, for which the solution $(x,0)$ can be easily calculated as x(t)=xe\^[\_0\^t 2d]{}=xe\^[t\^2]{},\[ex2-eat2\] whose Taylor expansion centered at $0$ is: x(t)=x\_[k=0]{}\^=x(1+ t\^2+[t\^42]{} + [t\^63!]{}+…),\[TE\] To illustrate the point, we now apply Theorem \[th:serie\] in order to calculate the first four coefficients of the power series in (\[TE\]). Note that we get the same quadratization (\[ex-newe\]), and frame (\[Vex\]), as in the previous example, provided $2t$ replaces $a$, therefore we have to calculate the sub–solution $x_1$ of solution $((x_1,1),0)$ of the quadratization. We have again ${\cal S}=\{2\}$, and $i=1$, thus, on account of Theorem \[th:S\], in (\[GEs\]) the summation in $1,i_1,\ldots, i_{s-1}$ has only one element, addressed by the string $1,2_{s-1}$, but the sum in $s$ has $k-1$ elements: c\_k(1)=(v\_[1,2]{}\^[k+1,2]{}+v\_[1,2,2]{}\^[k+1,3]{}+…+v\_[1,2\_[k]{}]{}\^[k+1,k+1]{}) |\_[t=0]{}x\_1.\[ex2-gf\] For $k=1$ we have $c_1(1)=x_1$ directly from (\[GE0\]). For any $k\ge 2$ by (\[Vkk\]) we readily have v\^[k,k]{}\_[1,2\_[k-1]{}]{}=(2t)\^[k-1]{}.\[vkk12\] By (\[ex2-gf\]), (\[vkk12\]), and (\[Vk2\]), (\[initM\]) we have: $$c_2(1)=(v_{1,2}^{3,2}+v_{1,2}^{2,2})\big|_{t=0}x_1=(\partial_t (2t) +2t)\big |_{t=0}x_1=(2+2t)\big|_{t=0}x_1=2x_1$$ thus the coefficient of $t^2$ is $2x_1/2=x_1$. By (\[Vks\]): $$v^{4,3}_{1,2,2}=v^{3,2}_{1,2}(2t)+\partial_t v^{3,3}_{1,2,2}=\partial_t(2t) (2t)+\partial_t (2t)^2=2(2t)+4(2t) =12t$$ therefore $$c_3(1)=(v_{1,2}^{4,2}+v_{1,2,2}^{4,3}+v_{1,2,2,2}^{4,4})\big|_{t=0}x_1=(\partial^2_t (2t) +12t+ (2t)^3)\big |_{t=0}x_1=0.$$ which is the coefficient of $t^3$ in (\[TE\]). Let us calculate just another coefficient. $$c_4(1)=(v_{1,2}^{5,2}+v_{1,2,2}^{5,3}+v_{1,2,2,2}^{5,4}+v_{1,2,2,2,2}^{5,5})\big|_{t=0}x_1$$ for which we need to calculate $v^{5,3}$, $v^{5,4}$: &&v\_[1,2,2,2]{}\^[5,4]{}=v\_[1,2,2]{}\^[4,3]{}(2t)+\_t v\_[1,2,2,2]{}\^[4,4]{}=12t(2t)+\_t(2t)\^3=24t\^2+6(2t)\^2=48t\^2\
&&v\_[1,2,2]{}\^[5,3]{}=v\_[1,2]{}\^[4,2]{}(2t)+\_t v\_[1,2,2]{}\^[4,3]{}=\^2\_t (2t)(2t)+\_t 12t=12and finally $$c_4(1)=(\partial_t^3(2t)+48t^2+12+(2t)^4)\big|_{t=0}x_1=12$$ thus the coefficient of $t^4$ is $12x_1/ 4! = x_1/2$.\

### The basic stationary Driver-type ODE {#sss:blons}

$\dot x=ax^2$, for $a\in \re^+$. We have $\cal S=\{1\}$, $V=a$. The solution $(x,0)$ can be easily calculated by variable separation, the result is x(t)=[x1-axt]{}.\[ex3-sol\] Let us calculate the solution through formula (\[cstat\]): all coefficients $c_k(1)$ have each only one term in formula (\[Vkstat\]) : $c_k(1)=v^{k+1}_{1,1_k}kv_{1,1}$, therefore, by successive substitutions: $$c_k(1)=v^{k+1}_{1,1_k}ka=v^k_{1,1_{k-1}}k(k-1)a^2=\ldots =k!a^k$$ replacing which in (\[cstat\]) yields $$x(t)=\sum_{k=0}^\infty c_k(1){t^k\over k!}=x\sum_{k=0}^\infty a^kx^kt^k={x\over 1-axt}$$ which is (\[ex3-sol\]). The convergence ray, $r$, is given by the ratio’s criterion: $${1\over r}=\lim_{k\to +\infty} {|a^{k+1}x^{k+1}|\over |a^kx^k|}=a|x|.$$ As expected, the convergence interval is $(-1/a|x|, 1/a|x|)$: if $x>0$ (resp: $x<0$) the solution goes to $+\infty$ on the right (resp: to $-\infty$ on the left) end point. The limit function is defined and analytic for $t\le -1/a|x|$ (resp: for $t\ge 1/a|x|$), thus, by analytic prolongation, it is the maximal solution, defined on ${\bf I}_{x,0}=(-\infty, 1/a|x|)$ (resp; on ${\bf I}_{x,0}=(-1/a|x|,+\infty)$).\

### The Airy equation {#sss:blons}

$x^{(2)}=tx$. This is a second order linear ODE for which is not known a ’closed’ form of the solution $(x,0)$ – i.e. given through a finite composition of ordinary transcendental and/or elementary functions – but the solution can be calculated by power series. The result[^11] is a function of two parameters, namely $p_1,p_2$: x(t)=p\_1+p\_2t+[13!]{}p\_1t\^3+[24!]{}p\_2t\^4+[146!]{}p\_1t\^6+[257!]{}p\_2t\^7+[1479!]{}p\_1t\^9+…\[ex4-solf\] and an application of the ratio criterion readily yields an infinite convergence ray. In order to find (\[ex4-solf\]) by Theorem \[th:serie\], we need first to quadratize the Airy equation. Let us introduce the two indeterminates $x_1,x_2$ and set $x_1=\dot x$, $x_2=x$, then the solution $(x,0)$ of the Airy equation, is the sub-solution $x_2$, of the solution $((x_1,x_2),0)$ of the linear system:

[l]{} x\_1=tx\_2\
x\_2=x\_1

\[ex4-tra\] Note that $(0,0)$ is an equilibrium point of (\[ex4-tra\]), and the zero solution $x_1\equiv x_2\equiv 0$ is the only one singular solution. Therefore, all solutions, but the trivial one, are regular, thus (cf: §\[ssc:sinreg\]) they are equal to some solution $((x_1,x_2),0)$ with $x_1\ne 0$, and $x_2\ne 0$. Let us apply the Driver formula (\[dri\]) to (\[ex4-tra\]). We have two Driver variables $Z_{1,1}={x_2\over x_1}$ and $Z_{2,1}={x_1\over x_2}$. The equation for $Z_{1,1}$ is directly derived by the associated monomial $x_2=x_1^0x_2^1$ from which (using the exponents as coefficients): Z\_[1,1]{}=(\^1\_[1,1]{}v\_1’Z\_1 +\^1\_[1,2]{} v\_2’Z\_2)Z\_[1,1]{}=(-tZ\_[1,1]{}+Z\_[2,1]{})Z\_[1,1]{}\[eqZ11\] The equation for $Z_{2,1}$, in place of the Driver formula, can be directly derived observing that, by definition, $Z_{2,1}=Z^{-1}_{1,1}$, thus, it is the component $W_{1,1}$ of the inverse Driver (cf §\[Inv-driv\]): Z\_[2,1]{}=(tZ\_[1,1]{}-Z\_[2,1]{})Z\_[2,1]{}.\[eqZ21\] By renaming $x_3=Z_{1,1}$, $x_4=Z_{2,1}$, and gathering (\[ex4-tra\]) – rewritten as a final ODE – (\[eqZ11\]), and (\[eqZ21\]) we obtain the (inclusive) quadratization:

[l]{} x\_1=tx\_3x\_1\
x\_2=x\_4x\_2\
x\_3=x\_3x\_2-tx\_3\^2\
x\_4=tx\_4x\_3-x\_4\^2

\[ex4-tra2\] whose frame is V=(

[cccc]{} 0&0&t &0\
0&0&0&1\
0&1&-t&0\
0&0&t&-1

) \[ex4-fr\] Thus $\cal S=\{2,3,4\}$. Now, let us verify that $x_2(t)$ calculated through (\[TE\]) agrees with (\[ex4-solf\]). We have $c_0(2)=x_2$, which is the zero-coefficient of (\[ex4-solf\]) for $p_1=x_2$. As for $c_1(2)$: $$c_1(2)=\sum_{i_1\in \{2,3,4\}} v^{2,2}_{2,i_1} x_2x_{i_1}=v^{2,2}_{2,2} x^2_2 +v^{2,2}_{2,3} x_2x_{3} +v^{2,2}_{2,4} x_2x_{4}$$ thus, recalling that $x_3={x_2\over x_1}=x_4^{-1}$, and observing, in (\[ex4-fr\]), that $v_{2,2}=v_{2,3}=0$, $v_{2,4}=1$, we have: $$c_1(2)=x_1$$ which agree with (\[ex4-solf\]) provided that $x_2=p_1$. Now, the coefficient of $t^2$ in (\[ex4-solf\]) is zero, so we expect to obtain $c_2(2)|_{t=0}=0$. As a matter of fact:

[rl]{} c\_2(2)&=\_[i\_1{2,3,4}]{} v\^[3,2]{}\_[2,i\_1]{}x\_2x\_[i\_1]{} +\_[i\_1,i\_2{2,3,4}]{} v\^[3,3]{}\_[2,i\_1,i\_2]{} x\_2x\_[i\_1]{}x\_[i\_2]{}\
&=v\^[3,2]{}\_[2,2]{}x\_2\^2+v\^[3,2]{}\_[2,3]{}x\_2x\_[3]{} +v\^[3,2]{}\_[2,4]{}x\_2x\_[4]{}\
&+ v\^[3,3]{}\_[2,2,2]{} x\_2\^3 +v\^[3,3]{}\_[2,2,3]{} x\_2\^2x\_[3]{}+v\^[3,3]{}\_[2,2,4]{} x\_2\^2x\_[4]{}\
&+ v\^[3,3]{}\_[2,3,2]{} x\^2\_2x\_3 +v\^[3,3]{}\_[2,3,3]{} x\_2x\^2\_[3]{}+v\^[3,3]{}\_[2,3,4]{} x\_2x\_3x\_[4]{}\
&+ v\^[3,3]{}\_[2,4,2]{} x\_2\^2x\_4 +v\^[3,3]{}\_[2,4,3]{} x\_2x\_4 x\_[3]{}+v\^[3,3]{}\_[2,4,4]{} x\_2x\^2\_[4]{}

\[ex4-tm\] from (\[Vk2\]): $v^{3,2}_{2,i_1}=\partial_t v_{2,i_1}=0$, for $i_1=2,3,4$, and from (\[Vkk\]) we have

[l]{} v\^[3,3]{}\_[2,2,i\_2]{}=v\_[2,2]{} (v\_[2,i\_2]{}+v\_[2,i\_2]{})=0, i\_2=2,3,4\
v\^[3,3]{}\_[2,3,i\_2]{}=v\_[2,3]{} (v\_[2,i\_2]{}+v\_[3,i\_2]{})=0 i\_2=2,3,4\
v\^[3,3]{}\_[2,4,2]{}=v\_[2,4]{} (v\_[2,2]{}+v\_[4,2]{})=0\
v\^[3,3]{}\_[2,4,3]{}=v\_[2,4]{} (v\_[2,3]{}+v\_[4,3]{})=t\
v\^[3,3]{}\_[2,4,4]{}=v\_[2,4]{} (v\_[2,4]{}+v\_[4,4]{})=1-1=0

\[ex4-tm2\] therefore $
c_2(2)=v^{3,3}_{2,4,3} x_2x_{4}x_{3}=t x_2
$ and $c_2(2)|_{t=0}=0$ as expected. Clearly, the more $k$ increases the more difficult is carrying out the calculation by hand for $c_k(2)$. It’s still meaningful though, and not too cumbersome, to perform the calculation of one more coefficient, $c_k(2)$, as the parameters $p_1,p_2$ have been fixed, and we expect $c_3(2)|_{t=0}=x_2$. We have c\_3(2)=\_[i\_1{2,3,4}]{} v\^[4,2]{}\_[2,i\_1]{}x\_2x\_[i\_1]{} +\_[i\_1,i\_2{2,3,4}]{} v\^[4,3]{}\_[2,i\_1,i\_2]{} x\_2x\_[i\_1]{}x\_[i\_2]{} +\_[i\_1,i\_2,i\_3{2,3,4}]{} v\^[4,4]{}\_[2,i\_1,i\_2,i\_3]{} x\_2x\_[i\_1]{}x\_[i\_2]{}x\_[i\_3]{} \[ex4-c3\] and by (\[Vkk\]): $v^{4,4}_{2,i_1,i_2,i_3}=v^{3,3}_{2,i_1,i_2}(v_{2,i_3} +v_{i_1,i_3}+v_{i_2,i_3})$. By (\[ex4-tm2\]) $v^{3,3}_{2,4,3}=t$ is the only non zero element in the vector $v^{3,3}_2$, moreover, $v^{4,2}_{2,i_1}=\partial^2_tv_{2,i_1}=0$, and since $v^{3,2}_2=0$, by (\[Vks\]) we have $v^{4,3}_2=\partial_t v^{3,3}_{2}$, which is zero but $v^{4,3}_{2,4,3}=1$. Therefore (\[ex4-c3\]) turns into $$c_3(2)=v^{4,3}_{2,4,3} x_2x_{4}x_{3}
+t\sum_{i_3\in \{2,3,4\}} (v_{2,i_3} +v_{4,i_3}+v_{3,i_3}) x_2x_{4}x_{3}x_{i_3}\\$$ thus $c_3(2)|_{t=0}=v^{4,3}_{2,4,3} x_2x_{4}x_{3}=x_2$, as expected. Worth to be stressed is the fact that, the method of coefficients calculation we have pointed out in the latest note, applies just [*after*]{} having fixed the initial time-point ($t_0$, which has been set to zero), whereas, through the new method here presented we calculate the general expression, as a function of $t$, of $c_k(2)$ first, and then set $t=t_0=0$. Therefore we can use the same formulas for calculating the Taylor expansion with any different initial point $t_0$.\

Connection with the 22nd Hilbert’s Problem\
and Further Developments {#sc:HPFD}
===========================================

A partial solution to a differential form of the 22rd Hilbert’s problem
-----------------------------------------------------------------------

The exact quadratization-based solution method described in the present paper can be interpreted as a kind of [*uniformization*]{} of a sub-class of analytic differential equation and as such it represents a (partial) solution, in the real case, of a [*differential version*]{} of the 22nd Hilbert’s problem [@Hi12], as below explained. In the original formulation, the 22rd Hilbert’s problem consists, given a set of complex polynomial (non differential) equations \^nF(x,v)=0,\[Fvx\] in $n$ complex indeterminates, $x\in \co^n$, with $m$ complex [*parameters*]{} (say $v$ the vector in $\co^m$ collecting them) in finding an open dense set in $\cal V\subset \co^m$ and a surjective analytic function $f:\cal V\to \co^n$, such that $f(v)$ solves (\[Fvx\]), i.e. $(f(v),v)$ describes the set of solutions of (\[Fvx\]) through $v$ (which means that this set is ’uniformized’). In the most general formulation, $x$ and $v$ are both ’indeterminates’: it is the uniformization $(f(v),v)$ that actually decides that $v$ is a set of ’parameters’. In other words: the uniformization is a parametrization of the solutions set such that (by the surjectivity of $f$) for any $x$ there is always a set of parameters $v$ such that $x$ solves (\[Fvx\]). We do not consider here this more general formulation, since throughout this paper we have adopted, the more general [*differential*]{} setting on the one hand, but a more restrictive one on the other hand, in which all quantities are real, but we always know in advance what is a parameter and what an indeterminate, and that an unique solution always exists whatever we choose a set of parameters. A complete solution of the (non differential) uniformization problem has so far not been obtained, except for the one dimensional case, whereas the solution proposed by Poincairé in 1907 [@Po07] is considered a partial one.

Now, provided we move to the real analytic case, and generalize the equations by allowing them 1) to be ODEs: \^nG(w,x,x,t)=0,\[Fvxt\] depending on a set $w$ of real parameters (which can be even analytic functions of $t$) 2) for any fixed $w$, $G$ to belong to the class of analytic $\sigma\pi$-reducible functions (thus, properly including the polynomials), of rank $n$ with respect to the $\dot x$ indeterminates, then Theorems \[th:serie\], and \[th:seriest\], provide a solution of this generalized uniformization problem. For, note first that for any given $w$, a solution $(x,t_0)$ of (\[Fvxt\]), which is a real analytic function defined on some open interval ${\bf I}_{x,t_0}\subset \re$, can be identified with the sequence of coefficients of its Taylor expansion: $c=\{c_k\}_{k=0}^\infty$, where $c_k=(c_k(1)(t_0),\ldots, c_k(n)(t_0))$. Let ${\cal S}_{t_0}$ the set of all sequences $c$ such that the corresponding Taylor series centered at $t_0$ has a convergence ray $r>0$. Now, as we have shown in §\[sc:imp\], (\[Fvxt\]) is quadratizable, thus $(x,t_0)$, for fixed $w$ is a sub-solution of some Driver in $d$ indeterminates ($d\ge n$), and such a Driver is characterized by a frame $V(w)\in \cal V$ – where $\cal V$ denotes the set of all $d\times d$ matrices – which is a function (analytic) of the original coefficients $w$. Theorem \[th:serie\] defines an analytic map $\cal V\to \cal S_{t_0}$, giving the solution in the form of a sequence $c\in \cal V$. Therefore the quadratizing map combined with the recursive equations (\[Vkk\])–(\[Vk2\]) defines a map $f:w\mapsto c$, such that $c=f(w)$ is the solution $(x,t_0)$ of (\[Fvxt\]). The latter is a kind of generalized (to the differential case) uniformization for the ODE (\[Fvxt\]).

Such [*differential*]{} uniformization, clearly does not imply that the problem is solved as well in the non differential case, but surely that, if an uniformization exists for the set of solutions of the real analytic equation (\[Fvx\]), for $F$ in the class of differentially $\sigma\pi$-reducible functions (i.e. such that the associated ODE is $\sigma\pi$-reducible), then the same uniformization is the differential uniformization of the ODE (\[Fvxt\]), for $G(v,x,0,t)=F(x,v)$, thus of the set of all [*constant solutions*]{} of the ODE, i.e. of its [*equilibrium points*]{}. By the theory developed in the previous section, the latter problem amounts to finding the $x$, if any, such that $c_k=0$ $\forall k>0$ in (\[Vkk\])–(\[Vk2\]), so that $c_0=x$ and $x(t)\equiv x$, where $x(\cdot)$ is the curve of the solution $(x,t_0)$ of the ODE.\

Control problems, and differential flatness
-------------------------------------------

We sketch here few guidelines for a possible development of the present research, by showing that the results of §\[sc:SDTheo\] can be used for finding a [*flat output*]{} for a control system. [*Differential flatness*]{} for control systems has been first introduced in [@FLMR92], [@FLMR93], in the context of differential algebra, and more recently defined in a more geometric Lie-Bäcklund framework relying on the paper [@Po93]. Most of the latter approach can be found in the textbook [@Le09]. We point out as well a further approach in which flatness can be described in terms of the notion of absolute equivalence defined by E. Cartan in [@Ca14], [@Ca15b], and where a nonlinear control system is viewed as a Pfaffian system on an appropriate space [@NiRM98]. The discussion below, mirrors the Lie-Bäcklund framework of [@Le09], but it is translated in the setting of the present paper.

Let $y=y^{(0)}$, and $y^{(k)}$, $k$ integer, be a countable set of indeterminates in $\re^n$, and consider the following [*countably infinite*]{} ODEs: \_t y\^[(k)]{}=y\^[(k+1)]{}, k.\[inf-lin\] We denote by $({\bf y},t_0)$, where $\bf y$ is the sequence $\{y^{(k)}\}_{k\in \na}$, the generic solution of (\[inf-lin\]) passing through $t_0\in \cal T\subset \re$, $\cal T$ an open interval, i.e. the function $t\mapsto {\bf y}(t)$, such that $y^{(k)}(t_0)=y^{(k)}$. Note that [*any*]{} map $t\mapsto {\bf y}(t)$, which is analytic on $\cal T$ in all components, determines a solution of (\[inf-lin\]), where $y^{(k)}$ is the $k$th derivative of $y$ for any k. By reason of that (\[inf-lin\]) is in fact said [*the trivial system*]{}. Now suppose that the ODE x=f(x,t),\[Ofl\] is a Driver, with domain $\cal D\times \cal T\subset \re^{n+1}$, of some pre-Driver (i.e. it is a quadratization of some $\sigma\pi$-system, cf. §\[SD:sec\]). In the most general situation, in turn, the latter $\sigma\pi$-system includes as sub-solutions all solutions of some original $\sigma\pi$-reducible system. Let $v_{i,j}$ the $n^2$ coefficients of (\[Ofl\]), all defined and analytic on $\cal T$, and suppose there exist analytic functions $h_{i,j}: \cal U{\rm (open)}\subset \re^m \to \re^{n^2}$, with $m\le n^2$, such that $v_{i,j}=h_{i,j}(u)$, where $u\in \cal U$, is a vector of [*controls*]{}, analytic time-functions defined on $\cal T$, that can be arbitrarily chosen.

A control problem can be set as follows, given a function $t\mapsto x(t)$, find $u(t)$ such that $x(t)$ is the solution of (\[Ofl\]) such that $x(t_0)=x$, for some $x\in \cal D$ and $t_0\in \cal T$. [^12] Now, let us rewrite (\[Ofl\]) by defining $g(x,u)=f(x,t)$, which makes explicit the dependence on $u$, define the new indeterminates $u^{(k)}$, $k\in \na$, with $u^{(0)}=u$, and consider the infinite dimensional ODE:

[l]{} x=g(x,u),\
\_t u=u\^[(1)]{},\
\_t u\^[(1)]{}=u\^[(2)]{}\
\

\[ext-sy\] whose solutions $((x,{\bf u}), t_0)$ – $\bf u$ denoting the sequence $\{u^{(k)}\}_{k\in \na}$ – include the solution $(x,t_0)$ of (\[Ofl\]) for $v_{i,j}=h_{i,j}(u)$ (and in fact it is the sub-solution of (\[ext-sy\]) in the first $n$ components). We denote $\cal D^*$ (resp: $\cal Y$) the domain of (\[ext-sy\]) (resp: of (\[inf-lin\])) given by[^13]: &&D\^\*=DU\^m\_,\[domDst\]\
&&Y=\^n\_.\[domcalY\] Note that, by the recursive equations (\[Vkk\])–(\[Vk2\]), there exist analytic functions $\phi^{k,s}$ such that v\^[k,s]{}=\^[k,s]{}(u,u\^[(1)]{},…,u\^[(k)]{}), s=2,…, k; k2\[vksu\] then Lemma \[lm:serie\] and Theorem \[th:serie\] prove that there exists a map : ((x,[**u**]{}), t\_0)([**y**]{}, t\_0),\[mapPsi\] such that $x=y$. As a matter of fact, such map is given by $y^{(k)}=q^{k+1}(x)$, $k\ge 0$, where $q^1(x)=x$, and for $k\ge 2$, the $q^k$’s are the polynomials given recursively by (\[qkrec\]). Then, formula (\[qkrec\]), the definition of $A_i$ given in (\[dA\]), and the differential formula (\[papk\]), shows that $p^{k+1}$ is a polynomial in $x$ and in $v, \partial_t v, \ldots \partial_t^k v$. Therefore, by (\[vksu\]) for any $k\in \na$ there exist, [*and can be calculated*]{}, analytic functions $\psi^k$, such that y\^[(k)]{}=\^k(x,u,u\^[(1)]{},…,u\^[(k)]{}) \[ykpsi\] and the map $\Psi$ in (\[mapPsi\]) is given by the infinite aggregate of the maps (\[ykpsi\]). The control problem is solved, as soon as one is able to show that the map (\[mapPsi\]) is invertible, and, in particular, it is a diffeomorphism from $\cal D^*$ onto its own image in $\cal Y$. Here we are assuming a setting similar as in [@Le09], and we consider $\cal D^*$ and $\cal Y$ as differentiable manifolds (in infinite dimension), with the underlying topology generated by sets of the type: $U\times \re^m_\infty$ – for the manifold (\[domDst\]) – with $U$ open set in $\cal D\times (\re^m)^N$ (for some positive integer $N$), and $U\times \re^n_\infty$ – for the manifold (\[domcalY\]) – with $U$ open set in $\cal D\times (\re^n)^N$. A diffeomorphism $\Phi$ between manifolds endowed with the topology defined above, is said a [Lie-Bäcklund isomorphism]{}: it has components depending on a finite set of arguments only [@Le09], thus, there exist analytic maps $\alpha,\beta$, and an integer $\rho$, such that &&x=(y,y\^[(1)]{}, …, y\^[()]{})\[x-alp\]\
&&u=(y,y\^[(1)]{}, …, y\^[()]{})\[x-bet\] Note that (\[x-alp\]) in our setting is always verified by $x=y$, which is given by the inverse of the first $n$ components of the map $\Psi$ in (\[mapPsi\]) (which is the identity in $\re^n$), thus, the solution to the control problem in open loop, amounts to finding $\beta$ in (\[x-bet\]) only. Generally speaking, Theorem \[th:serie\] suggests a principle way for calculating $\beta$: noticing that $c_k(i)=y_i^{(k)}$, one has to, first, find an integer $\rho$, and $\rho$ equations among (\[GE0\]), (\[GEs\]), that are solvable with respect to a finite number, upper bounded by $\rho$, of $v^{k,s}$, and next in solving – on account of (\[vksu\]) – a finite set of differential equations, extracted from (\[Vkk\])–(\[Vk2\]), with respect to $u$. It should be noted that, a necessary condition for the map $\Phi$ to be a Lie-Bäcklund isomorphism, is that the flat output $y$ has the same dimension of the input $u$ [@Le09]. This is verified in our scheme, where $x=y$, only in a situation (not usual for control systems) where $m=n$. As a matter of fact, since (\[Ofl\]) is a self-Driver, the [*original*]{} system has a lower number of indeterminates, say $l\le n$, and since (\[Ofl\]) is a quadratization of the original system, there are just $l$ out of the $n$ $x_i$’s which we are interested to, and typically $m\le l$ (the number of controls is lower than the number of states). Nevertheless, Lemma \[lm:serie\] holds [*component-wise*]{}, therefore we can even set $y\in \re^m$, and $y^{(k)}_{i_s}=q_{i_s}^{k+1}$, and obtain a trivial system as well, for a suitable choice of $i_1,\ldots,i_m$ such that $y^{(k)}_{i_s}=x^{(k)}_{i_s}$, $\forall k\in \na$, and $x_{i_s}$, $s=1,\ldots, m$ are $m$ out of the $l$ indeterminates associated to the solution of the original system. By using this setting, the jet-manifold (\[domcalY\]), i.e. the domain of the trivial-system, becomes $\cal Y=\re^m_{\infty}$, the necessary condition for a Lie-Bäcklund isomorphism can be verified, $\alpha$ in (\[x-alp\]) is no more an identity, but the problem is well-posed in any realistic situation.

We highlight that, as noted in [@Le09], flatness can be seen as a solution of another kind of generalized (in the framework of manifolds of jets of infinite order) uniformization problem of the 22rd Hilbert’s problem.

Conclusion
==========

The goal of this paper was twofold: first, to give new results (generalizations and/or improvements) and further insight on the topic of exact quadratization, first issued in the article [@Ca15], and second, which is the main result, to show that exact quadratization allows to write the solution of any analytic quadratizable ODE (roughly speaking: [*practically all*]{} analytic ODEs) as a power series, directly from the ODE coefficients.

Sections I and II has been devoted to the first task. In §I we have summarized the basic result of [@Ca15], but using a quite different terminology, which in our intention aimed to focus the [*strictly essential*]{} features of [*quadratizable*]{} ODEs. In §\[ssc:sinreg\] we have given the definitions, of singular and regular part of an analytic system. which are used to prove one of the two most important new results of the section, which is Theorem \[th:GEQ\] showing the globality of the exact quadratization of $\sigma\pi$-systems. Theorem \[th:GEQ\] relies on Theorem \[th:EQns\], which is a re-formulation, in the setting of the present paper – and where also a shorter proof is offered – of the basic result on exact quadratization for $\sigma\pi$-systems of [@Ca15]. The second most relevant result of the section is given in Theorem \[th:imp\], which is the generalization of exact quadratization to implicitly defined $\sigma\pi$-ODEs. A lot of top-up issues on exact quadratization – that were not given or not enough scrutinized in [@Ca15] – are given throughout the rest of the section, such as: definitions and properties of the [*inverse Driver*]{} (§\[Inv-driv\]) which allows to view exact quadratization as a multi-variable generalization of the [*Bernoulli-type*]{} ODEs ([*Ibid.*]{}), as well as the concept of inclusive quadratization (Definition \[df:inc-qua\]), Driver-type ODE, and self-Driver (given in §\[SD:sec\]),

In §III we have given the main result of the paper, stated in Theorem \[th:serie\]. The general recursive formula (\[Vkk\])–(\[Vk2\]) allows the calculation, directly from the coefficients $v_{i,j}$, of the coefficients $c_k(i)$, for any $k\in\na$, of the Taylor expansion (\[SPSD\]) of the $i$th component $x_i(t)$ of the solution of any [*Driver-type*]{} ODE, and, hence, by Theorem \[th:Drty\], and by the exact quadratization theorem, of the solution of any $\sigma\pi$-reducible ODE. Theorem \[th:seriest\] yields the stationary case, in which the recursive formula for the coefficients calculation, simplifies in (\[Vkstat\]). Theorem \[th:serie\].is proven on the basis of a preliminary result, given in Lemma \[lm:serie\], which has an importance in its own right, since formula (\[kder\]) allows to write any $\sigma\pi$-reducible ODE as a [*trivial system*]{}, a fact that we have highlighted in the final section, §\[sc:HPFD\], where the interplay of exact quadratization with the 22rd Hilbert’s problem has been discussed, as well as with the concept of [*differential flatness*]{}, used in control theory.

[14]{} \[1\][\#1]{} \[1\][`#1`]{} urlstyle \[1\][doi:\#1]{}

F. Carravetta,“Global Exact Quadratization for Continuous-Time Nonlinear Control Systems” [*SIAM Journal on Control and Optimization*]{}, vol. 53, pp.  235 - 261, 2015.

W. Boothby. An introduction to Differentiable Manifolds and Riemannian Geometry. Academic Press, New York, 1975.

M. Fliess, I. Kupka “A finiteness Criterion for Nonlinear Input-Output Differential Systems”. [*SIAM Journal on Control and Optimization*]{} vol. 21, pp.  721 - 728, 1983.

H. Cartan. Calcul Différentiel. Hermann, Paris, 1967.

E. Cartan. “Sur l’equivalence absolue de certains systèmes d’equations différentielles et sur certaines familles de courves”. [*Bull. Soc. Math. France*]{}, vol. 42, pp. 12–48, 1914. Also in [*Ouvres Complètes*]{}, part II, vol. 2, pp.1133–1168, CNRS, Paris, 1984.

E. Cartan. Sur l’intégration de certains systèmes indéterminés d’equations différentielles. [*J. für reine und angew. Math.*]{}, vol. 145: pp. 86–91, 1915. Also in [*Ouvres Complètes*]{}, part II, vol. 2, pp.1164–1174, CNRS, Paris, 1984.

D. Hilbert. “Uber den Begriff der Klasse von Differentialgleichungen” [*Math. Ann.*]{}, Vol 73, pp. 95 – 108, 1912. Also in [*Gesammelte Abhandlungen*]{}, Vol. 3, pp. 81 – 93, Chelsea, New York, 1965.

H. Poincaré. “Sur l’uniformisation des fonctions analytiques”, [*Acta Mathematica*]{}, Vol 31, pp. 1 – 63, 1907. Also in Oeuvres de Henri Poincaré, t. 4, pp. 70 – 139, Gauthier-Villars, Paris, 1950.

L. Bers, “On Hilbert’s 22nd problem”. In Browder F., editor, Mathematical Developments Arising From Hilbert Problems, [*Proceedings of Symposia in Pure Mathematics*]{}, pp. 559 – 609. American Mathematical Society, Providence, Rhode Island, 1976.

M. Fliess, J. Lévine, Ph. Martin, and P. Rouchon. “Sur les systèmes non linéaires différentiellement plats” [*C.R. Acad. Sci. Paris*]{}, I-315, pp. 619–624, 1992.

M. Fliess, J. Lévine, Ph. Martin, and P. Rouchon. “Linéarisation par bouclage dynamique et transformations de Lie-Bäcklund”. [*C.R. Acad. Sci. Paris*]{}, I -317 pp. 981– 986, 1993.

J.B. Pomet. “A differential geometric setting for dynamic equivalence and dynamic linearization”. In [*Workshop on Geometry in Nonlinear Control*]{}, Banach Center Publications, Warsaw, 1993.

J. Lévine “Analysis and Control of Nonlinear Systems: A Flatness-based Approach”, Springer, 2009.

M. van Nieuwstadt, M. Rathinam, and R. M. Murray. “Differential flatness and absolute equivalence”. [*SIAM Journal on Control and Optimization*]{}, vol. 36, pp. 1225–1239, 1998.

[^1]: We state in a moment in which sense we talk about a [*sub-solution*]{}.

[^2]: For any pair of points in $\cal A$ we can say they are equivalent if they are connected by the graph of a solution. Indeed, one readily verifies that this is an equivalence relation. This motivates the identification of the graph (and, hence, of the solution) with any point $(x,t)$ of the graph itself.

[^3]: In other words $\Psi$ preserves (locally) the equivalence relation between points of the same solution, i.e. maps (locally) solution of (\[S1\]) into solution of (\[Quad1\]). In a differential geometric framework (see for instance [@Bo75]) that is the well known property of [*relatedness*]{} of the vector fields associated to the two ODEs.

[^4]: In §II, it will be made clear that $\Phi$ is meromorphic on $\cal A$, but analytic on $\cal A\setminus \{x\in\re^n | x_i=0, i=1,\ldots,n\}={\rm Dom}\{\Phi\}$

[^5]: In [@Ca15], exact quadratization has been defined as a kind of [*systems immersion*]{}, a concept defined in [@FlKu83] in a framework of control systems

[^6]: Note that ${\bf S}={\bf C}$ if and only if all criticality indices are singularity indices

[^7]: Notice that the [*full*]{} system in the indeterminates $x_1,x_3$ associated to (\[ex-singpart\]), has the solution $x_1=0$, and thus a singular part in $x_3$, given by $\dot x_3=4+x_3^{1\over 2}$. However $x_1=0$ is not a sub-solution of the sub-system ${\bf S}\times I_k$.

[^8]: Recall that conventionally the zero-system has domain $\re^n$.

[^9]: We chose, in the paper [@Ca15], the names [*Driver*]{} and [*final stage*]{}, since a canonical quadratization can be depicted, in a similarity with electrical devices, as a [*cascade*]{} of systems, where the Driver actually ’drives’ the final stage yielding the signal representing the original ODE solution.

[^10]: We point out that the driver equations (\[dri\]) are quadratic differential equation of the Bernoulli type (Riccati type with no zero degree term). The ODE (\[inverse-driver\]) is in fact consistent with the basic properties of this kind of equations.

[^11]: A way for calculating the solution, which can be found in many exercise books of analysis courses, is to set it equal to $\sum_{k=0}^\infty a_kt^k$ (the general form of a power series), and then calculating the second derivative (term to term). By substitution, equating the terms with the same degree, we obtain a recursive equation for the coefficients $a_{k+1}=a_{k-1} (k+1)^{-1}(k+2)^{-1}$, relating every coefficient to the coefficient placed two times before, which allows the calculation of all coefficients from $a_0$ and $a_1$. By setting $p_1=a_0$, $p_2=a_1$ we obtain (\[ex4-solf\])

[^12]: We use here some abuse of terminology, in that what we call a ’control problem’ is actually a kind of [*inversion problem*]{}, whereas what is actually of interest in control theory is the [*tracking*]{} of a trajectory $t\to x(t)$. Moreover, we are considering just an [*open loop*]{} perspective.

[^13]: $\re^n_{\infty}$ denotes the infinite product $\re^n\times \re^n \times \ldots$
