The algorithm used appears in the paper 'Polynomial and rational solutions of holonomic systems' by Oaku-Takayama-Tsai (2000). The method is to combine isomorphisms of Bjork and Kashiwara with the restriction algorithm.
i1 : W = QQ[x, D, WeylAlgebra=>{x=>D}] o1 = W o1 : PolynomialRing |
i2 : M = W^1/ideal(x*D+5) o2 = cokernel | xD+5 | 1 o2 : W-module, quotient of W |
i3 : RatExt M 1 o3 = HashTable{0 => QQ } 1 1 => QQ o3 : HashTable |