当我阅读本文http://articles.adsabs.harvard.edu/cgi-bin/nph-iarticle_query?1976ApJ...209..214B&data_type=PDF_HIGH&whole_paper=YES&type=PRINTER&filetype=.pdf时 我尝试数值求解eq（49），这似乎是一个fokker-planck方程，我发现有限差分法不起作用，它很不稳定。 有人知道如何解决吗？
答案 0 :(得分：0)
computational science stack exchange是您可以提出并希望得到答案的地方。或者，您可以尝试使用其physics表弟。您引用的方程是整数微分方程，是相当非线性的...福克-普朗克方程。绝对不是典型的福克木板。
您可以尝试使用有限差分或有限元素离散化函数g(x,t)的空格部分。毕竟，0 < x < x_max有边界条件。您还必须离散化相应的集成。那么也许有限元可能更合适？有限元素意味着您可以将g（x，t）编写为一系列精心选择的基础，这些基础经过紧凑地支持的足够简单的简单函数Bj（x）：j = 1 ... N，范围为[0, x_max]
g(x,t) = sum_j=1:N gj(t)*Bj(x)
这会将您的函数变成一个（大）向量gj(t) = g(x_j, t), for j = 1, 1, ...., N。结果，您将获得一个非线性的ODE系统
dgj(t)/dt = Qj(g1(t), g2(t), ..., gN(t))
j = 1 ... N
之后，使用Runge-Kutta之类的东西对ODE系统进行数值积分。