The Asymptotic Expansion Method via Symbolic Computation

Algorithm 4

Implementation of the algorithm. Steps (3.1) and (3.2).

xdx [𝟶][𝟶][𝟶].add            Initialize 𝑋(0,0,0)=1.

(𝟷.,𝟶.,𝟶,𝟶.,𝟶.);

xdx [𝟷][𝟶][𝟶] = u.solvePVI       Compute 𝑋(1,0,0) as the solution to (5.2),

(a1, a0, t0, x0, dx0);       and 𝑋(0,1,0)=𝑑/𝑑𝑡(𝑋(1,0,0)).

xdx [𝟶][𝟷][𝟶] = xdx [𝟷][𝟶][𝟶].der ();