@flacjacket : For the example given in the Maple Paper, that is, `eq = x*log(x)*sin(f(x))*f(x).diff(x) + cos(f(x))*(1 - x*cos(f(x)))` I get a basis of 94 factors. If I assume the coefficients, to be constant, that is the factor itself to be `chi` it takes around 10 seconds to figure out that `chi` is `[cos(y)**2/(log(x)*sin(y))]` , and that is when I just substituted the factors in the PDE `chix + h*chiy - hy*chi` . It would increase obviously when I use solve and try to group for around 90+ factors and solve for the coefficients. And when I increase the degree such that the coefficients are bivariates. I thought there might be a workaround, but nothing is striking me. What are your thoughts?