level 1
无忧并无期
楼主
源代码如下：
Clear["Global`*"];
Clear[Derivative];
s = NDSolve[{D[rho[x, t], t] + D[rho[x, t] V[x, t], x] == 0,
D[V[x, t], t] + D[(1/2 V[x, t]^2 + rho[x, t]), x] ==
10 - rho[x, t] - V[x, t], rho[0, t] == rho[20, t],
rho[x, 0] == Sin[x]}, {rho, V}, {{x, 0, 20}, {t, 0, 20}}];
请问出了什么错，该如何解决，谢谢各位了
（顺便，这种周期性边界条件的pde用mma该怎呢解啊）
2021年06月19日 08点06分
1
Clear["Global`*"];
Clear[Derivative];
s = NDSolve[{D[rho[x, t], t] + D[rho[x, t] V[x, t], x] == 0,
D[V[x, t], t] + D[(1/2 V[x, t]^2 + rho[x, t]), x] ==
10 - rho[x, t] - V[x, t], rho[0, t] == rho[20, t],
rho[x, 0] == Sin[x]}, {rho, V}, {{x, 0, 20}, {t, 0, 20}}];
请问出了什么错，该如何解决，谢谢各位了
（顺便，这种周期性边界条件的pde用mma该怎呢解啊）
