Először meghívjuk a sympy csomagot.
from sympy import *
init_printing()
Ezután megadjuk az x változót, mint szimbólumot, az f-t, pedig mint függvényt. Ez utóbbit a function() függvény segítségével tudjuk megtenni.
f = Function("f")
x = symbols("x")
Majd megoldjuk a diffegyenletet. Közönséges diffegyenletnek a megoldását legegyszerűbb a dsolve() függvénnyel lehet, melynek működése hasonló a solve()-éhoz, csak a megoldása függvény:
dsolve(diff(f(x), x, 2) + 9 * f(x) - 1, f(x))