diff --git a/README.md b/README.md
index 94d9b9e..23082e0 100644
--- a/README.md
+++ b/README.md
@@ -89,20 +89,34 @@ Dose Estimation by Simulation of China Space Station
#### 俘获辐射带
1. 电子
- 根据文献，在轨道较高（$L>2.5$）的内带中，电子能谱中有$\sim64\%$符合指数分布，$\sim4\%$符合幂律分布。[^6]
+ 在轨道较高（$L>2.5$）的内带中，电子能谱中有$\sim64\%$符合指数分布，$\sim4\%$符合幂律分布。[^6]
- 但在本例中，轨道较低（$L\sim1.06$），因此使用指数分布$j = j_0e^{-\frac{E}{E_0}}$和幂律分布$j=j_0E^{-\alpha}$进行拟合，随后使用`Origin Lab`对两者进行拟合比较，最终选用指数分布模型$j_0=1.333680\times10^6,\ E_0=0.0824$。
+ 但在本例中，轨道较低（$L\sim1.06$），因此使用指数模型$j = j_0e^{-\frac{E}{E_0}}$和幂律模型$j=j_0E^{-\alpha}$进行拟合，随后使用`Origin Lab`对两者进行`BIC test`，最终选用指数模型$j_0=1.333680\times10^6,\ E_0=0.0824$。
- | Model | $R^2$ | BIC | BIC diff |
+ | Model | $R^2$ | BIC | BIC diff |
|:-----:|:-------:|:---------:|:--------:|
| Ex | 0.98107 | 402.50208 | 0 |
| Pow | 0.94279 | 423.51296 | 21.01088 |
+
+
2. 质子
+ 与电子保持一致，额外考虑幂律平移模型$j=j_0(E-\beta)^{-\alpha}$，同电子进行`BIC test`，最终选用幂律平移模型$j_0=559.76377,\ \alpha=-2.40873,\ \beta=-1.53424$。
+
+
+
+ | Model | $R^2$ | BIC | BIC diff |
+ |:-----:|:-------:|:--------:|:--------:|
+ | Ex | 0.9957 | 67.24594 | 0 |
+ | Pow-T | 0.99883 | 39.11969 | 28.12625 |
+
+
+
+
## 空间站结构
1. 尺寸与分区[^7]
diff --git a/debug.txt b/debug.txt
new file mode 100644
index 0000000..8824d44
--- /dev/null
+++ b/debug.txt
@@ -0,0 +1 @@
+22210
diff --git a/docs/ExFit-trapped-electron.png b/docs/ExFit-trapped-electron.png
new file mode 100644
index 0000000..1825a32
Binary files /dev/null and b/docs/ExFit-trapped-electron.png differ
diff --git a/docs/ExFit-trapped-proton.png b/docs/ExFit-trapped-proton.png
new file mode 100644
index 0000000..5c16236
Binary files /dev/null and b/docs/ExFit-trapped-proton.png differ
diff --git a/fit.opju b/fit.opju
index 04c8708..b1b9e33 100644
Binary files a/fit.opju and b/fit.opju differ
