5 Notation mapping between Saunders/Smith algorithm and original Kovacis algorithm

I have implemented the original Kovacis algorithm using Maple 2021 based on the original paper (1). The following are notation difference between the two algorithms and the implementation by Smith [3] that I found.

  1. Kovacis algorithm  uses α± defined as α±=12±121+4b for the case when O()=2. Smith algorithm uses e0 for the 1+4b part only. In both algorithms the b value is calculated in the same way. It is the coefficient of 1x2 in the Laurent series expansion of r at . But we do not need to find Laurent series expansion of r at to find b here. It can be found using b=lcoeff(s)lcoeff(t) where r=st and gcd(s,t)=1.
  2. Smith algorithm finds e1,e2, values for each pole. This is part b of step 1 for poles of order 2, these correspond to only the 1+4b part in Kovacis algorithm (this is part c2 of step1), where there it finds [r]c for each pole and αc±=12±121+4b where b is the coefficient of 1(xc)2 in the partial fraction decomposition of r. This b value is also the same for Smith algorithm in its es.

More mappings to be added next.