The DESIR Package

Differential Equations Singularities Irregular and Regular

The DESIR package contains commands in Maple that help you solve linear ordinary differential equations in the neighborhood of singularities in the complex plane.

Jean-Guillaume Dumas

Formal computation and manipulation of formal series

The first part consists in a new version of the DESIR-II package, which was written in Maple V [1, 5]. Both versions are improved forms of the DESIR-I package written in Reduce [2]. The purpose of this part is the computation of formal solutions of homogenous linear ordinary differential equations. In order to have access to all information the formal solutions contain, the internal data structure can be assigned to an optionally passed argument. The Gevrey caracteristics (type and order) of the divergent series are also stored in this data structure. Specific information can then be extracted by using the primitive functions (irregpart, regpart,param).

Jean-Guillaume Dumas

Numerical computation and visualization

The second part collects some functions for the numerical computation and the graphical visualization of the solutions in the complex plane.
The principle of the representation is the following [6]: it consists in plotting the image under the considered function f of a circle or a circular arc around the singularity, in general 0 (or infinity). The color is used to associate a point in the domain and its image: each point f(x) is plotted with a color corresponding to the argument of x.
As the studied functions are in general multi-valued, we consider them in the neighborhood of 0 as functions on the Riemann surface of the logarithm and points in the domain are represented by their Euler coordinates, whereas the image points are computed in cartesian coordinates.

Jean-Guillaume Dumas

Computation of Stokes matrices

The third part is new. It deals with the Stokes phenomenom, in the neighborhood of an irregular singular point. The main function is StokesMatrices, which describes this phenomenom as a list of objects [arg, M_arg]: arg is the argument of a Stokes ray and M_arg is the corresponding Stokes matrix. Indeed, the Stokes constants can be numerically computed by the analysis of the singularities of the Borel transform of the divergent series, and this for a large class of differential equations of single rank k >1: the function works under the hypothesis that the Borel transforms don't have aligned singularities in the Borel plane, and it allows  there any polar, ramified or logarithmic singularities [4]. The function monodromy computes the matrix of formal monodromy. Then it is possible to use the knowledge of these matrices to define a function (in Maple) which calculates automatically, for a particular solution f, a suitable linear combination, depending on the sector. This is the aim of the function combli_variable.

Jean-Guillaume Dumas



The directory desir_2016 contains two maple files,
                       maple.mla and,
that you have to download in your own directory "mydirectory".
In a maple session, you have first to add this directory to the global variable libname:
and then to load the package: with(Desir);

The command ?Desir will give you an overview of the package. All the functions that are accessible to the user are documented.
To run the package, you need: Maple at least 10.0

Jean-Guillaume Dumas


Françoise JUNG
Laboratoire Jean Kuntzmann
38401 Domaine universitaire de Saint Martin d'Hères, UGA, France.