MATLAB小程序一个，大家来试试
 

输入参数为ε1，ε2，ε3，εd
ε1e=ε1-εd
ε3e=ε3-εd
γ＝－ε3e/ε1e
Q1＝ε1e（3γ－1）/(3-γ)+ ε3e
Q3＝ε3e（3γ－1）/(3-γ)+ ε1e
ε1p=0.4 *Q3+78.4
ε3p=-0.2 *Q3-39.3
ε1e(1)=ε1e-ε1p
ε3e(1)=ε3e-ε3p
迭代到ε1e(n)与ε1e(n-1)及ε3e(n)与ε3e(n-1)接近为止1