Я пытаюсь решить следующий набор связанных дифференциальных уравнений (переменные p, n и psi) с помощью решателя bvp в python: -
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.
Я могу составить уравнения. Но, учитывая тот факт, что есть две разные области, в которых уравнения должны быть решены, я столкнулся с некоторыми проблемами. Любая помощь в этом отношении будет принята с благодарностью.