python - sympy を使用して測地線方程式を解く際のエラー

Question

私は一般相対性理論で遊んでいましたが、ほぼ完成しました。しかし、sympy.dsolve を使用して測地線方程式を解くと、常に "NotImplementedError" が報告されます。私のコードは次のようになります。

from gravipy.tensorial import *
from sympy import *

t, r, theta, phi = symbols('t, r, theta, phi')
tau = Symbol('tau')
coord = [s(tau) for s in symbols('t, r, theta, phi', cls=Function)]
chi = Coordinates('chi', [t, r, theta, phi])
M = Symbol('M')
Metric = diag(-(1 - 2 * M / r), 1 / (1 - 2 * M / r), r ** 2,
              r ** 2 * sin(theta) ** 2)
g = MetricTensor('g', chi, Metric)
pprint(g(All, All))
christoffel = Christoffel('gamma', g)
geodesics = Geodesic('geodesics', g, tau)
set_ = []
for i, j in zip(tuple(geodesics(All)), coord):
    pprint(solve(i, Derivative(j, (tau, 2)), dict=True)[0])
    set_.append(Eq(Derivative(j, (tau, 2)), solve(i, Derivative(j, (tau, 2)), dict=True)[0][Derivative(j, (tau, 2))]))

pprint(dsolve(set_))

誰もこの問題を解決する方法を知っていますか?