diff --git a/docs/source/Tutorial/dynamic/integ_converg.rst b/docs/source/Tutorial/dynamic/integ_converg.rst
deleted file mode 100644
index 14efb24d..00000000
--- a/docs/source/Tutorial/dynamic/integ_converg.rst
+++ /dev/null
@@ -1,13 +0,0 @@
-积分收敛性
-===================
-
-:Author: Zhu Dengda
-:Email:  zhudengda@mail.iggcas.ac.cn
-
------------------------------------------------------------
-
-通过输出积分过程文件，直观感受源点和场点深度接近时，积分收敛性的变化，以及峰谷平均法的作用。
-
-动态和静态情况类似，这里以动态解为例。
-
-TODO
\ No newline at end of file
diff --git a/docs/source/Tutorial/index.rst b/docs/source/Tutorial/index.rst
index e8e1fd2e..bd23bb88 100644
--- a/docs/source/Tutorial/index.rst
+++ b/docs/source/Tutorial/index.rst
@@ -82,7 +82,7 @@
    static/static_syn
    dynamic/strain_stress
    static/static_strain_stress
-   dynamic/integ_converg
+   integ_converg/integ_converg
 
 
 
diff --git a/docs/source/Tutorial/integ_converg/integ_converg.rst b/docs/source/Tutorial/integ_converg/integ_converg.rst
new file mode 100644
index 00000000..e3608961
--- /dev/null
+++ b/docs/source/Tutorial/integ_converg/integ_converg.rst
@@ -0,0 +1,205 @@
+积分收敛性
+===================
+
+:Author: Zhu Dengda
+:Email:  zhudengda@mail.iggcas.ac.cn
+
+-----------------------------------------------------------
+
+通过输出核函数文件，观察源点和场点深度接近时，积分收敛性的变化，以及峰谷平均法的作用。 **具体积分表达式以及分类详见** :ref:`gfunc_rst`。
+
+核函数文件
+---------------
+该文件记录了波数积分过程中不同波数对应的核函数值。以下用法展示了如何输出该文件。
+
+.. tabs::  
+
+    .. tab:: C 
+
+        .. literalinclude:: run/run.sh
+            :language: bash
+            :start-after: BEGIN DGRN
+            :end-before: END DGRN
+
+        输出的核函数文件路径会在  :rst:dir:`GRN_grtstats/milrow_{depsrc}_{deprcv}/` 路径下。
+
+    .. tab:: Python
+
+        .. literalinclude:: run/run.py
+            :language: python
+            :start-after: BEGIN DGRN
+            :end-before: END DGRN
+
+        输出的核函数文件路径会在自定义路径下。
+
+
+C和Python导出的核函数文件是一致的，底层调用的是相同的函数。文件名称格式为 ``K_{iw}_{freq}``，其中 ``{iw}`` 表示频率索引值， ``{freq}`` 表示对应频率(Hz)。文件为自定义的二进制文件， **强烈建议使用Python进行读取及后续处理**。这里还是给出两种读取方法。
+
+.. tabs::  
+
+    .. tab:: C 
+
+        :command:`grt.k2a` 程序可将单个核函数文件转为文本格式。
+
+        .. literalinclude:: run/run.sh
+            :language: bash
+            :start-after: BEGIN grt.k2a
+            :end-before: END grt.k2a
+
+        输出的文件如下，
+
+        .. literalinclude:: run/stats_head
+            :language: text
+
+        后续你可以选择习惯的方式读取和处理。
+
+    .. tab:: Python
+
+        .. literalinclude:: run/run.py
+            :language: python
+            :start-after: BEGIN read statsfile
+            :end-before: END read statsfile
+
+其中除了波数 ``k`` 外，每条结果的命名格式均为 ``{srcType}_{q/w/v}{m}``，与 :ref:`gfunc_rst` 部分介绍的积分公式中的核函数 :math:`q_m, w_m, v_m` 保持一致。
+
+
+可视化
+-------------
+以下将使用Python进行图件绘制。 **在Python函数中指定震源类型、阶数、积分类型，可自动绘制核函数、被积函数和积分值随波数的变化**，其中积分类型对应 :ref:`gfunc_rst` 部分介绍的4种类型。
+
+.. literalinclude:: run/run.py
+    :language: python
+    :start-after: BEGIN plot stats
+    :end-before: END plot stats
+
+.. image:: run/DC_20.png
+    :align: center
+
+
+可以通过指定 ``RorI=False`` 参数来指定绘制虚部，传入 ``RorI=2`` 表示实虚部都绘制。
+
+.. literalinclude:: run/run.py
+    :language: python
+    :start-after: BEGIN plot stats RI
+    :end-before: END plot stats RI
+
+.. image:: run/DC_20_RI.png
+    :align: center
+
+
+从图中可以看到，当波数增加到一定范围以上，积分收敛良好，无振荡现象。
+
+
+
+当震源和场点深度接近/相等时，积分收敛速度变慢
+----------------------------------------------
+
+假设其它参数不变，我们调整 **震源深度为0.1km**，再计算格林函数，读取核函数文件，绘制图像。
+
+.. tabs::  
+
+    .. tab:: C 
+
+        .. literalinclude:: run/run.sh
+            :language: bash
+            :start-after: BEGIN DEPSRC 0.0 DGRN
+            :end-before: END DEPSRC 0.0 DGRN
+
+    .. tab:: Python
+
+        .. literalinclude:: run/run.py
+            :language: python
+            :start-after: BEGIN DEPSRC 0.0 DGRN
+            :end-before: END DEPSRC 0.0 DGRN
+
+
+.. image:: run/DC_20_0.0_RI.png
+    :align: center
+
+从图中可以清晰地看到，相比震源深度2km时，积分收敛速度明显变慢，积分值振荡，这要求增加波数积分上限，但这必然降低计算效率。
+
+你可以尝试 **当震源和场点位于同一深度**，收敛振荡现象也很明显。要注意的是这和是否在地表无关，即使震源和场点都在地下，深度接近时积分收敛也会变慢。
+
+
+峰谷平均法
+-------------------------
+具体原理很简洁易懂，从以下图像中你也能了解大概。详见 :ref:`(Zhang et al., 2003) <zhang_2003>` :ref:`(张海明, 2021) <zhang_book_2021>`。
+
+在以上示例中，当震源和场点深度接近时，你会注意到输出的核函数文件增加了，这是因为程序内部设定 **当震源和场点深度差小于1km时，自动使用峰谷平均法（PTAM）**。具体流程为，程序中在进行完离散波数积分后，继续增加k（不同震中距使用不同dk），使用PTAM寻找足够数量的波峰和波谷（内部设定为36个），再对这些波峰波谷取缩减序列 :math:`M_i \leftarrow 0.5\times(M_i + M_{i+1})` ，得到估计的积分收敛值。取缩减序列的C代码如下，
+
+.. code-block:: C 
+
+    for(int n=n1; n>1; --n){
+        for(int i=0; i<n-1; ++i){
+            arr[i] = 0.5*(arr[i] + arr[i+1]);
+        }
+    }
+
+
+在 ``K_{iw}_{freq}`` 文件同级目录下，程序把 **PTAM过程中的核函数以及积分峰谷位置分为两个文件** 保存在 ``PTAM_{ir}_{dist}/`` 目录下（ ``{ir}`` 为震中距索引， ``{dist}`` 为震中距），其中  ``PTAM_{ir}_{dist}/K_{iw}_{freq}`` 为核函数文件（格式不变）， ``PTAM_{ir}_{dist}/PTAM_{iw}_{freq}`` 为峰谷文件，其中记录积分值的峰谷。
+
+.. note:: 
+
+    :command:`grt.k2a` 程序也支持将 ``PTAM_{ir}_{dist}/PTAM_{iw}_{freq}`` 文件转为文本格式，
+
+    .. literalinclude:: run/run.sh
+        :language: bash
+        :start-after: BEGIN grt.k2a ptam
+        :end-before: END grt.k2a ptam
+
+    输出的文件如下，
+
+        .. literalinclude:: run/ptam_stats_head
+            :language: text
+
+    记录了不同震源、不同积分类型的峰谷位置。
+
+
+故为了绘制完整的波数积分+PTAM过程，涉及三个文件。不过Python函数已经做了优化，由于程序输出的目录结构固定，文件命名格式固定，只需传入 ``PTAM_{ir}_{dist}/PTAM_{iw}_{freq}`` 文件，Python函数会自动读取对应的三个文件。
+
+.. literalinclude:: run/run.py
+    :language: python
+    :start-after: BEGIN plot ptam
+    :end-before: END plot ptam
+
+.. image:: run/DC_20_0.0_ptam_RI.png
+    :align: center
+
+
+红十字为选取的波峰波谷，再取缩减序列即可得到估计的积分收敛值。
+
+
+
+静态解的积分收敛性
+-------------------------
+以上部分是以动态解为例，静态解从积分类型、收敛特征、文件格式、绘图完全类似，只是不再有频率索引值。
+
+假设震源深度0.1km，场点位于地表，场点仅定义一个点(2,2)做示例，这里直接给出脚本。
+
+.. tabs::  
+
+    .. tab:: C 
+
+        .. literalinclude:: run/run.sh
+            :language: bash
+            :start-after: BEGIN SGRN
+            :end-before: END SGRN
+
+    .. tab:: Python
+
+        .. literalinclude:: run/run.py
+            :language: python
+            :start-after: BEGIN SGRN
+            :end-before: END SGRN
+
+
++ 只使用离散波数积分
+
+.. image:: run/DC_20_0.1_static.png
+    :align: center
+
++ 使用峰谷平均法
+  
+.. image:: run/DC_20_0.1_ptam_static.png
+    :align: center
+
diff --git a/docs/source/Tutorial/integ_converg/run/DC_20.png b/docs/source/Tutorial/integ_converg/run/DC_20.png
new file mode 100644
index 00000000..2a0fc806
Binary files /dev/null and b/docs/source/Tutorial/integ_converg/run/DC_20.png differ
diff --git a/docs/source/Tutorial/integ_converg/run/DC_20_0.0_RI.png b/docs/source/Tutorial/integ_converg/run/DC_20_0.0_RI.png
new file mode 100644
index 00000000..33f49b66
Binary files /dev/null and b/docs/source/Tutorial/integ_converg/run/DC_20_0.0_RI.png differ
diff --git a/docs/source/Tutorial/integ_converg/run/DC_20_0.0_ptam_RI.png b/docs/source/Tutorial/integ_converg/run/DC_20_0.0_ptam_RI.png
new file mode 100644
index 00000000..cb7dd9b3
Binary files /dev/null and b/docs/source/Tutorial/integ_converg/run/DC_20_0.0_ptam_RI.png differ
diff --git a/docs/source/Tutorial/integ_converg/run/DC_20_0.1_ptam_static.png b/docs/source/Tutorial/integ_converg/run/DC_20_0.1_ptam_static.png
new file mode 100644
index 00000000..049fa6f6
Binary files /dev/null and b/docs/source/Tutorial/integ_converg/run/DC_20_0.1_ptam_static.png differ
diff --git a/docs/source/Tutorial/integ_converg/run/DC_20_0.1_static.png b/docs/source/Tutorial/integ_converg/run/DC_20_0.1_static.png
new file mode 100644
index 00000000..aad68103
Binary files /dev/null and b/docs/source/Tutorial/integ_converg/run/DC_20_0.1_static.png differ
diff --git a/docs/source/Tutorial/integ_converg/run/DC_20_RI.png b/docs/source/Tutorial/integ_converg/run/DC_20_RI.png
new file mode 100644
index 00000000..a006b9ac
Binary files /dev/null and b/docs/source/Tutorial/integ_converg/run/DC_20_RI.png differ
diff --git a/docs/source/Tutorial/integ_converg/run/milrow b/docs/source/Tutorial/integ_converg/run/milrow
new file mode 100644
index 00000000..a99f4912
--- /dev/null
+++ b/docs/source/Tutorial/integ_converg/run/milrow
@@ -0,0 +1,9 @@
+0.2  3.4  1.7  2.3  9e10  9e10
+0.6  3.7  1.9  2.4  9e10  9e10
+0.5  4.2  2.1  2.4  9e10  9e10
+0.5  4.6  2.3  2.5  9e10  9e10
+0.7  4.9  2.8  2.6  9e10  9e10
+0.5  5.1  2.9  2.7  9e10  9e10
+6.0  5.9  3.3  2.7  9e10  9e10
+1.   6.9  4.0  2.8  9e10  9e10
+0.0  8.2  4.7  3.2  9e10  9e10
\ No newline at end of file
diff --git a/docs/source/Tutorial/integ_converg/run/ptam_stats_head b/docs/source/Tutorial/integ_converg/run/ptam_stats_head
new file mode 100644
index 00000000..eb3e342c
--- /dev/null
+++ b/docs/source/Tutorial/integ_converg/run/ptam_stats_head
@@ -0,0 +1,11 @@
+#     sum_EXP_00_k                        sum_EXP_00       sum_VF_00_k                         sum_VF_00       sum_DC_00_k                         sum_DC_00      sum_EXP_02_k                        sum_EXP_02       sum_VF_02_k                         sum_VF_02       sum_DC_02_k                         sum_DC_02       sum_HF_10_k                         sum_HF_10       sum_DC_10_k                         sum_DC_10       sum_HF_11_k                         sum_HF_11       sum_DC_11_k                         sum_DC_11       sum_HF_12_k                         sum_HF_12       sum_DC_12_k                         sum_DC_12       sum_HF_13_k                         sum_HF_13       sum_DC_13_k                         sum_DC_13       sum_DC_20_k                         sum_DC_20       sum_DC_21_k                         sum_DC_21       sum_DC_22_k                         sum_DC_22       sum_DC_23_k                         sum_DC_23
+    2.67622465e+01    6.66555330e+02-6.52355489e+03J    2.67622568e+01   -6.71825014e+01+1.37589835e+03J    2.67622465e+01   -2.66622132e+03+2.60942196e+04J    2.66053630e+01    1.35857785e+04+7.55803841e+02J    2.66053633e+01   -1.97736289e+03-2.98716114e+01J    2.66053630e+01   -5.43431139e+04-3.02321536e+03J    2.66053604e+01   -6.25363762e+02-3.80092044e+02J    2.65219514e+01    2.42729886e-11-5.31849265e-11J    2.67622542e+01   -1.03466177e+01+9.34738740e+00J    2.65154089e+01    6.92762231e-13-1.64735342e-13J    2.67622568e+01    6.71825014e+01-1.37589835e+03J    2.66205772e+01    1.83164486e-10+3.18094236e-12J    2.66053605e+01    1.38503321e+03+8.64167945e+02J    2.65098060e+01    0.00000000e+00+0.00000000e+00J    2.67622465e+01    1.33311066e+03-1.30471098e+04J    2.66045993e+01   -4.06156039e+02-2.45114486e+02J    2.66046025e+01    2.71917459e+04+1.23617601e+03J    2.67622465e+01   -1.80515555e+04+2.77535386e+04J
+    2.70764072e+01    2.52519538e+03-6.56643813e+03J    2.70764171e+01   -1.32379593e+02+1.37827970e+03J    2.70764072e+01   -1.01007815e+04+2.62657525e+04J    2.69195212e+01    1.26971000e+04+7.88692689e+02J    2.69195216e+01   -1.79715336e+03-3.51777813e+01J    2.69195212e+01   -5.07883999e+04-3.15477075e+03J    2.69195189e+01   -4.86255714e+02-3.83331419e+02J    2.66699769e+01    2.62038907e-11-5.32231922e-11J    2.70764150e+01   -1.00454716e+01+9.34041115e+00J    2.66063826e+01    6.77250543e-13-1.64858007e-13J    2.70764171e+01    1.32379593e+02-1.37827970e+03J    2.66898801e+01    1.86562483e-10+3.02839093e-12J    2.69195190e+01    1.60593887e+03+8.59028414e+02J    2.65883458e+01    0.00000000e+00+0.00000000e+00J    2.70764072e+01    5.05039077e+03-1.31328763e+04J    2.69187667e+01   -4.22517373e+02-2.44734408e+02J    2.69187695e+01    2.54141606e+04+1.30196624e+03J    2.70764073e+01   -1.21489937e+04+2.76171490e+04J
+    2.73905679e+01    6.73937049e+02-6.52446155e+03J    2.73905774e+01   -6.93843871e+01+1.37603618e+03J    2.73905679e+01   -2.69574820e+03+2.60978462e+04J    2.72336795e+01    1.35651038e+04+7.57392855e+02J    2.72336799e+01   -1.97297785e+03-3.01233966e+01J    2.72336795e+01   -5.42604152e+04-3.02957142e+03J    2.72336774e+01   -6.23162563e+02-3.80200555e+02J    2.67559430e+01    3.41902954e-11-5.32000344e-11J    2.73905757e+01   -1.03383271e+01+9.34715594e+00J    2.66773218e+01    6.90513962e-13-1.65285048e-13J    2.73905774e+01    6.93843871e+01-1.37603618e+03J    2.69233162e+01    2.01682508e-10+3.07569034e-12J    2.72336775e+01    1.38857526e+03+8.64017751e+02J    2.66668856e+01    0.00000000e+00+0.00000000e+00J    2.73905679e+01    1.34787410e+03-1.30489231e+04J    2.72329339e+01   -4.06424282e+02-2.45106147e+02J    2.72329362e+01    2.71503815e+04+1.23935525e+03J    2.73905680e+01   -1.80264891e+04+2.77511940e+04J
+    2.77047285e+01    2.51888998e+03-6.56562927e+03J    2.77047378e+01   -1.30349223e+02+1.37815606e+03J    2.77047285e+01   -1.00755599e+04+2.62625171e+04J    2.75478378e+01    1.27159979e+04+7.87267798e+02J    2.75478383e+01   -1.80122353e+03-3.49512546e+01J    2.75478378e+01   -5.08639916e+04-3.14907119e+03J    2.75478359e+01   -4.88323031e+02-3.83233078e+02J    2.69858115e+01    1.64622008e-11-5.31797570e-11J    2.77047363e+01   -1.00533787e+01+9.34063305e+00J    2.69057304e+01    7.38538662e-13-1.65293921e-13J    2.77047378e+01    1.30349223e+02-1.37815606e+03J    2.70777608e+01    1.97480539e-10+3.19880445e-12J    2.75478361e+01    1.60258550e+03+8.59167341e+02J    2.67454255e+01    0.00000000e+00+0.00000000e+00J    2.77047285e+01    5.03777996e+03-1.31312585e+04J    2.75471010e+01   -4.22260093e+02-2.44742526e+02J    2.75471029e+01    2.54519700e+04+1.29911541e+03J    2.77047287e+01   -1.21712548e+04+2.76192913e+04J
+    2.80188892e+01    6.79285408e+02-6.52518567e+03J    2.80188982e+01   -7.12624018e+01+1.37614758e+03J    2.80188892e+01   -2.71714163e+03+2.61007427e+04J    2.78619961e+01    1.35477811e+04+7.58676203e+02J    2.78619967e+01   -1.96918808e+03-3.03281499e+01J    2.78619961e+01   -5.41911242e+04-3.03470481e+03J    2.78619945e+01   -6.21215273e+02-3.80290105e+02J    2.71012124e+01    3.30180096e-11-5.30607034e-11J    2.80188970e+01   -1.03307817e+01+9.34694370e+00J    2.69926585e+01    6.94850926e-13-1.64914019e-13J    2.80188982e+01    7.12624018e+01-1.37614758e+03J    2.71405501e+01    1.98200699e-10+3.14969160e-12J    2.78619946e+01    1.39175621e+03+8.63888863e+02J    2.68239653e+01    0.00000000e+00+0.00000000e+00J    2.80188892e+01    1.35857082e+03-1.30503713e+04J    2.78612679e+01   -4.06671081e+02-2.45098278e+02J    2.78612694e+01    2.71157239e+04+1.24192286e+03J    2.80188894e+01   -1.80067770e+04+2.77492322e+04J
+    2.83330498e+01    2.51439617e+03-6.56497902e+03J    2.83330585e+01   -1.28607138e+02+1.37805527e+03J    2.83330498e+01   -1.00575847e+04+2.62599161e+04J    2.81761545e+01    1.27319178e+04+7.86107252e+02J    2.81761552e+01   -1.80476247e+03-3.47653851e+01J    2.81761545e+01   -5.09276711e+04-3.14442901e+03J    2.81761531e+01   -4.90162208e+02-3.83151174e+02J    2.71755016e+01    2.84552368e-11-5.29861787e-11J    2.83330576e+01   -1.00605833e+01+9.34083574e+00J    2.70775518e+01    7.18655762e-13-1.64877890e-13J    2.83330585e+01    1.28607138e+02-1.37805527e+03J    2.72057560e+01    1.95856080e-10+3.27609856e-12J    2.81761532e+01    1.59956258e+03+8.59287249e+02J    2.69025051e+01    0.00000000e+00+0.00000000e+00J    2.83330498e+01    5.02879234e+03-1.31299580e+04J    2.81754346e+01   -4.22023282e+02-2.44750129e+02J    2.81754358e+01    2.54838206e+04+1.29679353e+03J    2.83330500e+01   -1.21886443e+04+2.76210913e+04J
+    2.86472104e+01    6.83013242e+02-6.52577117e+03J    2.86472189e+01   -7.28827184e+01+1.37623912e+03J    2.86472104e+01   -2.73205297e+03+2.61030847e+04J    2.84903130e+01    1.35331158e+04+7.59729619e+02J    2.84903137e+01   -1.96587449e+03-3.04975494e+01J    2.84903130e+01   -5.41324632e+04-3.03891848e+03J    2.84903117e+01   -6.19473898e+02-3.80365318e+02J    2.73193464e+01    3.17829035e-11-5.29670298e-11J    2.86472181e+01   -1.03238981e+01+9.34675032e+00J    2.71675681e+01    7.25092543e-13-1.64995602e-13J    2.86472189e+01    7.28827184e+01-1.37623912e+03J    2.73247322e+01    2.02381773e-10+3.32872257e-12J    2.84903119e+01    1.39463388e+03+8.63777020e+02J    2.69810449e+01    0.00000000e+00+0.00000000e+00J    2.86472104e+01    1.36602648e+03-1.30515423e+04J    2.84896011e+01   -4.06898398e+02-2.45090950e+02J    2.84896021e+01    2.70863837e+04+1.24403039e+03J    2.86472106e+01   -1.79915090e+04+2.77475776e+04J
+    2.89613710e+01    2.51135736e+03-6.56445052e+03J    2.89613792e+01   -1.27096316e+02+1.37797185e+03J    2.89613710e+01   -1.00454294e+04+2.62578021e+04J    2.88044714e+01    1.27454566e+04+7.85147773e+02J    2.88044722e+01   -1.80787280e+03-3.46104256e+01J    2.88044714e+01   -5.09818265e+04-3.14059109e+03J    2.88044703e+01   -4.91814753e+02-3.83081845e+02J    2.75044851e+01    4.05272699e-11-5.30701150e-11J    2.89613787e+01   -1.00671646e+01+9.34102012e+00J    2.73243991e+01    7.48357530e-13-1.64621873e-13J    2.89613792e+01    1.27096316e+02-1.37797185e+03J    2.75113883e+01    1.92670301e-10+3.35712528e-12J    2.88044705e+01    1.59681879e+03+8.59391825e+02J    2.70595847e+01    0.00000000e+00+0.00000000e+00J    2.89613710e+01    5.02271471e+03-1.31289010e+04J    2.88037675e+01   -4.21804977e+02-2.44757180e+02J    2.88037682e+01    2.55109071e+04+1.29487395e+03J    2.89613713e+01   -1.22019695e+04+2.76226148e+04J
+    2.92755316e+01    6.85430224e+02-6.52624925e+03J    2.92755395e+01   -7.42947329e+01+1.37631540e+03J    2.92755316e+01   -2.74172089e+03+2.61049970e+04J    2.91186299e+01    1.35205919e+04+7.60606336e+02J    2.91186308e+01   -1.96294815e+03-3.06397871e+01J    2.91186299e+01   -5.40823677e+04-3.04242534e+03J    2.91186290e+01   -6.17902335e+02-3.80429446e+02J    2.75540770e+01    4.00135248e-11-5.30596745e-11J    2.92755392e+01   -1.03176014e+01+9.34657458e+00J    2.74979689e+01    7.12410051e-13-1.64503655e-13J    2.92755395e+01    7.42947329e+01-1.37631540e+03J    2.75966853e+01    1.96164020e-10+3.36759361e-12J    2.91186292e+01    1.39725398e+03+8.63679015e+02J    2.71381245e+01    0.00000000e+00+0.00000000e+00J    2.92755316e+01    1.37086045e+03-1.30524985e+04J    2.91179337e+01   -4.07108158e+02-2.45084173e+02J    2.91179342e+01    2.70613279e+04+1.24578436e+03J    2.92755318e+01   -1.79799669e+04+2.77461729e+04J
+...
diff --git a/docs/source/Tutorial/integ_converg/run/run.py b/docs/source/Tutorial/integ_converg/run/run.py
new file mode 100644
index 00000000..a518ce31
--- /dev/null
+++ b/docs/source/Tutorial/integ_converg/run/run.py
@@ -0,0 +1,116 @@
+# BEGIN DGRN
+import numpy as np
+import pygrt 
+
+modarr = np.loadtxt("milrow")
+
+depsrc = 2.0
+deprcv = 0.0
+
+pymod = pygrt.PyModel1D(modarr, depsrc=depsrc, deprcv=deprcv)
+
+# 通过statsfile参数自定义核函数文件的输出目录, statsidxs指定想输出的频率索引值
+distarr = [5,8,10]
+stgrnLst = pymod.compute_grn(
+    distarr=distarr, nt=500, dt=0.02,
+    statsfile=f"pygrtstats_{depsrc}_{deprcv}", statsidxs=[50,100]
+)
+# END DGRN
+
+
+
+# BEGIN read statsfile
+# 可使用通配符简化输入，因为对应索引值下只会有一个文件
+# 返回的是自定义类型的numpy数组
+statsdata = pygrt.utils.read_statsfile(f"pygrtstats_{depsrc}_{deprcv}/K_0050_*")
+print(statsdata.dtype)
+# [('k', '<f8'), ('EXP_q0', '<c16'), ('EXP_w0', '<c16'), ('VF_q0', '<c16'), ('VF_w0', '<c16'), ('HF_q1', '<c16'), ('HF_w1', '<c16'), ('HF_v1', '<c16'), ('DC_q0', '<c16'), ('DC_w0', '<c16'), ('DC_q1', '<c16'), ('DC_w1', '<c16'), ('DC_v1', '<c16'), ('DC_q2', '<c16'), ('DC_w2', '<c16'), ('DC_v2', '<c16')]
+# END read statsfile
+
+
+# BEGIN plot stats
+ir = 2
+dist=distarr[ir]
+srctype="DC"
+mtype="2"
+ptype="0"
+fig, ax = pygrt.utils.plot_statsdata(statsdata, dist=dist, srctype=srctype, mtype=mtype, ptype=ptype)
+fig.savefig(f"{srctype}_{mtype}{ptype}.png", dpi=100)
+# END plot stats
+
+# BEGIN plot stats RI
+ir = 2
+dist=distarr[ir]
+srctype="DC"
+mtype="2"
+ptype="0"
+fig, ax = pygrt.utils.plot_statsdata(statsdata, dist=dist, srctype=srctype, mtype=mtype, ptype=ptype, RorI=2)
+fig.savefig(f"{srctype}_{mtype}{ptype}_RI.png", dpi=100)
+# END plot stats RI
+
+
+
+# BEGIN DEPSRC 0.0 DGRN
+depsrc = 0.0
+deprcv = 0.0
+pymod = pygrt.PyModel1D(modarr, depsrc=depsrc, deprcv=deprcv)
+
+stgrnLst = pymod.compute_grn(
+    distarr=distarr, nt=500, dt=0.02,
+    statsfile=f"pygrtstats_{depsrc}_{deprcv}", statsidxs=[50,100]
+)
+
+statsdata = pygrt.utils.read_statsfile(f"pygrtstats_{depsrc}_{deprcv}/K_0050_*")
+
+dist=10
+srctype="DC"
+mtype="2"
+ptype="0"
+fig, ax = pygrt.utils.plot_statsdata(statsdata, dist=dist, srctype=srctype, mtype=mtype, ptype=ptype, RorI=2)
+fig.savefig(f"{srctype}_{mtype}{ptype}_{depsrc}_RI.png", dpi=100)
+# END DEPSRC 0.0 DGRN
+
+
+# BEGIN plot ptam
+ir = 2
+statsdata1, statsdata2, ptamdata = pygrt.utils.read_statsfile_ptam(f"pygrtstats_{depsrc}_{deprcv}/PTAM_{ir:04d}_*/PTAM_0050_*")
+
+srctype="DC"
+mtype="2"
+ptype="0"
+fig, ax = pygrt.utils.plot_statsdata_ptam(statsdata1, statsdata2, ptamdata, dist=distarr[ir], srctype=srctype, mtype=mtype, ptype=ptype, RorI=2)
+fig.savefig(f"{srctype}_{mtype}{ptype}_{depsrc}_ptam_RI.png", dpi=100)
+# END plot ptam
+
+
+
+# BEGIN SGRN
+import numpy as np
+import pygrt 
+
+modarr = np.loadtxt("milrow")
+
+depsrc = 0.1
+deprcv = 0.0
+pymod = pygrt.PyModel1D(modarr, depsrc=depsrc, deprcv=deprcv)
+
+xarr = np.array([2.0])
+yarr = np.array([2.0])
+static_grn = pymod.compute_static_grn(xarr, yarr, statsfile=f"static_pygrtstats_{depsrc}_{deprcv}")
+
+ir = 0
+distarr = np.sqrt(xarr**2 + yarr**2)
+statsdata1, statsdata2, ptamdata = pygrt.utils.read_statsfile_ptam(f"static_pygrtstats_{depsrc}_{deprcv}/PTAM_{ir:04d}_*/PTAM")
+
+srctype="DC"
+mtype="2"
+ptype="0"
+# 只使用离散波数积分的积分变化
+fig, ax = pygrt.utils.plot_statsdata(statsdata1, dist=distarr[ir], srctype=srctype, mtype=mtype, ptype=ptype, RorI=True)
+fig.savefig(f"{srctype}_{mtype}{ptype}_{depsrc}_static.png", dpi=100)
+
+# 使用了峰谷平均法的积分变化
+fig, ax = pygrt.utils.plot_statsdata_ptam(statsdata1, statsdata2, ptamdata, dist=distarr[ir], srctype=srctype, mtype=mtype, ptype=ptype, RorI=True)
+fig.savefig(f"{srctype}_{mtype}{ptype}_{depsrc}_ptam_static.png", dpi=100)
+
+# END SGRN
diff --git a/docs/source/Tutorial/integ_converg/run/run.sh b/docs/source/Tutorial/integ_converg/run/run.sh
new file mode 100755
index 00000000..bc02c574
--- /dev/null
+++ b/docs/source/Tutorial/integ_converg/run/run.sh
@@ -0,0 +1,43 @@
+#!/bin/bash
+
+# BEGIN DGRN
+# -S<i1>,<i2>,...表示输出对应频率点索引值的核函数文件
+grt -Mmilrow -D2/0 -N500/0.02 -OGRN -R5,8,10 -S50,100
+# END DGRN
+
+# BEGIN grt.k2a
+grt.k2a GRN_grtstats/milrow_2_0/K_0050_5.00000e+00 > stats
+# END grt.k2a
+
+head -n 10 stats > stats_head
+echo "..." >> stats_head
+
+
+# BEGIN DEPSRC 0.0 DGRN
+grt -Mmilrow -D0/0 -N500/0.02 -OGRN -R5,8,10 -S50,100
+grt.k2a GRN_grtstats/milrow_0_0/K_0050_5.00000e+00 > stats
+# 绘制图像部分见Python
+# END DEPSRC 0.0 DGRN
+
+# BEGIN grt.k2a ptam
+grt.k2a GRN_grtstats/milrow_0_0/PTAM_0002_1.00000e+01/PTAM_0050_5.00000e+00 > ptam_stats
+# END grt.k2a ptam
+
+head -n 10 ptam_stats > ptam_stats_head
+echo "..." >> ptam_stats_head
+
+
+
+
+
+
+
+# BEGIN SGRN
+# -S 表示输出核函数文件
+stgrt -Mmilrow -D0.1/0 -X2/2/1 -Y2/2/1 -S > grn
+
+# grt.k2a 也可以读取静态解输出的核函数文件，格式一致
+grt.k2a stgrtstats/milrow_0.1_0/K > static_stats
+
+# 绘制图像部分见Python
+# END SGRN
diff --git a/docs/source/Tutorial/integ_converg/run/stats_head b/docs/source/Tutorial/integ_converg/run/stats_head
new file mode 100644
index 00000000..42319c23
--- /dev/null
+++ b/docs/source/Tutorial/integ_converg/run/stats_head
@@ -0,0 +1,11 @@
+#                k                            EXP_q0                            EXP_w0                             VF_q0                             VF_w0                             HF_q1                             HF_w1                             HF_v1                             DC_q0                             DC_w0                             DC_q1                             DC_w1                             DC_v1                             DC_q2                             DC_w2                             DC_v2
+    4.18879020e-02   -2.76036558e-01+5.35711450e-01J   -1.01611310e+02-3.61631266e+01J   -1.22267892e-01-6.70379186e-02J    4.96588630e+00-1.37205887e+01J   -2.21729097e+01+7.37515638e+00J    9.73448946e-02+1.00798198e-01J    2.21734994e+01-7.37567188e+00J   -3.33840313e+00+1.99821238e+00J   -2.03210388e+02-7.23135866e+01J   -6.32239892e+01-2.31293889e+02J   -1.80498128e+00+1.32938528e+00J    6.32368482e+01+2.31302607e+02J   -9.28776671e-01+3.08929828e-01J    4.07757341e-03+4.22222504e-03J    9.28801369e-01-3.08951421e-01J
+    8.37758041e-02   -5.53773641e-01+1.06959884e+00J   -1.01564466e+02-3.62216360e+01J   -2.44884392e-01-1.34651089e-01J    4.97974884e+00-1.37129492e+01J   -2.21736340e+01+7.36677560e+00J    1.94765330e-01+2.01748618e-01J    2.21759999e+01-7.36883678e+00J   -6.68038934e+00+3.99067034e+00J   -2.03079982e+02-7.23925670e+01J   -6.31209463e+01-2.31301596e+02J   -3.61794113e+00+2.64994395e+00J    6.31723448e+01+2.31336782e+02J   -1.85761402e+00+6.17157549e-01J    1.63166221e-02+1.69016527e-02J    1.85781222e+00-6.17330226e-01J
+    1.25663706e-01   -8.34876166e-01+1.59980148e+00J   -1.01486568e+02-3.63185118e+01J   -3.68183235e-01-2.03423546e-01J    5.00272203e+00-1.37000275e+01J   -2.21748220e+01+7.35281473e+00J    2.92334892e-01+3.03005240e-01J    2.21801728e+01-7.35744894e+00J   -1.00294633e+01+5.97154881e+00J   -2.02862928e+02-7.25227932e+01J   -6.29490358e+01-2.31313427e+02J   -5.44662611e+00+3.95266224e+00J    6.30645283e+01+2.31393780e+02J   -2.78657032e+00+9.23981950e-01J    3.67358859e-02+3.80767614e-02J    2.78724272e+00-9.24564302e-01J
+    1.67551608e-01   -1.12093624e+00+2.12438694e+00J   -1.01377906e+02-3.64528030e+01J   -4.92468171e-01-2.73956193e-01J    5.03460210e+00-1.36815458e+01J   -2.21764452e+01+7.33328524e+00J    3.90123305e-01+4.04725463e-01J    2.21860261e+01-7.34151417e+00J   -1.33889696e+01+7.93488509e+00J   -2.02559714e+02-7.27021689e+01J   -6.27080448e+01-2.31327860e+02J   -7.29831179e+00+5.22816351e+00J    6.29129307e+01+2.31473650e+02J   -3.71569905e+00+1.22870374e+00J    6.53657871e-02+6.78124021e-02J    3.71730435e+00-1.23008251e+00J
+    2.09439510e-01   -1.41343288e+00+2.64131923e+00J   -1.01238944e+02-3.66231977e+01J   -6.17997403e-01-3.46873476e-01J    5.07508972e+00-1.36571277e+01J   -2.21784636e+01+7.30820477e+00J    4.88194640e-01+5.07071590e-01J    2.21935712e+01-7.32104012e+00J   -1.67620054e+01+9.87451894e+00J   -2.02171147e+02-7.29277930e+01J   -6.23977665e+01-2.31342773e+02J   -9.17955069e+00+6.46655569e+00J    6.27168941e+01+2.31576449e+02J   -4.64504655e+00+1.53062683e+00J    1.02247246e-01+1.06200825e-01J    4.64821068e+00-1.53331506e+00J
+    2.51327412e-01   -1.71368751e+00+3.14843202e+00J   -1.01070374e+02-3.68280460e+01J   -7.44967511e-01-4.22827939e-01J    5.12377698e+00-1.36263113e+01J   -2.21808268e+01+7.27759878e+00J    5.86605699e-01+6.10212065e-01J    2.22028232e+01-7.29603592e+00J   -2.01513244e+01+1.17840443e+01J   -2.01698457e+02-7.31960030e+01J   -6.20180894e+01-2.31355451e+02J   -1.10959030e+01+7.65731424e+00J    6.24755683e+01+2.31702227e+02J   -5.57464981e+00+1.82906007e+00J    1.47430092e-01+1.53363019e-01J    5.58017810e+00-1.83369383e+00J
+    2.93215314e-01   -2.02281428e+00+3.64340910e+00J   -1.00873179e+02-3.70653962e+01J   -8.73496901e-01-5.02501935e-01J    5.18013215e+00-1.35885678e+01J   -2.21834742e+01+7.24150249e+00J    6.85404670e-01+7.14322424e-01J    2.22138008e+01-7.26651155e+00J   -2.35592317e+01+1.36567765e+01J   -2.01143445e+02-7.35024415e+01J   -6.15691096e+01-2.31362617e+02J   -1.30516479e+01+8.78921387e+00J    6.21879084e+01+2.31851007e+02J   -6.50453437e+00+2.12331943e+00J    2.00971146e-01+2.09450274e-01J    6.51342659e+00-2.13065247e+00J
+    3.35103216e-01   -2.34166518e+00+4.12377419e+00J   -1.00648719e+02-3.73330466e+01J   -1.00360904e+00-5.86605825e-01J    5.24348343e+00-1.35433257e+01J   -2.21863362e+01+7.19996335e+00J    7.84630151e-01+8.19585973e-01J    2.22265269e+01-7.23247713e+00J   -2.69874682e+01+1.54857410e+01J   -2.00508641e+02-7.38421555e+01J   -6.10512612e+01-2.31360474e+02J   -1.50494788e+01+9.85032232e+00J    6.18526730e+01+2.32022766e+02J   -7.43471262e+00+2.41273088e+00J    2.62932087e-01+2.74645896e-01J    7.44818066e+00-2.42362635e+00J
+    3.76991118e-01   -2.67077029e+00+4.58689283e+00J   -1.00398809e+02-3.76286187e+01J   -1.13521600e+00-6.75872023e-01J    5.31300339e+00-1.34900004e+01J   -2.21893357e+01+7.15304371e+00J    8.84310647e-01+9.26194250e-01J    2.22410284e+01-7.19394215e+00J   -3.04370881e+01+1.72636875e+01J   -1.99797487e+02-7.42097363e+01J   -6.04654610e+01-2.31344778e+02J   -1.70901870e+01+1.08280697e+01J    6.14684235e+01+2.32217404e+02J   -8.36518249e+00+2.69663395e+00J    3.33377260e-01+3.49167006e-01J    8.38467018e+00-2.71205230e+00J
+...