Factorization of KdV Schr\"odinger operators using differential subresultants

We address the classical factorization problem of a one dimensional Schr\"odinger operator $-\partial^2+u-\lambda$, for a stationary potential $u$ of the KdV hierarchy but, in this occasion, a "parameter" $\lambda$. Inspired by the more effective approach of Gesztesy and Holden to the "direct" spectral problem, we give a symbolic algorithm by means of differential elimination tools to achieve the aimed factorization. Differential resultants are used for computing spectral curves, and differential subresultants to obtain the first order common factor. To make our method fully effective, we design a symbolic algorithm to compute the integration constants of the KdV hierarchy, in the case of KdV potentials that become rational under a Hamiltonian change of variable. Explicit computations are carried for Schr\"odinger operators with solitonic potentials.