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 [email protected]
To unsubscribe from this group, send email to [email protected]
For more options, visit this group at http://groups.google.com/group/sympy?hl=en
-~----------~----~----~----~------~----~------~--~---
