R:deSolve-偏微分方程

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

我在书中发现了一个类似的例子 Karline Soetart，Jeff_Cash，Francesca Mazzia:求解R中的偏微分方程，第9章，第179-181页， 其中是方程的解：

本书的在线版本：

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")

但是我不知道如何定义这个方程中的一阶偏导数 对于RUE和RUD的衍生品（更多的是“xgrid”，还是？）

提前感谢您的帮助。

您的PDE实际上是2D+时间（r_e为x，r_d为y）。因此，我认为您应该更深入地研究tran.2D（grid.2D和ode.2D也是如此）。看一看完整的手册cran.r-project.org/web/packages/ReacTran/ReacTran.pdf或cran.r-project.org/web/packages/ReacTran/vignettes/PDE.pdf PS：从未尝试过

r

，因此我不会冒险寻找答案

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