用mathematica求解具有边界条件的连续性方程(偏微分方程),要可执行的...答:bc = {n[t, 0] == 0, D[n[t, x], x] == 0 /. x -> 1};dirac[r_, a_] := Sqrt[a/Pi] Exp[-a r^2]ic = n[0, x] == 2 n0 dirac[x - 1, 100];(*注意这个 2 是必要的,因为你的DiracDelta设在了端点上*)nsol = NDSolveValue[{eq, ic, bc} /. n0 ->...