使用python中的bvp解算器在两个不同域中求解微分方程

Question

我正在尝试使用python中的bvp解算器求解以下耦合微分方程组（变量为p、n和psi）：-

1. d2n/dx2-(d2psi/dx2)n-(dpsi/dx)(dn/dx)=k1 in domain 1
2. d2p/dx2+(d2psi/dx2)p+(dpsi/dx)(dp/dx)=k2 in domain 1
3. d2psi/dx2= k*(p-n) in domain 1
4. d2psi/dx2= 0 in domain 2 where domain 1 is given by x<=L and domain 2 is given by L<x<=(L+t)

根据边界条件：

p(0)= p0, n(0)= n0, psi(0)= 0, p(L)= p0*exp(-psi/Vth), n(L)= n0*exp(psi/Vth), psi(L+t)= 5. 

我可以建立方程。但是考虑到有两个不同的区域需要解方程，我面临一些麻烦。如有任何相关帮助，将不胜感激