我有一个的主方程，其中都是离散的。数值积分是由这个 Matlab程序使用欧拉方法完成的。尽管欧拉方法很粗糙，但数值积分的结果与模拟结果一致，因此程序的正确性不是重点。$P(N_A^+,N_B^+,N_C^+,t)$$N_A^+,N_B^+,N_C^+$

我想看看如何随时间演变。由于除了时间之外还有 3 个坐标，因此很难将变化可视化。理想情况下，我想要类似的东西$P(N_A^+,N_B^+,N_C^+,t)$$P(N_A^+,N_B^+,N_C^+,t)$

但是我在 Matlab 中找不到任何可以做类似事情的函数。我能找到的最接近的是 Matlab 的contourslice函数。

通过运行上面提到的 Matlab 程序

ABC_driver(0.5,100,20,200,10,true);

你会得到类似下面的东西：

你可以猜到会发生什么，但它远没有原子轨道图那么明显。$P(N_A^+,N_B^+,N_C^+,t)$

我的问题是，如何以类似于原子轨道密度图的方式或者任何人都可以提出更好的方法来可视化的时间演变？$P(N_A^+,N_B^+,N_C^+,t)$$P(N_A^+,N_B^+,N_C^+,t)$