R：deSolve-偏微分方程

Question

有人可以用R软件包 deSolve / ReacTran 帮助我，该软件包可用于求解任何抛物型偏微分方程：

我在书中找到了类似的例子 Karline Soetaert，Jeff_Cash和Francesca Mazzia：用R解决偏微分方程，第9章，第179-181页， 方程的解在哪里：

图书的在线版本：Solving Partial Differential Equations in R

PDE和解决方案的定义如下：

Schrodinger <- function(t, u, parms) {
du <- 1i * tran.1D (C = u, D = 1, dx = xgrid)$dC +
1i * gam * abs(u)ˆ2 * u
list(du) }

N <- 300
xgrid <- setup.grid.1D(-20, 80, N = N)
x <- xgrid$x.mid

out <- ode.1D(y = yini, parms = NULL, func = Schrodinger,
times = times, dimens = 300, method = "adams")

有第二个偏导数的示例：

tran.1D (C = u, D = 1, dx = xgrid)$dC

但是我不确定如何在PDE中定义第一个偏导数 对于 r_e 和 r_d （更多'xgrid'或？）

在此先感谢您的帮助。