请勿删，评论区有文本
就相当于a=1，b=1，c=a＋bt，d=c＋at然后求e=c＋d

Result = {}; For[t = 0, t <= 10, t = t + 0.1, g = 1，\[CapitalOmega] = 1，n = 0，\[Theta] = 0， \[CapitalDelta] == Sqrt[8*(g^2)*(1 + 2*n) + \[CapitalOmega]^2]， p == (-\[CapitalOmega] + \[CapitalDelta])/2; v == (-\[CapitalOmega] - \[CapitalDelta])/2; Subscript[c, 1] == {(Cos[\[Theta]] + Sin[\[Theta]])*g*Sqrt[ n]*(e^(Ivt_) - e^(Ipt_))}/\[CapitalDelta]，Subscript[c, 2] == {-v*(Cos[\[Theta]] + Sin[\[Theta]])*e^(Ivt_) + p*(Cos[\[Theta]] + Sin[\[Theta]])*e^(Ipt_)}/ 2*\[CapitalDelta] + {(Cos[\[Theta]] - Sin[\[Theta]])/2}* e^(I\[CapitalOmega]t_)， Subscript[c, 3] == {-v*(Cos[\[Theta]] + Sin[\[Theta]])*e^(Ivt_) + p*(Cos[\[Theta]] + Sin[\[Theta]])*e^(Ipt_)}/ 2*\[CapitalDelta] + {(Sin[\[Theta]] - Cos[\[Theta]])/2}* e^(I\[CapitalOmega]t_)， Subscript[c, 4] == {(Cos[\[Theta]] + Sin[\[Theta]])*g*Sqrt[ n + 1]*(e^(Ivt_) - e^(Ipt_))}/\[CapitalDelta]， e == -(Abs[Subscript[c, 2][t_]]^2 + Abs[Subscript[c, 4][ t_]]^2 - \[Sqrt](Abs[Subscript[c, 1][t_]]^4 + 4 Subscript[c, 2][t_] Conjugate[ Subscript[c, 3][t_]] Subscript[c, 3][t_] Conjugate[ Subscript[c, 2][t_]] - 2 Abs[Subscript[c, 1][t_]]^2 Abs[Subscript[c, 4][t_]]^2 + Abs[Subscript[c, 4][t_]]^4)); AppendTo[Result, {t, e}]; ]ListPlot[Result]Export["c:/math/lzlx.txt", Result]

代码问题太多，以至于甚至无法理解出你想要干什么。
总之建议找本书学学吧，吧里有推荐书的贴子。

用nest不会超出递归深度