Может ли кто-нибудь помочь мне с пакетом R deSolve / ReacTran , который можно использовать для решения любого параболического уравнения в частных производных:
Я нашел похожий пример в книге
Карлайн Соетаерт, Джефф Кэш, Франческа Мацция: Решение дифференциальных уравнений в частных производных в R , глава 9, стр. 179-181,
где решение для уравнения:
Электронная версия книги: Решение дифференциальных уравнений в частных производных в 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 или?)
Заранее благодарю за любую помощь.