Hey, For the past week I have been trying to figure out how to implement solver for most basic types of PDEs and now it is finally time to show something:
In [1]: from sympy import * In [2]: from sympy.solvers.solvers import * In [3]: from sympy import Derivative as D In [4]: t,x,y,z = symbols('txyz') In [5]: a = Symbol('a', Real=True) In [6]: u = Function('u') In [7]: eq = Eq(D(u(x, t), t) + a*D(u(x, t), x)) In [8]: eq Out[8]: d d a⋅──(u(x, t)) + ──(u(x, t)) = 0 dx dt In [9]: pdesolve(eq, u(x, t)) Out[9]: [x - a⋅t] I have pushed the preliminary current code to my github repository into 'char-method-first-order' branch. Comments, suggestions, patches are more than welcome :) Cheers, Priit :) --~--~---------~--~----~------------~-------~--~----~ You received this message because you are subscribed to the Google Groups "sympy" group. To post to this group, send email to sympy@googlegroups.com To unsubscribe from this group, send email to sympy+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sympy?hl=en -~----------~----~----~----~------~----~------~--~---