python - ODE のシステム simpy dsolve

Question

複数の ODE を解くにはどうすればよいですか? sympy.dsolve は同じ統合定数を返すため、解決できません。

import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
from sympy.interactive import printing
printing.init_printing(use_latex=True)
from sympy import *

x,l,F,k_0,q_0,n,a,C1,C2,C3,C4,Φ1,Φ2,D1,D2=symbols('x l F k_0 q_0 n a C1 C2 C3 C4 Φ_1 Φ_2 D1 D2')
n=5

q=0
k1=k_0
k2=n*k_0

eq1=Function('eq1')(x)
eq2=Function('eq2')(x)

w1=dsolve(k1.diff(x)*eq1.diff(x)+k1*eq1.diff(x,2)+q).rhs
w2=dsolve(k2.diff(x)*eq2.diff(x)+k2*eq2.diff(x,2)+q).rhs

display(w1)
display(w2)

コードは次を返します。

1+2*

1+2*

私は次のようなものが欲しいのですが：

1+2*

3+4*