def prde_spde(a, b, Q, n, DE):
    R, Z = list(zip(*[gcdex_diophantine(b, a, qi) for qi in Q]))

    A = a
    B = b + derivation(a, DE)
    Qq = [zi - derivation(ri, DE) for ri, zi in zip(R, Z)]
    R = list(R)
    n1 = n - a.degree(DE.t)

    return (A, B, Qq, R, n1)
Special Polynomial Differential Equation algorithm: Parametric Version.
dataset/ETHPy150Open sympy/sympy/sympy/integrals/prde.py/prde_spde