我是一名物理专业的学生，​​在编程方面几乎没有经验。这个学期我们应该做一个热力学的计算项目。我必须解决这两个耦合的 diff eqns：

$$\begin{aligned} ω (p,T ) &= \frac{p ^2}{2m} + \frac{2}{N}\sum_q f(p − q)\,n (q) − \frac{1}{N^2}\sum_{s,t} f(s − t)\,n(s)\,n(t) \quad\text{and} \\ n(p) &= \frac{1}{\exp\left[\cfrac{ω (p,T ) − \mu}{kT}\right] − 1} \end{aligned}$$

$\omega$是每个玻色子的能量。

$p$是玻色子的动量。

$n(p)$是具有动量的状态下的玻色子数。$p$

$f$是

$$f(p)=\frac{1}{2} \left[\epsilon_0- \frac{p^2}{2m}\right]$$

$\epsilon_0$是 0 K 的基本激发能。

$N$是玻色子的总数。

$T$是温度，是常数。$k$

有人可以指导我使用任何简单的方法来为这个问题生成一些粗略的解决方案吗？基于一篇论文：“基于非线性理论的 0 - 2.1 K 之间超流氦的比热评估”作者：Shosuke Sasaki (arXiv:0807.1361v1 [cond-mat.other] 2008 年 7 月 9 日