PDE solving is still limited in sympy. Here is how to solve first order pde
Solve \( u_t(x,t)=u_x(x,t)\)
import sympy as sp x,t = sp.symbols('x t') u = sp.Function('u') pde = sp.Eq( sp.diff(u(x,t),t) , sp.diff(u(x,t),x)) sol = sp.pdsolve(pde) sp.latex(sol)