Wolfram mathematica Can'；t在Mathematica中求解库仑力微分方程

我试图解决一个问题：我必须用Mathematica找到电子在石墨烯晶格中的轨迹。我试着用NDSolve解库仑力方程，并绘制每个方向的结果，但我得到了一个白色的图。有人能帮我吗？先谢谢你。以下是x方向的代码：

coordx = {0.6327, 1.88058, 3.03927, 4.28716, 5.44584, 6.69373, 
   7.85241, 9.10029, 1.9728, 3.22069, 4.37937, 5.62726, 6.78594, 
   8.03382, 9.19251, 10.4404, 3.3129, 4.56079, 5.71947, 6.96736, 
   8.12604, 9.37393, 10.53261, 11.7805, 4.653, 5.90089, 7.05956, 
   8.30746, 9.46614, 10.71403, 11.87271, 13.1206};
me = 9.01*10^-31;
pi = 3.14159;
epsilon0 = 8.854*10^-12;
q = -1.6*10^-19;
Q = 1.6*10^-19;
step = 0.01;
Forzax[p_, r_] := 
 Sum[(Q*q)/(4 pi*epsilon0*Norm[r - p[[i]]]^2), {i, Length[p]}]
Forzax[coordx, {x[t]}];
NDSolve[{x''[t] == Forzax[coordx, {x[t]}]/me , x[0] == 0, 
  x'[0] == 0}, {x[t]}, {t, 0, 1500}]
Show[ParametricPlot[Evaluate[{x[t]} /. %], {t, 0, 1500}, 
  PlotRange -> All]]

---

我不知道你想画什么，但这几处修改似乎可以画出你的函数

sol = NDSolve[{x''[t] == Forzax[coordx, {x[t]}]/me,
    x[0] == 0, x'[0] == 0}, {x}, {t, 0, 1500}];
f = sol[[1, 1, 2]];
Plot[f[t], {t, 0, 1500}, PlotRange -> All]

---

据我所知，您的代码运行正常，但绘图为空，因为您只使用一个函数调用了

ParametricPlot

。从文档中，您可以这样调用

参数化绘图

：

ParametricPlot[{fx[t], fy[t]}, {t, tmin, tmax}]

由于您仍在解决1D问题，并且只有

x[t]

，

ParametericPlot

无法绘制任何图形；它缺少轨迹的

y

坐标。一旦您进行了二维计算，

paramatriclot

应该能够为您提供所需的图形。如果要进行三维计算，应使用

paramatriclot3d

---

一个问题是：你打算如何对石墨烯中的电子进行三维计算？电子在第三维度的运动完全不会遵循牛顿力学，因为它在那个方向受到了很大的限制。事实上，我对在石墨烯中使用牛顿第二定律持谨慎态度，因为石墨烯中的电子表现得像无质量的粒子。我把结果的解释权留给你们，但我身上的物理学家忍不住要加上这句话。

我正试图绘制石墨烯晶格中电子的轨迹。我有我的方程：mx'=Fcx；我的''=Fcy，mz''=Fcz。我必须解这些方程，并在三维x，y，z图中绘制电子的轨迹。实际上我只是想画出一个微分方程的结果（第一步）。我怎样才能解决我的问题？我是数学新手。谢谢你的建议。