我正在尝试解决以下等式：

$$\dfrac{\partial }{\partial \epsilon}\left[B\left(\epsilon\right)f_e\left(1-f_e\right)+D\left(\epsilon\right)\dfrac{\partial f_e}{\partial \epsilon}\right] = R\left(f_e, f_h, E\right)$$

到现在为止，没有什么特别的，我可以解决这个问题（用一个简单的$R$）。

---

但是，这个等式中最重要的部分是一种重组形式$R$：

$$R \approx \dfrac{1}{A-B\cdot f_e}$$

$A,B$- 是已知的常数和$(A-B\cdot f_e)$应该永远是积极的！

但是我使用的数值方案对此一无所知，它只是跳过了这个边界，这意味着没有包含特殊的重组。

的典型值$(A-B\cdot fe) \in [ 10^{-14},10^{-7}]$

有什么方法可以解决吗？主要问题是重组的共振特性$R$.

---

我试图改变$A$在每次迭代中等待收敛，在改变边界参数（这是一个精确改变的参数）之后再次求解。我执行这个程序直到$|A - A_\text{const}| = 1e-4$.

（由于牛顿迭代 ~ 1e-6 中的松弛系数，此方法效果不佳。）