The Asymptotic Expansion Method via Symbolic Computation

Algorithm 3

Implementation of the algorithm. Steps (1) and (2).

MQ xdx [𝙼+𝟷][𝙼+𝟷][𝙾𝚁𝙳𝙴𝚁+𝟷];          Define a three-dimensional array of quasi

                 polynomials𝑋(𝜌₁,𝜌₂,𝑞), where 𝜌₁,𝜌₂,𝑞∈

                 𝑁, 0≤𝜌₁,𝜌₂≤𝑀 and 0≤𝑞≤𝑄, being

                 𝑄 the order of the asymptotic expansion.

MQ *s;                 Define an array of quasipolynomials 𝑥(𝜌),

s=new MQ [𝙾𝚁𝙳𝙴𝚁+𝟷];          where 0≤𝜌≤𝑄.

int 𝚒,𝚓,𝚔,𝚙,𝚚,𝚗𝚞;         Definition of auxiliary variables.

MQ U;