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^2*D^2) o2 = cokernel | x2D2 | 1 o2 : W-module, quotient of W |
i3 : PolyExt(M) 2 o3 = HashTable{0 => QQ } 2 1 => QQ o3 : HashTable |