多变量Scipy-odeint函数

Question

我想解一个二阶微分方程，其中目标函数是两个变量的函数。在

等式如下

del^2 (p) = - exp(p)

where p is a function of R and z.

按照scipy.integrate.odeint文档中的过程，我重新编写了一组四个一阶微分方程。在

q = dp/dR 

s = dp/dz

dq/dR = - exp(p) - ds/dz - q/R

ds/dz = - exp(p) - (q+R dq/dR)/R

使用y = [p, q, s]我可以设置初始条件，比如y0 = [p0, 0, 0]。在

我的问题是如何建立导数函数。在1变量的情况下，你会得到

^{pr2}$

但是，如果一个人用一个双变量的情况下

dydr = [....]

包含关于不同事物的导数。 有什么明显的方法可以解决这个问题，还是我的方法完全错误？在