using polynomial power series for solving PDEs

[Foad]

Hello everyone, 

While ago I was trying to use polynomial power series to solve a system of partial differential and algebraic equations, when realized there is no implementation of the idea. There is only Mathematica's AsymptoticDSolveValue which is just for ODEs. So I decided to implement it myself. Thanks to the Sympy community we now have some progress. I have one implementation over here by Nicolas CELLIER and some ideas over here by me and some of SymPy developers. I was wondering if we could join forces to come of with a general algorithm, then implementionation on diffrent languages shouldn't be that difficult. I was wondering if you could take a look at this question and help us out. 

FYI I have made a issue on SymPy's github over here which you can follow. 

Best,

Foad