Introduction
A multivariate formal series f∈Q[[x1,…,xn]] can be written as f(x1,…,xn)=∑i1,…,in≥0c(i1,…,in)xi11⋯xinn. When 1≤k≤n denote fk the partial derivative of f with respect to the variable xk, and when 1≤k1,…kt≤n denote fk1⋯kt the order t partial derivative with respect to the variables xk1,…,xkt. Since two formal power series are equal if all the coefficients of their monomials are equal, a partial differential equation for f (formal PDE) is equivalent to a system of polynomial equations for the coefficients c(i1,…,in).
Example of mixed exponential type
Consider n=2 variables, and the formal power series f(x1,x2)=∑d≥1N(d)edx1x3d−12(3d−1)!. This is an example of formal series of "mixed exponential type", because its presentation involves exponentials of some of the variables. By expanding the exponential and renaming indices, it can be written as above with c(i1,i2)=(i2+13)i1N(i2+13)/(i1!i2!) if i2≡2(mod 3) and c(i1,i2)=0 otherwise. The PDE f222=f2112−f111f122 is equivalent to the polynomial equations for d≥2 N(d)=\sum_{d_1+d_2=d}N(d_1)N(d_2)\left( d_1^2d_2^2 \binom{3d-4}{3d_1-2} -d_1^3d_2\binom{3d-4}{3d_1-1}\right) and this sequence is uniquely determined by the value of N(1). As discovered by Kontsevich, if N(1)=1 then N(d) is the number of rational curves of degree d through 3d-1 generic points in the complex projective plane.
Question
Can one use SageMath to solve formal PDEs? By solving I mean the following: given the formal PDE and a degree bound b \in \mathbb{N}, compute all numbers c(i_1,\ldots ,i_n)\in\mathbb{Q} with i_1+\ldots +i_n\leq b for the solution (assuming that it exists and is unique). In practice, it would be nice if the proposed method allowed to work with formal series of "mixed exponential type" too. A good test would be to see if in the example above the method computes, given b\in\mathbb{N}, all numbers N(d) with d \leq b as listed in OEIS A013587 (cannot post link due to low karma).
