def dX_dt(X, t):
    vals = dict(x=X[0], y=X[1], t=t)
return [eq.evalf(subs=vals) for eq in systemOfEquations]
SymPy/SciPy: решение системы обыкновенных дифференциальных уравнений с разными переменными