> <\body> |~>>|~>>|~>|~>>|~>>>||algebraic differential equations>||<\author-address> Dépt. de Mathématiques (bât. 425)Université Paris-Sud91405 Orsay CEDEXFrance >|>>> <\abstract> In our PhD. we have given an algorithm for the algebraic resolution of algebraic differential equations with real transseries coefficients. Unfortunately, not all equations do admit solutions in this strongly monotonic setting, even though we recently proved an intermediate value theorem. In this paper we show that the algorithm from our PhD. generalizes to the setting of weakly oscillatory or complex transseries. Modulo a finite number of case separations, we show how to determine the solutions of an arbitrary algebraic differential equation over the complex transseries. We will show that such equations always admit complex transseries solutions. However, the field of complex transseries is not differentially algebraically closed. In , we have studied the asymptotic behaviour of solutions to algebraic differential equations in the setting of strongly monotonic or real transseries. We have given a theoretical algorithm to find all such solutions, which is actually effective for suitable subclasses of transseries. More recently, we have proved the following ``differential intermediate value theorem''. > Let be the real field of grid-based transseries in and let be a differential polynomial with coefficients in . Then, given transseries g\\> with 0> and 0>, there exists a \> with h\g> and .> This theorem implies in particular that any algebraic differential equation of odd degree, such as +e>*f*f+\(log x+1)*f+log (e+\(\(x))=0,> has at least one real transseries solution. This theorem is striking in the sense that it suggests the existence of theories of ordered and/or valuated differential algebra.\ However, a main drawback of the setting of real transseries, is that not every algebraic differential equation can be solved; actually, even an equation like +1=0> has no solutions. In order to get a better understanding of the asymptotic behaviour of solutions to algebraic differential equations, it is therefore necessary to search for a complex analogue of the theory of real transseries. This paper is a first contribution in this direction. The first problem is to actually define complex transseries. The difficulty is that it is not clear whether an expression like >> should be seen as an infinitely large or an infinitely small transmonomial. Several approaches can be followed. A first approach, based on pointwise algebras, was already described in chapter 6 of . However, this approach has the drawback that it is not easy to compute with complex transseries. A second more computational approach is described in section . Roughly speaking, it is based on the observation that all computations with complex transseries can be done in a similar way as in the real setting, except for testing whether a monomial like >> is infinitely large or small. Now whenever we have to make such a choice, we will actually consider both cases, by applying the automatic case separation strategy (see ). We implicitly reject the case when >> is bounded, which is ``degenerate'', but which deserves to be studied later. The last approach, which is described in section , is more structural and really allows us to define a complex transseries in a not too difficult way. The underlying idea is analogue to the concept of a maximal ideal. Intuitively speaking, we assume the existence of some ``god'', who has decided for us which monomials like >> are infinitely large and which ones are infinitely small. It turns out that all possible choices lead to isomorphic fields of transseries. However, the geometric significance of these fields is hard to grasp. In section , we introduce parameterized complex transseries, which are necessary to express generic solutions to differential equations. Indeed, such solutions may involve integration constants. As usual, our approach is based on the automatic case separation strategy. The remaining sections deal with the resolution of asymptotic algebraic differential equations with complex transseries coefficients. Our approach is similar to the one followed in , but we have made a few simplifications and we corrected an error (see section ). Our main results are stated in sections and . We show that there exists a theoretical algorithm to express the generic solution to an algebraic differential equation by means of parameterized complex transseries and we give a bound for the logarithmic depth of the generic solution. We also show that an algebraic differential equation of degree admits at least complex transseries solutions when counting with multiplicities. As a consequence, each linear differential equation admits a full system of solutions. However, our fields of complex transseries are not differentially algebraically closed and several interesting problems still need to be solved (see section ). The reader should be aware of a few changes in notations , which are summarized in the following table: >|>|>|>|>|>|>|>|>|>|>|>|>|>|>|>|>|>>>>>> In all what follows, let be a and its complexification. This means that has the structure of a totally ordered field and functions R>, which are compatible with this ordering. More precisely, we assume that admits an inverse function with domain >> and that the function restricted to /2,\/2[> admits a totally defined inverse. Here , where /2)> and =4*arctan 1>. Furthermore, ||>|||>>>> for all R>. Finally, for each \> resp.\>>> and R>, we require that |>|>*x+\+>*x;>>||>|>*x+\>*x->*x.>>>>> form a real trigonometric field.> <\proof> The functional equations are classical. The inequality for was first proved in . As to the inequality for , we have )|(2*k)!>*x=n>|(4*k)!>*1-|(4*k+1)*(4*k+2)>\0,> if 4*n>. Otherwise, )|(2*k)!>*x=1-|2>+|(4*k)!>*1-|(4*k+1)*(4*k+2)>\-1\cos x,> since 1>. and for expansions at order instead of .> and may naturally be extended to , numbers in may naturally be written in polar form, etc.> and are no longer required to be total and the functional equations resp. inequalities are only required to hold, whenever they make sense. For instance, if dom exp>, then we require that dom exp> and .> Let > be a totally ordered monomial group (or set) with -powers. Then we recall that the field >> of grid-based power series is naturally totally ordered by 0\c\0>, for all 0>. This ordering is compatible with the multiplication: 0\g\0\f*g\0>. More generally, if > is only partially ordered, then we define an ordering on >> to be with the asymptotic ordering on >, if g\g\0\f+g\0> for all R>>. In what follows, we are rather interested in the complexification >> of >>. Obviously, this -algebra can not be given an ordering which is compatible with the multiplication. Nevertheless, it is interesting to consider orderings on >> which are only compatible with the -algebra structure of >>. Such an ordering is again said to be compatible with the asymptotic ordering on >, if () holds for all R>>. Assuming that such orderings on > and >> are total, the condition () implies that =\\sign f=sign g> for all non zero C>>. Consequently, the ordering on >> is totally determined by the sets >={c\C\|c*\\0},> where > runs over >. Each >> is actually the set of strictly positive elements of a total ordering on , which is compatible with the -module structure of . Therefore, each >> is characterized by an angle >\[>,\)> and a direction >\{,1}>, via >={c\C\|(Re (c*e>>)\0)\(Re (c*e>>)=0\Im(\>*c*e>>)\0)}.> This situation is illustrated in figure . |Shape of the region >> for >=1> resp. >=>.> In is also possible to consider complex powers of monomials: a is a monomial group > with -powers, with an asymptotic ordering > which is compatible with the expo-linear -vector space structure of >. For instance the formal group >*z>> is a monomial group with -powers for the ordering >*z>\1\(Re \\0)\(Re \=0\Re \\0)>. This group is not totally ordered, since > and are incomparable. We may make the ordering total by deciding that >*i>\(log z)>*i>\1>. Now consider a totally ordered grid-based algebra of the form =C>>, where > is a totally ordered complex monomial group, and where the ordering on is assumed to be compatible with the asymptotic ordering on >. Assume that we also have a partial logarithmic function on , such that <\description> coincides with the usual logarithm on >>; If dom log>, then dom log> and . We say that is a if the following conditions are satisfied: <\description> \>\|c\R>}>; \\>>, for all \\>; For all \\>>, we have )=l\\>, where >)|k>\C[[z]]>, as well as the following conditions for the logarithm: <\description> >=\*log \> for all \\> and \C>; \\\log \\log \> for all ,\\\>; \\> for all \\\\{1}>. In , we write g>, if and only if \log \> for all \>>. In view of , this means that g\f>\g>> for all ,\\C>>. until replace the requirement that should be a logarithmic function in the definition of real fields of transseries. It can indeed be checked that our conditions are to the usual conditions on in this case.> <\remark> The domain of the logarithm on may be further extended, by setting +log (f/c)> for all \>>, where > is defined using the function on . Of course, such an extension of the logarithm to involves a choice of a principal determination. Furthermore, such an extension cannot satisfy both the properties and . On the other hand, the partial inverse of may be extended canonically in such a way that the equation admits a solution for each \>>. Indeed, it suffices to extend via exp(f)*exp(c)> for all dom exp> and C>. In what follows, we will always assume that the partial inverse of has been extended in this way. > Consider the formal -vector space =C*log z\C*log log z\\> generated by the formal symbols >. Given angles ,\,\\[-\,\)> and directions ,\,\\{,1}>, we define a total ordering on > as explained in section . Then the formal exponential =z*(log z)*\> of > is a complex monomial group for the asymptotic ordering > defined by \\\log \\log \>, where >*(log z)>*\)=\*log z+\*log log z+\>. In order to avoid confusion, we will sometimes write ,\>> instead of >. Assume from now on that > and > were chosen such that >*log z\R>*log log z\\\0>. Given a non-zero grid-based series \=C>> with \R>>, we define its logarithm by +log \+l\ (1+\).> We may extend the total ordering on \\> to in a similar way as in section , by extending the angle and direction families > resp. > into larger families |^>> resp. |^>>. It is easily verified that the field |^>>,|^>>>\\> with this ordering is a pre-field of complex transseries. Actually, the structure of |^>>,|^>>>> does not really depend on the choices of >, >, |^>> and |^>>, modulo rotations and conjugations. Indeed, assume that > and > are a second family of angles and directions with indices in }>. Then we define an increasing isomorphism > between ,\>> and ,\>> by :\{log z, log log z,\}>f>*\\\{log z, log log z,\}>*e>>*\>*\>>f>*e>>*\,> where (z)>||>|(z)>||>,>>>>> for all C>. We infer that |\>:\,\>\\,\>; exp f\exp \(f)> is an isomorphism of complex monomial groups. Now if |^>> and |^>> are families of angles resp. directions with indices in >, and which extend > and >, then we define an increasing isomorphism |^>> between |^>>,|^>>>> and |^>>,|^>>>> by |^>:\\,\>>f>*\\\\,\>>*e>>*\>*\>>f>*e>>*|\>(\).> We notice that |^>> extends > if and only if |\>=Id>, which is again equivalent to the condition that for each \{log z,log log z,\}> we have >=\>\ \ \>*\>=\>*\>.> In this case, we say that ,\)> and ,\)> are . We say that ,\)> and ,\)> are if the relation holds for all =log z> with sufficiently large \>. > Assume now that we are given a complex field of transseries =C>>, which is not stable under exponentiation (modulo the extension of the exponentiation as described in remark ). Let > and > be the associated families of angles and directions. Now consider the formal complex monomial group =exp \>,> whose asymptotic ordering is given by exp g\f\g>, for all \>>. Given extensions > and > of > and > to families indexed by monomials in >, we may totally order =C>> as explained in section . It is easily verified that > is a pre-field of complex transseries, which we call the of , relative to > and >. In cases of confusion, we will write ,\>> instead of . Notice that the exponential of any series in is defined in >. Again, the structure of ,\>> does not really depend on the choice of ,\)>. Indeed, if |^>>,|^>>)> and |~>>,|~>>)> are two different such choices, then :\\>f>*\\\\>*e>>*\>*\>>f>*e>>*\> is an increasing isomorphism between |^>>,|^>>>> and |~>>,|~>>>>. Starting with from the previous section, we may now consider the iterated exponential extensions =C>=\>, =C>=\>, =C>=\> of . The union =C=\\\\\\\=C\\\\\\>> of these fields is called a field of complex transseries in . Of course, the construction of depends on the successive choices of angles > and directions > for >, with indices in >. The angles > and directions > for coincide with these choices on each >. We will write ,\>> instead of whenever confusion may arise. We claim that ,\>> and ,\>> are isomorphic as soon as the restrictions of ,\)> and ,\)> to }> are compatible. We have already shown (see formulas () and ()) that there exist isomorphisms :\,\>\\,\>> for each . Now let > be such that z>=\ z>> for each l>. Then we observe that (log z)=\(log z)> for all \> and i\l-l>. By induction over l>, it then follows that for (\)=\(\)> for all \\>> and i>. Given \\>, this shows that the value of (\)> does not depend on the choice of , for sufficiently large . In other words, the > can be glued together into an isomorphism between ,\>> and ,\>>. It is possible to slightly generalize the construction of pre-fields of complex transseries in , when starting with =C*log z\C*log z\\> instead of =C*log z\log z\\>. Notice that is not necessarily a monomial when adopting this generalization.> Actually, in our construction of pre-fields of complex transseries in , it is reasonable to require that require that x>=0> for all sufficiently large \>, thereby eliminating all ambiguity (up to isomorphism) in the construction of . More generally, a pre-field of complex transseries is a , if it satisfies the following axiom: <\description> For each \\>, there exists an \\>, such that for all i> we have <\itemize> (log \)=log \(log \)>. (log \)>=0>. Then up to isomorphism, we have constructed field of grid-based complex transseries in . Actually, the same procedure of exponential extensions and direct limits can be used to close any field of complex transseries under exponentiation. Again, this closure is unique up to isomorphism. <\remark> In this paper we restrict our attention to grid-based complex transseries. Nevertheless, the results of this sections can easily be generalized to the case of Noetherian complex transseries. In this case, we recommend to replace the axiom by the following more complicated but better axiom: <\description> Let )\>> be a sequence of monomials in >, such that \supp log \>. Then there exists an \\>, such that for all i> we have <\itemize> \\> for all \log supp \>. >=0>. This axiom allows the resolution of certain functional equations like +f(log z)>,> which admits natural solutions of the form +e+e>>>,> which are called . > Consider a tuple =(\,\,\)> of non zero complex transseries in with \\\\\>. We call > a if the following conditions are satisfied: <\description> =log z> for some \>, which is called the of >. \C;\;\>> for each 1>. \\\\> (i.e. >\\\\>>). Such a transbasis generates a >. We say that C> can be expanded w.r.t. > if C>>. If , then we say that > (and any C>>) is . The following is proved in a similar way as in the case of real transseries: > be a transbasis and C> a complex transseries. Then can be expanded w.r.t. a super-transbasis |^>\\> of >.> We define a strong derivation w.r.t. on =C> in the usual way: we take >*\*(log z)>)=|z>+\+|z*log z*\*log z>z>*\*(log z)>> for all monomials >*\*(log z)>\\>. This yields a derivation on > through extension by strong linearity. Given a derivation on >, we define =f*exp f,> for all monomials exp (\)>=\>. This again yields a derivation on > through extension by strong linearity. By induction over , we thus obtain a derivation on . We recall that a derivation on is said to be resp. if the following conditions are satisfied: <\itemize> g\f\g>, for all \> with 1>; 1\(f\0\f\0)>, for all \>. Contrary to the case of real transseries, our derivation on cannot be strictly positive. Indeed, either \1> or \1>, say \1>. Then we have )=-e>, so either )\sign e> or )\sign (e)>. On the other hand, the following may be proved in the usual way: > is strictly valuated.> Actually, the proof involves upward shiftings of transseries: given C>, its upward (resp. downward) shifting is defined by =f\exp> (resp. =f\log>). Contrary to the case of real transseries, this transseries does not necessarily live in the same field of transseries as : if C,\>>, then we have \C\,\\>>, where \\>=\>> and \\>=\>> for all transmonomials \\>. In the case of downward shiftings, one may have to consider the generalized fields of complex transseries in from remark . It is more difficult to extend functional composition from the real to the complex setting due to possible incompatibility between the angles and directions. For instance, if =0>, then the transseries +e+\> can not be composed on the right with . In general, right composition with a given transseries is only defined on a certain subfield of >. Contrary to the case of real transseries, certain functional equations like +f(\*z),> with \C> seem to fall outside the scope of the theory of complex transseries, unless someone comes up with some really new ideas to incorporate the solutions to such equations inside this theory. One of the main ideas behind the construction of fields of complex transseries is that we do not longer require the ordering on the constant field to be compatible with the multiplication. Indeed, we just need the compatibility with the addition (or multiplication with reals), in order to obtain ordered monomial groups via exponentiation. The above idea may be used to generalize the results from this section to other circumstances. Consider for instance the set > of -adic complex numbers, where 2>. Then it is classical that there exists a partial logarithm on , which is defined for all C> with z=0>. By Zorn's lemma, there exists a total ordering on the -vector space . The theory of this section may now be adapted in order to construct the field > of complex -adic transseries. A first change concerns the condition , which should now become \>\|val c=0}={f\\>\|c\dom log}.> Furthermore, it is not as easy as before to characterize the total orderings on , which are compatible with the -vector space structure. Consequently, there is no natural analogue to the condition and we have to satisfy ourselves with the construction of pre-fields of complex transseries. Also, the exponentiation on > is not total. Notice that it seems to be possible to take itself for the indeterminate in the construction of >. This would yield a field of transseries which contains > and such that the logarithm is defined for all non zero elements. In practical computations with complex transseries the angles > and directions > are not known in advance and we have to choose them (or more precisely, to put constraints on them) as the computation progresses. This can be done by introducing a closed interval >\R/(2*\*\)> for each transmonomial >, which corresponds to the constraint \\\>:\|\-\>\|\\/2> on >>. Given such sets >>, we will work with which are in the ``intersection'' of all ,\>> such that > and > satisfy the above constraints. Actually, it is convenient to always work w.r.t. generic complex transbases, which we will introduce now. Let =(\,\,\)> be an -tuple of symbols. Assume that each > comes with closed interval \\/2*\*i> modulo >, such that \(\+\)=\>. Then we may order the monomial group =\*\*\> by >*\*\>\1\arg \\\,> for each non zero monomial >*\*\>> with \0>. We call > a of the scale >. Such a basis is called a , if <\description> =log z> for some \>, which is called the of >, and \>. > is a regular, infinitely large transseries in ;\;\>> for each 1>. >\\\\>>. An important question is whether the asymptotic constraints on the > determine a non empty region of the complex transplane (see chapter 6 of ). This question will be addressed in a forthcoming paper. ,e/(1-z)>)> is a transbasis, for the constraints z\e\e/(1-z)>>. Computations with respect to this transbasis are valid in regions of , where /(1-z)>\|\\|e\|\\|z\|\1>. This is for instance the case for , such that +\> in a region where +\\>\k+-\> for some small \(0,)> and \>.> A is an element of ;\;\>> for some complex transbasis ,\,\)>. It can be shown that two transbases which have a non empty region of definition in common can be merged together. In the remainder of the paper we will follow an easier approach, which consists of working with respect to a , which may be enlarged and on which we may impose additional asymptotic constraints during computations with complex transseries. By construction, all ring operations can already be carried out in an algebra of the form ;\;\>>. In order to invert a complex transseries, we first have to be able to compute its dominant monomial. In principle, both or > might be ``the'' dominant monomial of a transseries like >. Nevertheless, given a transseries C;\;\>> with dominant monomials ,\,\>, then we may always separate cases \\\\\\\\\\\\>>|\\\\\\\\\\\\>>|>>|\\\\\\\\\\\\>>>>>,> in each of which has only one dominant monomial. This case separation technique is explained in detail in . In the present context, the imposition of a constraint >*\*\>\1> with \0> reduces to the insertion of > in the interval >>. If the length of the new interval exceeds \>, then () can not be satisfied, so that the corresponding case does not need to be considered. <\remark> In order to be really complete, we should also consider the cases when several dominant monomials are asymptotic. For instance, in the case of the series >, we should consider the cases e> and e>, but also e>. However, in the present paper, we argue that the situation when e> is ``degenerate'' in the sense that it corresponds to a single ``direction'' 0[\]> among a continuous number if possibilities. As a consequence, we notice that the process of ``regularization'' of a complex transseries is much easier than in the case of multivariate transseries studied in . Indeed, in the case when one has to consider the possibility that e>, one also has to consider the possibility of cancellation =0> or \1>. This would necessitate refinements of the coordinates and rewriting of the series in ;\;\>>. <\example> Modulo cases separations, we may thus carry out all field operations. For instance, the inverse of > is either given by >=1-e+e+\,e\1,> or >=e-e+e+\,e\1,> Consider a non zero complex transseries C;\;\>>. Modulo case separations, we may assume that is regular, so that we can write >*\*\>*(1+\),> with ,\,\\C> and \1>. Consequently, *log \+\+\*log \+c+log(1+\).> If =0>, then this series is already in ;\;\>>. Otherwise, it still is, modulo the insertion of a new element =log \=log z> in front of the transbasis, subject to the constraint \>. Since > is a new symbol, this constraint is non contradictory with the existing expo-linear constraints on the >. The relation \\>> is automatically verified, since \>\z>. Consider a complex transseries C;\;\>>. Modulo case separations, we may assume that is regular. In order to compute the exponential of , we distinguish three cases: is bounded. We may write >, with C> and \1>. Hence, =e*e>>, with C> and >\C;\;\>>. \\\log \> for some i\n> (where we understand that the left resp. right hand side relation is verified if resp. ). We decompose +f>, where =[\*\*\] f\C;\;\>> and \1>. Inserting >> into > by \(\,\,\,e>,\,\,\)>, we then have =e>*e>\C;\;\;e>;\;\;\>>. \log \> for some . We may write *log \+g>, with \C>> and f>. Then =\>*e> and we compute > using the same algorithm. The computation of > cannot give rise to infinite loops, since the transbasis > would remain invariant in such a loop, while the index would strictly decrease. Consider the complex ``exp-log function'' e+i*z>+e>> and let us show how to expand it generically with respect to a generic complex transbasis. We start with \(z)> and recursively expand all subexpressions of . >.>In order to expand >, we fall into the second case of the exponentiation algorithm, since z> and 1>. Consequently, we insert > into > using \(z,e>), so that > expands as >. and +i*z>.>Since >> is a ring, we immediately have +i*z\C>>. Since the expansions of sums and products do not present any problems, we will omit them in what follows. +i*z>>.>In order to expand +i*z>>, we first have to determine the dominant of +i*z>. Two cases need to be distinguished for this, namely >={0}>, which corresponds to \1>, and >={\}>, which corresponds to \1>. In the first case, \i*z\log e>, so that +i*z>> needs to be inserted into >. In the second case, \i*z>, so we rewrite +i*z>=e*e>=e+e+*e+\\C>>. >>.>In the case when >={0}>, we have \log e+i*z>>, so we rewrite >=(e+i*z>)*e\C;e+i*z>>>. In the other case, when >={\}>, the argument > is bounded, so that >=1+i*e-*e+\\C>>. .>We first have to determine the dominant monomial of +i*z>+e>>. If >={0}>, then we separate the cases +i*z>>={-|4>}> in which +i*z>\e>>, and +i*z>>={|4>}> in which +i*z>\e>>. In the first case, we obtain ||+i*z+log (1+e-i*z>)>>|||+i*z+e-i*z>+>*e-2*i*z>+\\C;e+i*z>>.>>>>> In the second case, we get ||+log (1+e+i*z>)>>|||+e+i*z>+>*e+2*i*z>+\\C;e+i*z>>.>>>>> If >={\}>, then +i*z>+e>=e+1+e+i*e+\>, so we separate the cases >=[|2>,\]> in which \1>, and >=[\,|2>]> in which \1>. If >=[|2>,\]>, then ||>-1)+e*e>)>>|||+e+(i-1)*e->*e+\\C>.>>>>> Otherwise, we obtain ||>-1)+e*e>)>>|||-e+(1+i)*e->*e+\\C>.>>>>> |A plot of the function e+i*z>+e>>, which illustrates the four possible asymptotic behaviours of on non degenerate regions. The ``rows'' of singularities correspond to the borders between regions of different types.> In order to deal with integration constants when solving differential equations, we need to consider parameterized transseries. As in the case of generic transseries, if will often be necessary to distinguish several cases as a function of the values of the parameters. Again, this can be done by putting constraints on the parameters. Let =(\,\,\>)> be a >-tuple of complex parameters. We call a subset > of \>> a , if > is the set of solutions of a system of polynomial equations or inequations )>||>|)>|>|>>>> where C[\,\,\>]>, and ``rational function inequalities on the real parts'' (\)|c(\)>\0,> where ,c\C(\,\,\>)> and > does not vanish on >. Notice that > may be seen as a special kind of semi-algebraic set, under the isomorphism \>\R>>. The polynomial algebra ,\,\>]> will also be called the or . Given a non empty region \C\>>, let =(\,\,\)> be an -tuple of symbols. Assume that each > comes with a finite set =\>={\,\,\>}\P> of , such that > does not vanish on > for all i\r>, j\d>, and such that \\Cp>Re (\)|\>(\)>\0\\,> for all i\r> and j,j\d> with \j>. In the case when =0>, the directions > correspond to the extremal angles in the intervals > from the previous section. For each i\r>, there exists a natural partial ordering > on the -vector space , which is generated by the relations \0> for all . Indeed, the constraints () in an arbitrary point \\> guarantee the absence of relations *\+\+\>*\>=0,> with ,\,\>)\(R)>\\(0,\,0)>. Consequently, we may define a natural neglection relation > on the asymptotic scale =\*\*\> by >*\*\>\1\\\ 0,> for each non zero monomial >*\*\>> with \0>. We say that > is a , if <\description> =log z> for some \>, which is called the of >, and \>. > is a regular, infinitely large transseries in ;\;\>> for each 1>. >\\\\>>. A parameterized transseries is an element of ;\;\>> for some transbasis ,\,\)>. A regular parameterized transseries P>> is said to be , if either , or >(\)\0> for all \\>. In this section we prove that any parameterized transseries P>> can be uniformly regularized modulo case separations. We notice that a uniformly regular parameterized transseries on a region > remains uniformly regular on any subregion of >. Let \\*\*\> be a monomial. Then, modulo case separations, we may assume that either \1>, =1> or \1>.> <\proof> Write =\>*\*\>>, with ,\,\\P> and separate the following cases : <\itemize> For some i\n>, we have \0,\=0,\,\=0>; For some i\n>, we have \0,\=0,\,\=0>; =\=\=0>. In the cases , we have \1>. In the cases we have \1>. In case , we have =1>. Notice that the imposition of the constraints of the form =0> or \0> may involve a reduction of the region > and/or the insertion of new directions into >. Indeed, =0> is an additional algebraic constraint on >. In order to impose \0>, we first impose the constraints \0> and |\>\0> on >, for all j\d>. Next we insert > into >. Let ,\,\> be infinitesimal monomials in an arbitrary monomial group > with -powers, such that =\>*\*\>>, for certain ,\,\\\>. Then there exist infinitesimal monomials ,\,\\\>, such that \{\,\,\}>> for all i\k+1>.> <\proof> Since \1>, we may assume without loss of generality that \0>, modulo a permutation of indices. We will prove the lemma by induction over . For the lemma is trivial. So assume that 1> and let =\>*\*\>>. Then we have either \1>, =1> or \1>. If \1>, then there exist ,\,\\\>, such that ,\,\,\\{\,\,\}>>, by the induction hypothesis. Consequently, ,\,\\{\,\,\,\}>>. If =1>, then =\>>, whence ,\,\\{\,\,\}>>. If \1>, then there exist ,\,\\\>, such that >,\,\>,\>\{\,\,\}>>, by the induction hypothesis. Hence ,\,\\{\,\,\,\*\>}>>. The lemma follows by induction. Any P>> can be uniformly regularized modulo case separations.> <\proof> Let \1,\,\\1,\,\,\\\*\*\> be such that {\,\,\}>*{\,\,\}>. By lemma , we may assume without loss of generality that either \1>, =1> or \1> for each , modulo some case separations. Without loss of generality, we may therefore assume that admits a Cartesian representation in ,\,\>, i.e. f\{\,\,\}>*\>*\*\>> for certain ,\,\\\>. Choosing minimal, we will prove the theorem by induction over . If , then , and we have nothing to prove. So assume that 0>. We will first show how to regularize modulo case separations. So let ,\,\> be the set of dominant monomials of . By repeated application of lemma , and modulo reordering, we may assume that \\\\>. If all these inequalities are strict, then we are done, since > will be the only dominant monomial. Otherwise, we have =\> for certain j>, which yields a non trivial relation >*\*\>=1> for certain ,\,\\\>. Then lemma implies that we may find a Cartesian representation for in variables only, and we are done again, by the induction hypothesis. In order to make uniformly regular modulo case separations, we use the following algorithm: <\itemize> Regularize modulo case separations and let > be its dominant monomial (if 0>). If , or >(\)\0> for all \\>, then we are done. Separate the cases when >=0> and >\0> and go back to step 1. We have to show that this algorithm terminates. Assume the contrary and let ,\,\> be the successive dominant monomials of in step 1 on smaller and smaller subregions \\\\> of >. Ultimately, for each , there exists a \\> with >(\)=0> in step 2, and the next region is given by ={\\\\|f>(\)=0}> in step 3. Now the numerators of all coefficients >,f>,\> belong to the Noetherian polynomial ring ,\,\>]>. Consequently, the increasing chain of ideals >)\(f>,f>)\\> is stationary and so is the decreasing chain \\\\> of subregions of >: contradiction. Using the tool of uniform regularization, we may compute with parameterized complex transseries in a similar way as explained in sections , and . Of course, it may happen that we need to exponentiate or to take logarithms of parameterized constants in . Nevertheless, this can only happen a finite number of times, so that we may see these exponentials resp. logarithms as new parameters. Furthermore, we will show that it is never necessary to exponentiate or take logarithms of parameterized constants during the resolution of algebraic differential equations. <\example> Consider the expansion of the function +\*z>+e*e>).> \1>.>We insert into >> and get +\*z>+e*e>=e*z>+e+1)*z>+\+1+\*e+\.> We thus have to determine whether *z>\1> and *z>\1>, which leads to the following cases and expansions for : |+|2>>*e+-+|4>-|8>>**e+\>|>|=0,\>={});>>|*z>+\+\*e+\>||\0,\>={,>});>>|*z+e+\+e*z>+\>||\0,\>={,\}).>>>>>> \1>.>We insert into >> and next need to determine whether +\*z>\e*e+\*\*z>> or +\*z>\e*e+\*\*z>>. This leads to the following cases and expansions for : |+\*z+e*\*z>*e-1)*(e+\*z)>+\>|>|\\1,\>={},\+\*z>>={1-\};>>|*e+e*\*z>*e)*(e+\*z)>+\>||\0,\>={},\+\*z>>={\-1}).>>>>>> In the last exceptional case when =1>, we get +\*z>*(1+e*z>)),> so that we need to determine whether e*z>> or e*z>>. This leads to the following final cases and expansions for : |+\*z+log 2>|>|=0,\=1,\>={1});>>|+e*z>->*e*z>+\>||\0,\=1,\>={1,-\});>>|+\*z+e*z>+\>||\0,\=1,\>={1,\}).>>>>>> In the remainder of this paper, we will be concerned with the resolution of asymptotic algebraic differential equations like (f\\),> where \[f,f,\,f]> is a differential polynomial with transseries coefficients and \\> a transmonomial. In this section, we describe the differential Newton polygon method, which enables us to compute the successive terms of solutions one by one. In the next sections, we will be concerned with the transformation of this transfinite process into a finite algorithm. In sections , and the transseries in are assumed to be as in section . In section , we will consider parameterized transseries solutions. Except for the usual asymptotic relations ,\,\,\,\,\,\> and >, we will also need the flattened relations ,\,\> and there variants >,\>,\>>, where is an infinitely large or small transseries. These relations are defined by g>|>|\\h:f*\\g;>>|g>|>|\\h:f*\\g;>>|g>|>|g\g\f;>>|>g>|>|\\log h:f*\\g;>>|>g>|>|\\log h:f*\\g;>>|>g>|>|>g\g\>f.>>>>> Notice that >g\f\g>, >g\f\g> and >g\f\g>. > The differential polynomial is most naturally decomposed as P(f)= > P>*f>> Here we use vector notation for tuples =(i,\,i)> and =(j,\,j)> of integers: <\eqnarray*> \|>||>|\\<\|\|\>>||+\+i;>>|\\>|>|\j\\\i\j;>>|>>||i>*(f)>*\*(f)>;>>||\>>|||i>*\*|i>.>>>> The -th of is defined by =\\<\|\|\>=i>P>*f>,> so that P.> along orders> Another very useful decomposition of is its : P(f)= > P]> f]>> In this notation, > runs through tuples =(\,\,\)> of integers in ,r}> of length deg P>, and ]>=P(1)>,\,\(l)>]>> for all permutations of integers. We again use vector notation for such tuples <\eqnarray*> \|>||>|\\<\|\|\>>||+*\*+\\|>;>>|\\>|>|\|=\|\\|\\\\\\\\\|>\\\|>;>>|]>>||)>*\*f\|>)>;>>||\>>|||\>*\*\|>|\\|>>.>>>> We call \>\|\| the of > and P\<\|\|\>=max\|P]>\0>\<\|\|\>\\<\|\|\>> the of . > It is convenient to denote the successive logarithmic derivatives of by >>||/f;>>|i\>>||\\\> times)>.>>>>> Then each > can be rewritten as a polynomial in >,\,fi\>>: ||>|>||>*f;>>|>||>)+f\>*f>)*f;>>|>||>)+3*f\>*(f>)+(f\>)*f>+f\\>*f\>*f>)*f;>>||>|>>>> We define the of by =(i,\,i)>P\\>*f\\>,> where \\>=f>*(f>)>*\*(fr\>)>.> Now consider the lexicographical ordering > on >, defined by \\>|>|\j)\>>|||=j\i\j)\>>|||\>>|||=j\\\i=j\i\j).>>>>> This ordering is total, so there exists a maximal for > with \\>\0>, assuming that 0>. For this , we have )\P\\>*f\\>> for all , whose dominant monomial is sufficiently large. Given a differential polynomial and a transseries it is useful to define the and > and h>> of w.r.t. and the > of as being the unique differential polynomials, such that for all , we have (f)>||>|h>(f)>||>|(f\)>||.>>>>> The coefficients of > are explicitly given by P>= \\> |\>*h-\>*P>.> The coefficients of h>> are more easily expressed using decompositions along orders: Ph,[\]>= \\>|\>*h- \]>*P]>.> The coefficients of the upward shifting (or compositional conjugation by >) are given by (P\)]>=\\>s,\>*e\\<\|\|\>*z>*(P]>\),> where the ,\>> are generalized Stirling numbers of the first kind: |||||,\>>|=>|,\>*\*s\|>,\\|>>;>>|>|| s*x*f(log z).>>>>>> Given a differential polynomial with transseries coefficients, its > is defined by =max,\> \>>.> and its (or coefficient) \C[c,c,\,c]> by =>P,\>*c>.> The following theorem shows how > looks like after sufficiently many upward shiftings: Let be a differential polynomial with purely exponential coefficients. Then there exists a polynomial C[c]> and an integer >, such that for all \<\|\|\>P\<\|\|\>>, we have >=Q*(c)>>.> <\proof> Let > be minimal, such that there exists an > with \\<\|\|\>=\> and \)]>\0>. Then we have (D\)=e*z>> and >(c)=\\<\|\|\>=\>\\>s,\>*D]>*c]>,> by formula (). Since >\0>, we must have \\<\|\|\>D\<\|\|\>>. Consequently, D\<\|\|\>\\=\<\|\|\>D>\<\|\|\>\\<\|\|\>D\>\<\|\|\>\\>. Hence, for some >P\<\|\|\>>, we have D>\<\|\|\>=\<\|\|\>D>\<\|\|\>>. But then () applied on > instead of yields >=D>>. This shows that >> is independent of , for \<\|\|\>P\<\|\|\>>. In order to prove the proposition, it now suffices to show that >=D> implies >=Q*(c)>> for some polynomial C[c]>. For all differential polynomials of homogeneous weight >, let >=([c*(c)>]R)*c*(c)>.> Since >>=D>>, it suffices to show that if and only if >=0>. Now >=0> implies that (z)=0>. Furthermore, () yields \=e*z>*D.> Consequently, we also have (e)=e*z>*(D\)(e)=e*z>*(D(z))\=0>. By induction, it follows that (exp z)=0> for any iterated exponential of . We conclude that =P=0>, by (). Given an arbitrary differential polynomial , the above proposition implies that there exists a polynomial C[c]> and an integer >, such that >=Q*(c)>> for all sufficiently large . We call =Q*(c)>> the of . More generally, given a monomial >, we call \>>> the of >. Returning to the asymptotic differential equation (), we call \\> a , if \>>> admits a non trivial root (C)>>, where > stands for the algebraic closure of . If C\>>, then the corresponding term > is called a . The of (and of >) is the differential valuation of >,+c>>, i.e. the least such that >,+c,i>\0>. The of () is the largest possible degree of \>>> for monomials \\>. >| is a regular, non-zero transseries solution to )>. Then > is a potential dominant term.>> A potential dominant monomial > is said to be if \>>> is non homogeneous, and if \>>\C[c]>. A potential dominant monomial, which is both algebraic and differential, is said to be . Notice that () implies (P\>)\>>\*\(P),> if the coefficients of and > are purely exponential. The algebraic potential dominant monomials correspond to the slopes of the Newton polygon in a non differential setting. However, they can not be determined directly as a function of the dominant monomials of the >>, because there may be some cancellation of terms in the different homogeneous parts during multiplicative conjugations. Instead, the algebraic potential dominant monomials are determined by successive approximation: <\proposition> Let j> be such that \0> and \0>. <\enumerate-alpha> If is purely exponential, then there exists a unique purely exponential monomial >, such that (P\>)=\(P\>)>. Denoting by > the monomial > in )>, there exists an integer \<\|\|\>P\<\|\|\>>, such that for all k> we have ,i,j>=\,i,j>\>. There exists a unique monomial >, such that +P)\>>> is non homogeneous. <\proof> In (), let =(\,\,\)> be a purely exponential transbasis for the coefficients of . We prove the existence of > by induction over the least possible , such that we may write (P)/\(P)=\>*\*\>>. If , then we have =1>. Otherwise, let \>> with =\/(j-i)>>. Then (Q)\>\(P)*\\>\(P)*\\>\(Q),> so that (Q)/\(Q)=\>*\*\>> for some k> and ,\,\>. By the induction hypothesis, there exists a purely exponential monomial >, such that (Q\>)=\(Q>)>. Hence we may take =\*\>. As to the uniqueness of >, assume that =\*\>*\*\>> with \0>. Then (P\>)\>\(P\>)*\>\>\(P>)*\>\>\(P\>).> This proves (). With the notations from proposition , we have already shown that D\>\<\|\|\>\\<\|\|\>D>\<\|\|\>> and that equality occurs if and only if >=cD>\<\|\|\>>*(c)D>\<\|\|\>>>. Because of (), we also notice that D,\*e*z>>\<\|\|\>=\<\|\|\>D>\<\|\|\>> for all >. It follows that D,\\>\<\|\|\>\\<\|\|\>D\,\\,i,j>>\<\|\|\>\\> and similarly for > instead of >, since we necessarily have \,i,j>=\,i,j>\*e*z>> for some >. We finally notice that D,\\>\<\|\|\>=\<\|\|\>D\,\\,i,j>>\<\|\|\>> and D,\\>\<\|\|\>=\<\|\|\>D\,\\,i,j>>\<\|\|\>> imply that ,i,j>=\\>, since D>*(c)>)ex>>>\<\|\|\>=0\\=\<\|\|\>D>*(c)>>\<\|\|\>> whenever \0> and \0>. Consequently, D\,\\,i,j>>\<\|\|\>> and D\,\\,i,j>>\<\|\|\>> stabilize for k> with \<\|\|\>P\<\|\|\>>. For this , we have (). With the notations from (), ,i,j>\> is actually the unique monomial > such that +P)\>\>=D\>\>+D\>\>> is non homogeneous for all sufficiently large . Now +P)\>>=D+P)\>\>> for sufficiently large . This proves () for purely exponential differential polynomials , and also for general differential polynomials, after sufficiently many upward shiftings. The unique monomial > from part () of the above proposition is called an or the -equalizer for . An algebraic potential dominant monomial is necessarily an equalizer. Consequently, there are only a finite number of algebraic potential dominant monomials and they can be found as described in the proof of proposition . Notice that, given a transbasis =(\,\,\)> for the coefficients of , all equalizers for belong to P\<\|\|\>> \)*\*(log \)*\>. In order to find the differential potential dominant monomials, it suffices to consider the homogeneous parts > of , since \>,i>=N\>>>, if \|N\>>> and \>,i>\0>. Now we may rewrite > as i>> times a differential polynomial > of order 1> in >>. We call > the -th Ricatti equation associated to . Since solving (f)=0> is equivalent to solving (f>)=0>, we are entitled to expect that finding the potential dominant monomials of w.r.t.> is equivalent to solving (f>)=0> ``up to a certain extent''. <\proposition> The monomial \\> is a potential dominant monomial of w.r.t. (f)=0> if and only if the equation >>(f>)=0f>\>>> has strictly positive Newton degree. <\proof> We first notice that ,i>=(R\)e>> for all and . We claim that the equivalence of the proposition holds for and > if and only if it holds for > and \>. Indeed, > is potential dominant monomial w.r.t. (), if and only if > is a potential dominant monomial w.r.t. \(f\)=0> and () has strictly positive Newton degree if and only if >>\(f>\)=0f>\\*z*log z**\>>> has strictly positive Newton degree. Now the latter is the case if and only if >>\)e>(f\>)=0f\>\>>> has strictly positive Newton degree. But >>\)e>=(R\)>\,\e>=(R\)e,+\\>>=R,i,+\\>>.> This proves our claim. Now assume that > is a potential dominant monomial w.r.t.)>. In view of our claim, we may assume without loss of generality that and > are purely exponential and that \>>=D\>>>. Since > is homogeneous, we have \>>=\*(c)> for some \C\>> and >>>=\*ci>.> Since >>> is purely exponential, it follows that >,\z>>> has degree , so that the Newton degree of () is at least . Similarly, if > is not a potential dominant monomial w.r.t. (), then \>>=\*c> and >>>=\*> for some \C\>>. Consequently, >,\\>>=\> for any infinitesimal monomial >, and the Newton degree of () vanishes. Now we know how to determine potential dominant terms of solutions to (), let us show how to obtain more terms. A is a change of variables together with an asymptotic constraint +(\),> where \\>. Such a refinement transforms () into >()=0(\).> We call the refinement , if () has strictly positive Newton degree. Let > be the dominant term of > and assume that =\>. Then the Newton degree of )> is equal to the multiplicity > of as a root of \>>>.> <\proof> Let us first show that ,\\>>\> \ for any monomial \\>. Modulo replacing by \>> we may assume without loss of generality that =1>. Modulo a sufficient number of upward shiftings, we may also assume that =D>, that ,\\>>=D,\\>>>, and that , > and > are purely exponential. The differential valuation of =D>> being >, we have in particular (P,>)=\(P>)>. Hence, (P,\\,i>)\>\(P,\\,i>)*\\>\(P,\\>)*\>\>\(P,\\,>)> for all >. We infer that ,\\>>\>. At a second stage, we have to show that ,\\>>\>. Without loss of generality, we may again assume that =1>, that =D>, and that and > are purely exponential. The differential valuation of =D>> being >, we have (P,i>)\\(P>)> for all >. Taking =z>, we thus get (P,\\,i>)\>\(P,i>)\>\(P>)=\(P,>)\>\(P,\\,>)> for all >. We conclude that ,\\>>\>. Consider the algebraic differential equation P(f)=f+f*f-(f)=0.> Let us start by computing the potential dominant monomials of . We first have to find the -equalizer relative to (). Since >\c>*(c)>>, we cannot have >=P>, so we have to compute =f+e*(-f*f+f*f-(f)).> In order to ``equalize'' > and >, we have to conjugate multiplicatively with >: e>=e*(f-2*f-f*f+f*f-(f)).> At this point, we observe that e>\>=c-2*c\C[c]>, so we have found the -equalizer, which is =e\=z>. Since \>>=c-2*c>, the corresponding algebraic potential dominant term of is =*z>. As to the differential potential dominant monomials, we have >||>|>||>.>>>>> Clearly, > has no roots and (f>)=0> has all constants \C> as its solutions modulo )>. Consequently, *z>> is a potential dominant monomial of for all \C>, such that *z>\1>. The corresponding differential potential dominant terms are of the form ,\>=\*e*z>>, with \0> and *z>\1>. In order to find more terms of the solution to (), we have to refine the equation. First of all, consider the refinement +(\\),> which transforms () into 2*-2*z*+>*z*+*-()=0(\z).> Since =0>, we first observe that *z> is actually a solution to (). On the other hand, since > is a potential dominant term of multiplicity 1 of , the Newton degree of () is one. The only potential dominant monomials of > therefore necessarily correspond to solutions modulo )> of the Ricatti equation >+>*z*((>)+f'>)=0.> These solutions are of the form >=+\> and >=+\>, which leads to the potential dominant monomials and >, from which we remove >, since \z>. Expanding one term further, we see that the generic solution to () is =\*z+|2>>,> with \C> and where the case =0> recovers the previous solution. In other words, >*z+\*z+|2>>> is the first type of generic solution to (). As to the second case, we consider the refinement ,\>+(\\,\>),> which transforms () into \*e*z>+(\*f-2*\*f+f)*\*e*z>+f+*-()=0(\e*z>).> Again, this equation has Newton degree one. On the one hand, we observe that the linear part of this equation only admits solutions with dominant monomial *z>> or *z>>. Consequently, () admits at most one solution. On the other hand, we will show in the next section that quasi-linear equations ( of Newton degree one) always admit at least one solution. In our case, this leads to the following second type of generic solution to (): *e*z>->>+*\>>*e*z>.> For the present example, we actually even found exact solutions. Of course, the expansions are infinite in general. Let =(\,\,\)> be a purely exponential transbasis. A linear operator on =C[z];\;\>\C;\;\>> is said to be if its > L=>|\>> is grid-based. For all transseries \> we have (supp> L)*(supp f).> In particular, the differentiation > on is grid-based with > \\{z}\supp \>\\\supp \>.> Consequently, any linear differential operator +L*\+\+L*\r>> with coefficients in ;\;\>> is also grid-based, since > L\supp L\(supp L)*(supp> \)\\\(supp L)*(supp> \).> We will now show that also admits a so called L>, which is linear and grid-based. Here a to the equation is a solution , such that for all other solutions >, we have -f>>=0>. Distinguished solutions are clearly unique. We say that > is a distinguished left inverse of , if g> is a distinguished solution to for each . In what follows, we will often consider linear differential operators as linear differential polynomials. In this case, you should keep in mind that > denotes the coefficient of > in and not the -th homogeneous part. We will also denote \\\supp L> for any linear differential operator as above. <\theorem> Let +L*\+\+L*\r>> be a linear differential operator with coefficients in ;\;\>> and \0>. Then admits a distinguished linear left inverse > on ;\;\>>. This left inverse is grid-based and > L\\*\>,> where |>|=>|>*|\(L\>)>\\\;>>|>|=>|>\z>*\\>\>|\(L\>)>{1}.>>>>>> <\proof> Let =z>*\> and ={\\|h\\,L h=0}>. There exists a unique strongly linear operator :C[[\\\\]]\C[[\]]>, such that \=\>> for all \\\\\>. The operator > admits a natural left inverse :C[[\]]\C[[\\\\]]>, which is constructed as follows. Let *\\\>, where > is purely exponential. By proposition (), there exists a purely exponential monomial > with (L\>)=\>. Let > and respectively denote the dominant monomial and dominant part of \>>. Let =*i!|D\>,j>*(i+j)!>*z,> where is minimal with \0>. Then we observe that \> \\(D\>> \)*\\>>\z*\,> so that \=z*\>. Consequently, we may take (z*\)=\> and extend > to the whole of ]]> by strong linearity. Let >. By construction, the operator > is strictly extensive, and the operator )*\> coincides with on \\\]]>. Now consider the functional (f,g)=g-R*\ f.> By the implicit function theorem from , there exists a linear operator =(Id+R*\)=Id-R*\+(R*\)+\,> such that (\(g),g)=\(g)> for all C[[\]]>. Consequently, =\*(Id+R*\):C[[\]]\C[[\\\\]]> is a strongly linear left inverse for . In order to prove that > is actually a grid-based operator, we first notice that, by construction, > \>|>|;>>|> (R*\)>|>|.>>>>> Since > \*(Id+R*\)\(supp> \)*(supp> (R*\))>,> it therefore suffices to prove that > and > are grid-based. But this follows from theorem when considering *L\>> as a generic transseries in ,\,\>, for =\>*\*\>>. Indeed, there exist a finite number of regions, each on which *L\>> is uniformly regular. Consequently, /\(L\>)> can only take a finite number of values and \\>\>|\\>>>> is contained in the union of the supports of the generic transseries *L\>> on each of the regions. <\remark> Each \\> induces a +L L \\C>> to . This canonical solution satisfies >=1> and >=0> for all \\\\{\}>. Actually, the canonical solutions are polynomials of degree r> in with coefficients in ;\;\>>. In order to see this, let >||**\\z>*\\|i\Card {\\\\|\\\\>\}\z*\\\};>>|>||\L)(\).>>>>> Then we observe that maps ]]> into ]]> and that > maps ]]> into ]]>. Let > be a subset of a monomial group. The notion of operator support can be extended to strongly -linear operators ]]\C[[\]]> by > M=,\\)\\>,\,\)|\*\*\>.> More generally, if :C[[\]]\C[[\]]> is a Noetherian operator, then we define its by > \=\>supp> \,> where > stands for the -th homogeneous part of >. We have (f)\(supp> \) (supp f)>> for all C[[\]]>. We say that > is , if > \> is grid-based. Let =(\,\,\)> be a purely exponential transbasis and a differential polynomial with coefficients in ;\;\>>. Notice that we may naturally consider as a grid-based operator on ;\;\>>. The equation () is said to be if its Newton degree is one. A solution to such an equation is again said to be if we have (-f)>=0> for all other solutions > to (). Assume that the equation )> is quasi-linear. Then it admits a distinguished transseries solution.> <\proof> Without loss of generality, we may assume that =1> and =1>. We prove the proposition by induction over . If , then we must have =0>, so that is the distinguished solution to (). So assume that 0> and let ,\\>1>P,\>*\*f\>> be the dominant part of w.r.t. >. By the induction hypothesis, there exists a distinguished solution to the quasi-linear equation )=0(\\1).> We first proceed with the refinement +(\\),> so that ,0>=0>, and a sufficient number of upward shiftings, so that >> is purely exponential. We next decompose >> as >=\+R-,> where >||,1>;>>|>||,0>;>>|||>-P,0>-D,1>.>>>>> Let ={\\z>*\\|\\>1}>. Since ;\;\]]*\>> is stable under > and > for each \C>, the operator > is strictly extensive on ]]>. Consequently, the implicit function theorem from implies that the operator > can be inverted, like in the proof of theorem : )=Id-R*\+(R*\)+\.> In particular, =\ (Id+R*\) > is a solution to >()=0>. Furthermore, we have > (Id+R*\)\(supp> R*\)>,> so that \(supp> R*\)>*(supp g)>.> We claim that +> is the distinguished solution. Indeed, let \f> be another solution and let =\-f>>. If \>1>, then |^>=\>1>>*\> is a solution to (), so that >=\>=0>. If \>1>, then let =\>\>(-f)>*\.> Since )-P(f)=0>, we have \=0>, so that =\>> is the dominant monomial of a solution to the homogeneous equation h=0>. Consequently, >=0>, since \Im \>. <\remark> By induction over , it also follows that we need at most upward shiftings in order to express the distinguished solutions. In other words, if has coefficients in ;\;\>>, where > is purely exponential, then C z;\;z;\;\;\>>. Actually, if is the order of , then the number of upward shiftings we need is also bounded by . Indeed, denoting ={\\\\|\\>1}> and using a similar argument as in remark , we first observe that > is bijective on >> if *h=0> admits no solutions in >>. Moreover, if h=0> admits such a solution, then > has a root with the same dominant part w.r.t. >. The same observations recursively hold for all > involved in the resolution of (). Now if is the distinguished solution of (), then the linear equation h=0> admits at most solutions. Hence, there are at most transbasis elements > for which we need to make an upward shifting. Consider the linear differential equation L f=f-e*z>*f+e+\)*z>*f=1,> under the assumptions 1> and *z>\e*z>\1>. Then has coefficients in *z>>> and e*z>>=e*z>*(f+2*\*f+\*f-e*z>*f-\*e*z>*f+e+\)*z>*f)> for each \\>, so that (Le*z>>)>||+\+\)*z>;>>|e*z>> >|>|+\+\)*z>*{1,e*z>,e+\)*z>}.>>>>> Hence, theorem implies that () has a distinguished solution in *z>>> with f\z>*e+\)*z>*{e*z>,e+\)*z>}>>. Actually, it is easily seen that f\e+\)*z>*{e*z>,e+\)*z>}>> and the first terms of are given by =e+\)*z>+(\+\)*e+2*\)*z>+(\+\)*(\+2*\)*e+3*\)*z>-(\+\)*e+2*\)*z>+\.> In order to find all solutions to (), we have to solve the Ricatti equation associated to the linear part of (): g+g-e*z>*g+e+\)*z>=0.> This equation has two potential dominant terms *z>> and *z>> of multiplicities one. Consequently, we get quasi-linear equations when refining *z>+(\e*z>)> or *z>+(\e*z>)>. When setting =e*z>*h> resp. =e-\)*z>*h>, these equations are conveniently rewritten as e-\)*z>*h+h-e*z>*h-\*e*z>*h+1+\*e*z>=0(h\e-\)*z>)> and e-\)*z>*h-h+2*e-\)*z>*h+e*z>*h+(2*\-\)*e*z>*h+1+\*e*z>=0(h\e-\)*z>).> By theorem , these equations admit distinguished solutions >||-\)*z>-\*e*z>+\;>>|>||-\)*z>+\*e*z>+\.>>>>> More precisely, in the proof of theorem , and for (), we would have =h>, and -\)*z>*h-e*z>*h-\*e*z>*h+\*e*z>>. It follows that > R*\\{e*z>,e-\)*z>}>>, so that \\[[e*z>,e-\)*z>]]>. Similarly, \\[[e*z>,e-\)*z>]]>. Returning to (), we obtain the following solutions: >||*z>-e*z>-e-\)*z>-\*e*z>+\;>>|>||*z>+e-\)*z>+2*e-2*\)*z>+\*e-\)*z>+\,>>>>> which yield a basis >||>*e*z>->*e*z>--\>*e-\)*z>+|2*\>**e*z>+\>;>>|>||>*e*z>+-\>*e-\)*z>+-2*\>*e-2*\)*z>+|\-\>*e-\)*z>+\>.>>>>> of solutions to =0>. It is interesting to study the solutions +\*\+\*\> to () from an analytical point of view. Indeed, the asymptotic conditions \1> or \1> and \1> or \1> divide complex space into four non degenerate regions. However, each of these regions has infinitely many ``bounded connected components''. When moving from one connected component to another one, a ``generalized Stokes phenomenon'' occurs. Consequently, a specific formal solution to () only makes sense on a bounded connected component from the analytical point of view. Nevertheless, it is possible to give an asymptotic meaning to the generic formal solution to () on each region, by associating a ``generalized Stokes matrix'' to each connected component of the region. This issue will be detailed in a forthcoming paper. An interesting remaining question is the asymptotic behaviour of the Stokes matrices. Actually, the generalized Stokes phenomenon might be qualified as , since the Stokes phenomena occur with respect to several generalized sectors of different types. Equation () is one of the simplest examples which exhibits this multi-Stokes phenomenon. Theorem together with propositions , and suggest that the solutions to an arbitrary asymptotic algebraic equation () can be expressed using the field operations, exponentiation, logarithm and distinguished solutions of quasi-linear equations. This is indeed so, if the Newton degree decreases at each refinement in proposition . The remaining case, when the Newton degree repeatedly does not decrease in proposition , occurs when there are ``almost multiple solutions''. In order to ``unravel'' these solutions, we have to find their greatest common part. More precisely, consider an asymptotic algebraic differential equation () of Newton degree . Then an (or ) is a refinement +(\),> such that <\description> The Newton degree of >()=0(\)> equals . For any |~>\>, the Newton degree of +|~>>(|~>)=0(|~>\\(|~>))> is d>. Clearly, the series >, which is also called an , may be replaced by any other series of the form +\> with \>. From a theoretical point of view it is possible to prove a certain number of facts about unravellings. First of all, any unraveller admits a truncation >>, which is a , in the sense that >+(\>)> is an unravelling for all >> with >\>\>, and that a similar property does not hold for any proper truncation of >>. Secondly, it is possible to construct the so called > by transfinite induction: having constructed the first > terms of >, say of sum >, one looks at the equation >()=0(\\)> of Newton degree . If this equation has an algebraic potential dominant term > of multiplicity , then this term is unique, and we take it to be the next term of >. However, in what follows, we are interested in more constructive ways to obtain unravellings. For this purpose, we recall that in the more classical context of algebraic equations, multiple roots are usually found by solving the derivative (or a higher derivative) of the equation with respect to the indeterminate. In the next sections, we will describe a similar strategy in order to find the almost multiple solutions to asymptotic algebraic differential equations. The price to be paid is that we will need a sequence of so called partial unravellings (and adjusted partial unravellings) in order to construct a total unravelling. Consider an asymptotic algebraic differential equation () of Newton degree . Given a monomial \\> such that \>>> admits a root of multiplicity , we define by <\itemize> If \>>=\*(c-\)d-k>*(c)> with d>, then d-1> P\>/\ fd-1-k>*\ (f))\>>; If \>>=\*(c)>, then d-1> P\>/\ (f)d-1>)\>>. Now let \\\\>, =P>> and =Q>> be such that <\description> )=0>. The Newton degree of ()=0(\)> is . For any |~>\> with (|~>)=0>, the Newton degree of |~>>(|~>)=0(|~>\|~>)> is d>. Then the refinement +(\),> is said to be a with > as its . Notice that the equations )=0(\\\)> and (|~>)=0(|~>\)> are quasi-linear. Partial unravellings are constructed as follows. Let > and be as above. Then there exists a \\> which satisfies the conditions , and .> <\proof> We construct sequences ,\,\> and \\\> of approximations of > and >, such that all > and > satisfy the conditions and . We let > be the distinguished solution to the equation )=0(\\\)> and =\(\)>. As long as > and > do not satisfy the condition , there exists a \\> with >(\)=0>, such that the Newton degree of +\>()(\\)> is . Hence we may take =\+\> and =\(\)>. We claim that the sequences \ ,\,\> and ,,\> \ are of length at most , so that we may take their last elements for > and >. Indeed, for each , the series -\> with i> are solutions to the quasi-linear equation >(\)=0(\\\)>. Consequently, the dominant monomials of these series, which are pairwise distinct, are all dominant monomials of solutions to the homogeneous linear differential equation ,1>(h)=0>. But there are at most linearly independent solutions to this equation. <\proposition> Consider a partial unravelling )> as above, followed by a refinement =|~>+|~>(|~>\),> such that the Newton degree of |~>(|~>)=|~>>(|~>)=0(|~>\)> is equal to . Then, for =\(|~>)>, we have |>\log |>.> <\proof> Without loss of generality, we may assume that >, >, |~>>, > and > are purely exponential, that =1> and that ()=\()=1>. From it follows that > is neither a potential dominant monomial for (\)=0(\\)>, nor for (\)=0(\\)>. Proposition , applied to ,1>> and the ``non potential dominant monomial'' , therefore yields (R>,1,0>)=\(R>,1>).> Consequently, (|~>)\R>,1>|~>|>>\\(R>,1>)\\(>).> On the other hand, we have (>)|\()*>\>\log ,> so that (|~>)|>\log .> Now recall that (|~>)> is the coefficient of d-1-k>*(f)> in |~>(f)> for some . It follows that |~>>>\>>,> Now assume that > is a monomial with \>>> (so that \>>>). Then we have (|~>\,d-1>>)\>>\(|~>>>)*\\>>\\>>\(|~>\,d>>).> We conclude that the degree of |~>\>>> can not exceed . If > is chosen such that () has Newton degree , it thus follows that \>>,> which completes the proof. Proposition shows that by taking sequences of partial unravellings, we rapidly approach a total unraveling. The only problem which still remains to be solved is the appearance of highly iterated logarithms. We will first solve this problem in the particular case when the Newton degree of coincides with its normal degree. In the next subsection, we will show that the general case can be reduced to this case. In the sequel, we assume that () is an asymptotic differential equation of degree and Newton degree , such that the following additional conditions are satisfied for a certain purely exponential transbasis =(\,\,\)>: <\description> > has coefficients in ;\;\>>. ,\,P> have coefficients in ;\;\>>. (f)=0(f\\)> admits only potential dominant monomials in >. The two first conditions can clearly be met after a sufficient number of upward shiftings. In section , we will show that this is also the case for the last condition. <\proposition> Let > be a potential dominant term of multiplicity for (). Then <\enumerate-alpha> Modulo the insertion of new elements into >, we have =c*z\>*\\C*z>*\>. There exists a unique \z*C[z]*\>, such that either <\enumerate-roman> +(\)> is a total unravelling and > is >-minimal in >. The Newton degree of +(\z*\)> is . <\proof> Let us first prove (). If > is differential, then implies that > is purely exponential, so \C*\> after a suitable extension of >. If > is algebraic, then (\)> is the d)>-equalizer for each d>, since > has multiplicity . Proposition () implies that there exists a unique purely exponential monomial =e*z>*(\\)\e*(\\)> with \\>\>, such that (P\\>)=\(P\\>)> for all d>. More precisely, in the algorithm in proposition (), > is chosen such that (P\\\>)*\\>\(P\\\>)>, whence (P\\\>)=\(P\\\>)*e*z>> for some \\>, and > satisfies =\>. In particular, for , this yields \\>. We claim that \\\>. If =0>, then we have (P\\>)=\(P\\>)> and (P\\>)=\(P\\>)> for all d>. Since \\>=P\\>\>, this can only happen if \>>=D\>>>. Hence \> is the -equalizer w.r.t. for all d> and \\\>. If \0>, then implies that > is not a potential dominant monomial for \>. Consequently, the coefficients of > and > in \>>> both do not vanish. It follows that \> is the -equalizer w.r.t. and again \\\>. We prove the existence of > in () by induction over >; the uniqueness of > follows from . If =0>, then =0> clearly satisfies assumption . If \0>, then we refine +(\\),> and remark that >> satisfies the hypothesis , and , due to part (). Now consider the equation >()=0(\\)> of Newton degree . If this equation admits a potential dominant term |~>> of multiplicity with \|~>\\>, then the induction hypothesis implies that there exists a |~>\z*C[z]*\>, which satisfies the assumption or , and we may take =\+|~>>. If there does not exist such a potential dominant term |~>>, then there either do not exist potential dominant terms of multiplicity at all for (), so that holds for =\>, or such potential dominant terms do exist, and we have for =\>. Given a potential dominant term =c*z\>*\> of multiplicity , let > be as in proposition (). In case , we say by convention that +(\)> is an . In case , let =|~>+|~>(|~>\)> be a partial unraveling relative to the equation >()=0(\),> and with > as its associated monomial. Then we say that +|~>+|~>(|~>\)> is an . Notice that a partial unravelling like () always exists, by propositions and (). Notice also that we necessarily have |~>\C[z];\;\>>. Indeed, consider the differential polynomial with |~>)=0> in . Since , this differential polynomial is actually linear. Furthermore, since \\>, the coefficients of > are in ;\;\>> and the coefficients of > in ;\;\>>. We conclude that all solutions to )=0>, and in particular =|~>>, are in ;\;\>>. A consequence of our observation is that +|~>>> again satisfies the hypotheses , and , so that we may consider sequences of adjusted partial unravellings. <\proposition> Any sequence of adjusted partial unravellings >||+f(f\\);>>|>||+f(f\\);>>||>|>>>> is finite, say of length , and its composition +\+\+f(f\\)> is a total unravelling. <\proof> Let \\{\}> denote the length of the sequence of adjusted partial unravellings. For each i\l>, let =z>*\=\(\)\z>*\>. For each i\l-1>, let > denote the exponentiality of /\>. Given i\l-1>, proposition implies that |\>\log |\>,> Since \\> and \\>, this yields |\>\log |\>,> By induction, it follows that \\\\\0>. We conclude that \+1>. The composition of the sequence of adjusted partial unravellings is clearly a total unraveling. Let us now return to the case of a general asymptotic differential equation () of Newton degree . Assume that () is a partial unravelling with =\>> and that |~>> is a potential dominant monomial of multiplicity for ()=P>()=0(\).> Modulo a sufficient number of upward shiftings, we may assume that ,,|~>> and the coefficients of can be expanded w.r.t. a purely exponential transbasis =(\,\,\)>. Let > be the transbasis element such that /|~>\\>, and consider the dominant part > of |~>>> with respect to >: =\>\(|~>>)>|~>,\>*\> On the one hand, since \(\)>>=deg N\(\)>>=d>, we have (,\|~>,i>)\>\(\(\),i>)*|~>|\(\)>>\>\(\(\),d>)*|~>|\(\)>>\>\(|~>,d>),> for all d>, so that =d>. Consequently, > satisfies the conditions and from the previous section; we will see in section that it also satisfies , modulo some additional upward shiftings. On the other hand, the following proposition reduces the problem of determining the unravellings for () to a similar problem for >. In view of the previous section, this completes the effective construction of unravellings. |<\proposition> With the above notations, a refinement =|~>+|~>(|~>\)> with |~>\|~>> is a total unravelling w.r.t.)>> if and only if |~>||~>>+\||~>>> is a total unravelling w.r.t. the equation (g)=0g\||~>>.> > <\proof> Modulo a multiplicative conjugation, we may assume without loss of generality that |~>=1>. Now if () is an unravelling, then proposition implies that >\log \\log \,> so that \>1>. Actually, in the proof of proposition we showed that |~>,d-1>\>>1,> so that |~>,d-1>\0>. We infer that |~>,\\>>\d-1> for all \1> with \\>. In other words, for () to be an unravelling, it is again necessary that \>1>. The above argument shows that it suffices to prove the equivalence under the assumption that \>1>. Now we notice that for each transseries \>1> and each monomial \>1>, the dominant parts of ,\\>> and ,\\>> w.r.t.>> coincide. Consequently, ,\\>>=N,\\>>> for such > and >. In particular, we have |~>,\\>>=N|~>,\\>>> for all \> sufficiently close to >, so that the Newton degrees of |~>>(|~>)=0(|~>\)> and |~>>()=0(\)> coincide. Hence holds for () if and only if it holds for (). Similarly, for all \>, such that \>1>, the Newton degrees of >(|~>)=0(|~>\\>)> and >()=0(\\>)> coincide. Furthermore, for a similar reason as above, the Newton degrees of () and () are both bounded by if \>1>. In other words, holds for () if and only if it holds for (). We conclude that () is a total unravelling w.r.t. () if and only if () is a total unravelling w.r.t. (). One of the easiest examples which illustrates the importance of unravellings is P(f)=f2>--1>>*f+-1>>=e,> where z\e>. This equation admits as it's unique potential dominant term of multiplicity . However, the refinement (\1)> leads to the equation 2>--1>|1-z-1>>*+-1>|1-z-1>>=e(\1),> which again admits a unique potential dominant term > of multiplicity . Continuing like this leads to an infinite sequence (\1)>, =z-1>+|~>(|~>\z-1>),\> of refinements, and we do not succeed in separating the two solutions. Therefore, we rather construct a total unravelling. In order to do so, we first compute the distinguished solution =(1-z-1>)> to the quasi-linear differential equation )= P|\ f>(\)=2*\--1>>*\=0> and then refine +(\1).> This refinement (and partial unravelling with associated monomial ) is actually already the total unravelling we were looking for and the equation () transforms into =e(\1).> This time, the new equation admits two potential dominant terms > and > of multiplicities 2>, which allows us to compute the solutions to () by recursion. In general, total unravellings can only be achieved via successions of partial unravellings. An important example which illustrates this phenomenon is the following: f2>+2*f+>+*log z>+\+*log z*\*log z>=0> This equation becomes purely exponential after upward shiftings: f2>+*e>*\*exp z>*f+ z>+ z*exp z>+\+*e>*\*exp z>=0.> This new equation admits a unique potential dominant monomial z> of multiplicity . Indeed, this is easily seen when substituting z)*g> for : z>*g+*\*exp z>*g-2*g+1+ z>+\+*\*exp z>=0.> Now (\1)> is a partial unravelling w.r.t. () and z+(\exp z)> a partial unravelling w.r.t. (). However, the partial unravelling transforms () into an equation of the same form as (), but with decreased by . Consequently, we need a succession of unravellings || z>>+ff\ z>>;>>|>|| z*exp z>>+ff\ z*exp z>>;>>||>|>|>||*\*exp z>>+ff\*\*exp z>>.>>>>> in order to attain the total unravelling z>+ z*exp z>+\+*\*exp z>+f\*\*exp z>.> Notice that the ratios z,\,e> of the successive potential dominant terms indeed satisfy proposition . An open question is whether there exist examples which essentially need the technique of adjusted partial unravellings in order to limit the appearance of iterated logarithms. If not, this would allow some major simplifications in the algorithm in the next section. Actually, the whole computation process of total unravellings needs to be better understood in order to generalize it to other fields of transseries (see section ). In this section, we give algorithms to solve an asymptotic differential equation (), where the coefficients of and > can be expanded with respect to a parameterized transbasis =(\,\,\)>. Actually, we show how to solve such an equation modulo a transmonomial >. Here a >> to () is a transseries \\>, such that the Newton degree of >()=0(\\)> is strictly positive. In our algorithms, we will use the following conventions without further mention. Our algorithms are non deterministic in the sense that we allow automatic case separations, in a similar way as sections and . The algorithms are constructed in such a way that each solution modulo > to () occurs in precisely one case. In other words, for each specialization of the initial parameters and, for each specialization of the directions (which corresponds to a choice of a field of transseries as in section , which satisfies the constraints on the >>) and for each solution > in , there exists precisely one case and precisely one specialization of the newly introduced parameters, which yields > as a solution. The termination of our non deterministic algorithm follows from the termination of each of its branches by a similar Noetherianity argument as in the proof of theorem . For more details about the automatic case separation strategy, see chapter 8 in . We assume throughout our algorithms that all computations are done w.r.t. a purely exponential parameterized transbasis =(\,\,\)>. In order to do this, we are allowed to insert new elements into > and to ``shift the whole computation upwards'' whenever necessary. Upward shifting can for instance be implemented via an exception, which starts over the whole non deterministic computation, after shifting the input upwards. A better strategy, which consists of associating an ``exponential level'' to each transseries, is explained in more detail on page 276 of . It is useful to allow a few non standard transmonomials in our algorithms, like the formal infinitely large and small monomials >> resp.>>>, as well as monomials of the form z*exp z*\*z*log z*\)>. Of course, our algorithms are not really effective, as long as we do not provide algorithms to compute the expansions of parameterized transseries and to test them for being zero. In this paper, we assume that we have oracles for deciding such questions. Similarly, in our algorithms for computing distinguished solutions in the next section, we will merely give infinite formulas for the results. Nevertheless, we notice that for certain classes of transseries, the oracles may be replaced by real algorithms, in a similar way as described in chapter 12 of . In order to get a parameterized version of theorem , we have to make sure that the operator =\*(Id+R*\)> is well defined. This can be done as follows: <\algorithm|linear> A linear differential operator and transseries The distinguished solution to <\body> [Introduce generic monomial] <\indent> Let =(\,\,\)> denote the current purely exponential transbasis Let ,\,\> be temporary new parameters in and set \\>*\*\>> [Compute the distinguished solution] <\indent> Uniformly regularize *L\>> Let \ (Id+R*\) g>, with the notations from the proof of theorem [Clean up] <\indent> Destroy the parameters ,\,\> by projection of the regions Return In the algorithm, the uniform regularization of a differential operator or polynomial means that we uniformly regularize all its coefficients and that we make the corresponding dominant monomials pairwise comparable for > using theorem and lemma . In the parameterized context, the equation () has a Newton degree , if each specialization of the directions and parameters leads to an equation of Newton degree . We will show below how to compute the Newton degree modulo case separations. The distinguished solution to a quasi-linear equation (i.e. of Newton degree ) is computed inductively as in the proof of theorem . The dominant part of w.r.t. ,\,\> is computed by the usual formula after the uniform regularization of . <\algorithm|quasi_linear> ,k\n)> {0,\,n}> (with as default value), a differential polynomial with coefficients in ;\;\>> and a monomial \\*\*\>, such that () is quasi-linear> The distinguished solution to () <\body> [Normalize] <\indent> If \1>, then return > times \>,1,k)> Uniformly regularize If (P)\1> then return (P)*P,1,k)> [Recurse] <\indent> Compute the dominant part of w.r.t.>> Let \(D,1,k-1)> [Return] (Id+R*\) >, with the notations from the proof of theorem > The -equalizer of a differential polynomial is computed similarly as in the proof of proposition , by uniformly regularizing > and > at each recursion. <\algorithm|equalizer> A differential polynomial and integers i\j\deg P> The -equalizer for <\body> [regularize] > and >> [equalize] <\indent> Let \(P)/\(P)|j-i>> If \1> then return *(P\>,i,j)> [shift upwards] <\indent> If (P)=\(P))> and +P>\C[c]*(c)>>, then return Shift upwards and return to step 2 The analogue of the Newton polygon associated to an algebraic differential equation in the differential case is the determination of -equalizers, which occur as potential dominant monomials, and which are extremal in the sense that is maximal. These equalizers correspond to the slopes of the Newton polygon and the and to the first coordinates. <\algorithm|Newton_polygon> A differential polynomial \\\i=deg P> and potential dominant monomials \\\\> for , such that > is the ,i)>-equalizer for for each > <\body> [Initialize] <\indent> val P> 0> [Insert vertex] <\indent> \j> If then return ,\,i> and ,\,\> k+1> [Search edge] <\indent> Compute >=(P,j,j)|j-j>> for all \j> with >\0> Let \min> {\>}> and let > be minimal with =\>>. Set j> and go to step 2 The Newton degree of () can easily be read from the Newton polygon: <\algorithm|Newton_degree> )> A differential polynomial and a monomial > The Newton degree of (f\\)> ,\,\> by |,\,i>>> Let {0,\,k}> be maximal, such that \\> Return > Given the Newton polygon associated to , let us now show how to determine the potential dominant terms of solutions to () and their multiplicities. We separate a case for each edge and each vertex of the Newton polygon, and determine the corresponding algebraic or mixed resp. differential potential dominant monomials and terms. In order to determine the differential potential dominant monomials, we recursively have to solve Ricatti equations modulo )>. The algorithm which does this will be specified in the next section. <\algorithm|pdt> )> [non deterministic] A differential polynomial and a monomial > A potential dominant term > for () <\body> [Determine Newton degree] <\indent> ,\,\> by |,\,i>>> Choose a {0,\,k}>, such that or \\>, and set i> If then go to step 3 Separate two cases and go to step 2 resp. 3 [Algebraic and mixed potential dominant terms] <\indent> Let \\> Let be a new parameter in Impose the constraint \>>(c)=0> (as an algebraic constraint, since =0>) Return > [Differential but non mixed potential dominant terms] <\indent> Let (R,\>,(x*log x*log log x*\))> Let \ exp g>, where the integral is computed using If k> then impose the constraint \\> Otherwise impose the constraint \\> If 0> then impose the constraint \\> Let be a new parameter in and impose the constraint 0> Return > <\algorithm|multiplicity> )> A differential polynomial and a term > The multiplicity of )> as a root of \(\)>>> <\body> Repeat <\indent> Uniformly regularize \(\)>> If \ \(\)>>\C[c]*(c)>> then shift upwards Until \(\)>>\C[c]*(c)>> Return the multiplicity of >> as a root of \(\)>>> We can now state the main resolution algorithm for solving asymptotic algebraic differential equations () modulo monomials >. <\algorithm|solve> ,\)> A differential polynomial and monomials > and > A solution to (f\\)> modulo > <\body> [Initialize] <\indent> \0> \> [Are we done?] <\indent> If >,\)\0>, then separate two cases and respectively <\indent> 1. Return > 2. Proceed with step M3 [Compute potential dominant term] <\indent> Let (P>,\)> Let \(P>,\)> If \\> then kill this process [Does > have multiplicity d>?] <\indent> If (P,\)\d> then <\indent> \\+\> \\(\)> \> Go to step M2 [Dispatch on ] <\indent> If => then go to step H1 If => then go to step H3 If => then go to step U4 [Prepare first step unravelling loop] <\indent> \P> \\(\)> \> \> [Prepare partial unravelling loop] <\indent> While ,\\>>\C[c]*(c)>> shift upwards If ,\\>>\C[c]*(c)> with d> then (\d-1> \,\\>/\ fd-1-k>*\ (f))\,-\>> Otherwise (\d-1> \,\\>/\ (f)d-1>)\,-\>> [Partial unravelling] <\indent> If (Q>,\)=1> then <\indent> \(Q+\>,\)> \\+\+\> \\(\)> Go to step M2 [Dispatch on ] <\indent> If => then go to step T1 If => then go to step T2 [Prepare other steps unravelling loop] <\indent> Let be such that \\/\(\)> Uniformly regularize ,\\(\)>> Compute the dominant part > of ,\\(\)>> w.r.t.>> and set \\(\),-\>> \z> \> [Prepare next step unravelling loop] <\indent> If there is no purely exponential monomial > in > with (\)\\r}>*\>, then set \z> Let > be a purely exponential monomial in >, such that (\)\\,r}>*\> \> [Adjust partial unravelling] <\indent> If \\r}>*\> then <\indent> \\+\> \\(\)> Go to step M2 If \\> then \\(\)> \> Go to step H3 The algorithm gradually constructs a solution > modulo > to () via a succession of refinements. Each time we get back to the main entry of the loop, we actually have to solve the equation >(f)=0(f\\)>. Given a potential dominant term > for this equation, the next refinement (and value of >) depends on the variable. The core of the algorithm consists of steps in which case =>. As long as we do not hit a potential dominant term > of maximal multiplicity , the algorithm only executes steps . Given a potential dominant term > of multiplicity d>, we can simply take =\+|~>(|~>\\)> for our refinement, which corresponds to the assignments \\+\> and \\(\)> in step M4. The steps correspond to the first partial unravelling when we hit a potential dominant term > of maximal multiplicity . As long as we do not enter , we will have =P>, => and =>. We start by computing the differential polynomial from section (modulo an additive conjugation by >). We then keep refining \\+\+|~>(|~>\\)> as far as possible in step H3, where > is the distinguished solution to the quasi-linear equation +\>(\)(\\\)>. If the steps do not lead to a complete unravelling, we apply the theory from sections and in steps . We start by computing once and for all the differential polynomial > in , with this change that we apply a multiplicative and additive conjugation to it in order to make it ``compatible'' with . The transseries >, which is initialized with , may become an iterated exponential z> as a result of upward shiftings. As long as the current potential dominant term > has not yet the required form in order to start a partial unravelling, we have =>, and we keep on adjusting in step T3. The termination of is guaranteed by propositions and , modulo the hypothesis that the resolution process requires only a finite number of upward shiftings. An upper bound for will be given in the next section. The following main theorem describes the general form of solutions to asymptotic algebraic differential equations () with parameterized complex transseries coefficients. <\theorem> Consider an asymptotic algebraic differential equation )> with transseries coefficients. Then, modulo case separations, there exist a finite number of parameterized transseries solutions ,\,f> to )>, with the following properties: <\enumerate-alpha> The logarithmic depths of ,\,f> do not exceed the logarithmic depths of the coefficients of by more than a fixed constant \d*(4*w)>, which only depends on the Newton degree , the order and the weight P\<\|\|\>> of )>. For each specialization of the parameters occurring in the coefficients of , for each specialization of the directions, and for each solution to )> after these specializations, there exists exactly one > and exactly one specialization of the remaining parameters on which depends >, for which > specializes to . <\proof> In view of the algorithm from the previous section, we only have to prove (). We prove the bounds for > by a double induction over and . For , we necessarily have , and clearly =0>. So assume that 0>. If , then () has no solutions, so that =0>. Assume therefore that 0>. We first observe that the number of upward shiftings needed to compute a potential dominant term is bounded by \max {w,U}> by propositions (), and the induction hypothesis. We have to estimate the maximal number of upward shiftings which may occur in the main loop of , before we reach a lower Newton degree. Now the first partial unravelling requires at most +r> upward shiftings, in view of remark . By the induction hypothesis and proposition , the second loop of adjusted partial unravellings requires at most +U+1> upward shiftings. Finally, the decisive refinement, which decreases the Newton degree, again needs at most > upward shiftings. Altogether, we obtain \U+3*T+U+r+1.> Consequently, \d*(3*T+U+r+1)> In particular, for , we obtain \U\d*(3*d+2).> For 1>, we obtain \d*(4*U+r+1).> By induction, we finally notice that \(4*w)*(w+1),> which implies our bound. <\remark> Of course, it is possible to improve the bound for > for particular values of and . First of all, in the case when , it is easily checked that upward shifting is sufficient in proposition (), so that \1.> This observation implies the sharper bounds >|>|>|>|>|>>>>> for >. A careful analysis of the differential Newton polygon method will probably lead to even sharper bounds for small values of . Similarly, it is possible to improve the bounds for small values of , by using the fact that the weight of > is bounded by r>> for d>. Although the above theorem describes the general form of solutions to (), it does not claim the actual existence of such solutions. We say that > is a solution of multiplicity > of (), if the differential valuation of >> equals >. The following theorem stipulates the existence of solutions to () of a very special form. Consider an asymptotic algebraic differential equation )> of Newton degree and whose coefficients can be expanded w.r.t. a transbasis ,\,\)>. Then there exist at least solutions to )> when counting with multiplicities. Moreover, these solutions can be expanded w.r.t. \,\,log \,\,\,\)> for some .> <\proof> Without loss of generality, we may assume that =e>. Let us prove the theorem by induction over . For we have nothing to prove. For , the equation is quasi-linear and the distinguished solution can be expanded w.r.t. z,\,z,\,\,\)>. Assume therefore that 1>. If there exists only one algebraic potential dominant term with multiplicity , then consider the unravelling +(\)> we obtain by executing , but where we always choose the unique algebraic potential dominant term in . Since this branch only involves the computation of equalizers and solutions of quasi-linear equations, > can be expanded w.r.t. z,\,z,\,\,\)>. Modulo replacing by >>, we may thus assume without loss of generality that () admits no algebraic potential dominant terms of multiplicity . If there exists a mixed potential dominant monomial >, then > is a potential dominant term of multiplicity d> for each 0>, and the coefficients of >> can be expanded w.r.t. P\<\|\|\>> z,\,z,\,\,\)>. By the induction hypothesis each equation >()(\\)> admits at least one solution which can be expanded w.r.t. z,\,z,\,\,\)> for some . Hence, there exists an infinity of solutions with the required properties. In what follows, we therefore assume that all potential dominant monomials are algebraic, but not mixed. Now let i\\\i\d> and \\\\> be such that > is the ,i)>-equalizer for each {1,\,s}>. For each {1,\,s}> the Newton polynomial \>>> is a polynomial with valuation > and degree >, which has -i> roots (when counting with multiplicities). These root induce at least > potential dominant terms, which can be expanded w.r.t. P\<\|\|\>> z,\,z,\,\,\)>, and whose multiplicities are d>. By proposition and the induction hypothesis, this leads to at least > solutions of the required form, when counting with multiplicities. The theorem now follows from the fact that is a solution of multiplicity >. We recall that a differential field is said to be , if for any pair of differential polynomials over , such that the order of is strictly larger than the order of , there exists a root of in , which is not a root of .\ Let be a field of complex transseries as in section . Unfortunately, theorem is not sufficient for to be differentially algebraically closed. Indeed, the only transseries solutions for the elliptic equation +(f)+f=0> are , and . Consequently, there are no transseries solutions to this equation, which are not solutions of the equation of lower order +f=0.> Nevertheless, theorem is sufficient for the following application. <\theorem> Let be a linear differential operator of order with coefficients in . Then <\enumerate-alpha> can be completely factored over . There exist linearly independent solutions to in . <\proof> By theorem the Ricatti equation associated to has at least one solution \\>. Consequently, we may factor *(\-\),> for some linear differential operator > of lower order and coefficients in . Part () now follows by induction over . Now consider a factorization -\)*\*(\-\),> with ,\,\\\> and let >||\>;>>|>||-\) e\>;>>||>|>|>||-\)*\*(\-\)] e\>,>>>>> where > stands for distinguished integration. Then =\=L h=0>. Moreover, by the distinguished properties of the left inverses -\),\,[(\-\)*\*(\-\)]>, we have (h)>=0> for all j>. This guarantees the linear independence of ,\,h>. Indeed assume that there exists a relation *h+\+\*h=0> with \0>. Then *h+\*h+\+\*h)(h)>=(\*h)(h)>=\\0>. This contradiction completes the proof of (). in such a way that \\\h>, we even obtain the canonical basis of solutions of in the proof of theorem .> <\remark> In the case of real transseries, it can be shown that each linear differential operator may be factored as the product of a transseries and operators of the form |\ z>+a> or |\ z>+2*a-|b>*|\ z>+a+b+a-|b>=|\ z>+a-b*i+|b>*|\ z>+a+b*i.> \ Although the algorithm provides us with the generic solution to (), it is not clear priori> that the number of new parameters on which the solution depends does not exceed . In this section we sketch a proof of the fact that the number of such integration constants is indeed bounded by . We first notice that the only place where we introduce (continuous) integration constants is in step 3 of . Each integration constant can therefore be ``attached'' to a solution of a Ricatti equation of the form \>>. Given an arbitrary moment during the algorithm , we actually search solutions of the form *e\+\*e\+>>>,> where ,\,\> are the ``active integration constants''. The idea is now to set =\*e\+\*e\+>>>,F=\*e\+\*e\+>>>,\,F=\*e\+>> and to consider as a differential polynomial of order in >, with coefficients in ;\;\,F,\,F>>. In other words, we consider ,\,F> as new monomials and we give *\*\*FC>*\*FC>> the natural ``pointwise'' quasi-ordering > (see chapter 6 of ). The only obstruction for the computation with coefficients in ;\;\,F,\,F>> instead of coefficients in ;\;\>> is when the uniform regularization of a transseries in ;\;\,F,\,F>> is not possible. Now this obstruction corresponds to the imposition of an algebraic constraint on an active integration constant >, when performing the same computation in ;\;\>>. In order to solve this problem, an ``error handler'' is installed each time that we introduce a new continuous integration constant >. Whenever we impose an algebraic constraint on >, we go back to the error handler and reperform the same computations while assuming that > either did or did not (non determinism) satisfy the algebraic constraint right from the start. In all branches of the new resolution process, the order of the asymptotic differential equation, when rewritten as an equation in >, does not exceed . Consequently, r> at the end of each branch of the process. The reader may have noticed a certain number of changes with respect to the treatment of algebraic different equations in . Although the results of this paper were stated in the context of grid-based transseries, they may easily be adapted to the well-ordered context from , except for the results about parameterized transseries, which become more complicated. The algorithm may still be applied in the well-ordered context, except that the introduction of new parameters should then be interpreted as a new source of (continuous) case separations. During a careful reexamination of our previous work, we noticed that proposition 5.7() in does not hold for all . Consequently, our previous treatment of almost double solutions in section 5.5.1 does not work. The present, more complicated, treatment using unravellings corrects this error. When calling a refinement occurring in our construction of a total unravelling a , the proof of theorem 5.2 in remains correct (except for the bound for the maximal length of a chain of privileged refinements, which may have to be replaced by a larger bound). Some other changes with respect our previous work are the following: <\itemize> In view of theorem 3.3 in it is no longer necessary to develop the theory from section in the purely exponential setting first (as we did in ). We simplified and improved the construction of distinguished solutions to linear and quasi-linear equations, through a new application of the generalized implicit function theorem from . In comparison with the effective asymptotic resolution \ of algebraic differential equations in chapter 12 from , we noticed that we actually never need to impose exponential constraints on the parameters. After correcting the error related to privileged refinements, we therefore no longer need to assume the existence of an oracle to determine the consistency of first order systems of exp-log constraints in theorem 12.4. In the corrected version, we also consider the case when \>>=\*(c-\)*(c)> with k\d>. We forgot that case in the original version. In this paper, we have generalized the transseries technique for solving algebraic differential equation as far as reasonably possible. Three main problems remain to be solved. First of all, one has to show that the complex transseries solutions to algebraic equations have a genuine analytic meaning. This problem, which will be treated in a forthcoming paper, can actually be subdivided in two parts: <\itemize> We have to show that a consistent system of asymptotic constraints on the directions corresponds to a non empty asymptotic region of the complex plane. In general, this region does not need to be connected. We have to give an analytic meaning to our transseries solutions on regions as above. This analytic meaning should be compatible with the asymptotic relations, which have in particular to be formalized on disconnected regions. We have already remarked that our fields of complex transseries are not differentially algebraically closed. In other words, we still miss most of the solutions to algebraic differential equations in our formalism. In order to get a full understanding of the asymptotic behaviour of solutions to algebraic differential equations, two approaches may be followed: <\itemize> In order to solve an equation like +(f)+f=0> one may start with studying the solutions in the neighbourhoods of singularities other than >. This can be done by performing a change of variable >, which transforms the equation into an equation which does admit a solution space of dimension . More generally, for a general asymptotic algebraic differential equation (), the above trick leads to transseries solutions in >\>>, where > is the original variable, and >||+\*z;>>||>|>|>||+\*z.>>>>> It is not yet clear to us how to alternate usual refinements with substitutions of the form >. Assuming for simplicity that \ has constant coefficients, one may also start with studying the singularities of the dynamical system associated to the algebraic differential equation. For instance, one may use the theory from chapter 10 in to desingularize as a polynomial in ,\,f>. This leads to a better understanding of the behaviour of the dynamical system for different subregions of ``the ,\,f)>-space''. We next apply the asymptotic tools from this paper to obtain full solutions on these regions. Finally, one has to study how the solutions globally glue together. In any case, a purely local treatment seems not to be possible in order to describe all solutions to an algebraic differential equation. A good combination of a more global theory with our local results might lead to the resolution of interesting questions, such as to admit a natural boundary somewhere on its Riemann surface?> For Liouvillian functions and, in view of the theorem , for functions which are obtained via the repeated resolution of linear differential equations, the answer seems to be negative. In relation to a joint project with Aschenbrenner and Van den Dries, which aims to describe the model theory of ``real differential fields'' and ``valuated differential fields'', it seems important to better understand our technique of unravellings in cases where the transseries in do not necessarily admit finite logarithmic depths. A typical example of an equation which is ``hard to unravel'' is 2>+2*f+>+*log x>+*log x*log log x>+\=0.> A good question is whether there are essentially different examples of equations which are hard to unravel. Another question is whether we may avoid the adjusted partial unravellings from section . <\bibliography|bib|alpha|~/publs/all.bib> <\bib-list|[99]> J. Écalle. . Hermann, collection: Actualités mathématiques, 1992. J. van der Hoeven. . PhD thesis, École polytechnique, France, 1997. Joris van der Hoeven. A differential intermediate value theorem. Technical Report 2000-50, Univ. d'Orsay, 2000. Joris van der Hoeven. Operators on generalized power series. , 2000. Submitted. <\initial> <\collection> <\references> <\collection> > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > <\auxiliary> <\collection> <\associate|bib> Ec92 vdH:phd VdH:ivt vdH:phd vdH:phd vdH:phd vdH:phd D82 vdH:phd vdH:phd vdH:phd VdH:noeth VdH:noeth vdH:phd vdH:phd vdH:phd vdH:phd vdH:phd vdH:phd vdH:phd VdH:ivt VdH:ivt vdH:phd VdH:noeth vdH:phd vdH:phd <\associate|figure> Shape of the region |P>> for |\>=1> resp. |\>=|-1>>.|> |f=log e+i*z>+e>>, which illustrates the four possible asymptotic behaviours of |f> on non degenerate regions. The ``rows'' of singularities correspond to the borders between regions of different types.|> <\associate|toc> |math-font-series||1Introduction> |.>>>>|> |math-font-series||2Complex transseries> |.>>>>|> |2.1Real trigonometric fields |.>>>>|> > |2.2Series with complex coefficients and monomials |.>>>>|> > |2.3Pre-fields of complex transseries |.>>>>|> > |2.4Logarithmic complex transseries in |z> |.>>>>|> > |2.5Complex transseries in |z> |.>>>>|> > |2.6Fields of complex transseries |.>>>>|> > |2.7Extra structure on the field of transseries in |z> |.>>>>|> > |2.8Further generalizations |.>>>>|> > |math-font-series||3Generic complex transseries> |.>>>>|> |3.1Generic complex transbases |.>>>>|> > |3.2Case separations and the field operations |.>>>>|> > |3.3Logarithms of complex transseries |.>>>>|> > |3.4Exponentials of complex transseries |.>>>>|> > |Case 1: |.>>>>|> > |Case 2: |.>>>>|> > |Case 3: |.>>>>|> > |3.5A worked example |.>>>>|> > |Expansion of |e>. |.>>>>|> > |Expansions of |i*z> and |e+i*z>. |.>>>>|> > |Expansion of |e+i*z>>. |.>>>>|> > |Expansion of |e>>. |.>>>>|> > |Expansion of |f>. |.>>>>|> > |math-font-series||4Parameterized complex transseries> |.>>>>|> |4.1Definition of parameterized complex transseries |.>>>>|> > |4.2Uniform regularization |.>>>>|> > |4.3Computations with parameterized complex transseries |.>>>>|> > |Case |e\1>. |.>>>>|> > |Case |e\1>. |.>>>>|> > |math-font-series||5The differential Newton polygon method> |.>>>>|> |5.1Notations |.>>>>|> > |5.1.1Asymptotic relations |.>>>>|> > |5.1.2Natural decomposition of |P> |.>>>>|> > |5.1.3Decomposition of |P> along orders |.>>>>|> > |5.1.4Logarithmic decomposition of |P> |.>>>>|> > |5.1.5Additive and multiplicative conjugations and upward shifting. |.>>>>|> > |5.2Differential Newton polynomials |.>>>>|> > |5.3Potential dominant monomials and terms |.>>>>|> > |5.3.1Algebraic potential dominant monomials |.>>>>|> > |5.3.2Differential potential dominant monomials |.>>>>|> > |5.4Refinements |.>>>>|> > |5.5A worked example |.>>>>|> > |math-font-series||6Distinguished solutions> |.>>>>|> |6.1Distinguished left inverses of linear differential operators |.>>>>|> > |6.2Distinguished solutions of quasi-linear equations |.>>>>|> > |6.3A worked example |.>>>>|> > |math-font-series||7Unravellings> |.>>>>|> |7.1Total unravellings |.>>>>|> > |7.2Partial unravellings |.>>>>|> > |7.3Adjusted partial unravellings |.>>>>|> > |7.4Construction of total unravellings |.>>>>|> > |7.5Worked examples |.>>>>|> > |math-font-series||8Computation of parameterized solutions> |.>>>>|> |8.1About the algorithms in this section |.>>>>|> > |Automatic case separations. |.>>>>|> > |Automatic upward shiftings. |.>>>>|> > |Non standard monomials. |.>>>>|> > |Effective computations with transseries. |.>>>>|> > |8.2Distinguished solutions |.>>>>|> > |8.3Determining the Newton polygon |.>>>>|> > |8.4Computing potential dominant terms |.>>>>|> > |8.5Solving the differential equation |.>>>>|> > |math-font-series||9Main theorems and final remarks> |.>>>>|> |9.1Complex transseries solutions to algebraic differential equations |.>>>>|> > |9.2Linear differential equations |.>>>>|> > |9.3Bounding the number of integration constants |.>>>>|> > |9.4Comparison with previous work and errata |.>>>>|> > |9.5Conclusion and perspectives |.>>>>|> > |Analytic counterpart. |.>>>>|> > |Differentially algebraic closure. |.>>>>|> > |Unravellings. |.>>>>|> > |math-font-series||Bibliography> |.>>>>|>