在Approximate Molecular Orbital Theory by Pople John A., Beveridge David L. 的182頁出現了BINTGS函數，這個函數只有一個作用就是計算積分

用與AINTGS一樣的分部積分法

與之對應的遞推程序是

static double []b=new double[17]; public static void BINTGS120( double ρ, int k, double []b ) throws IOException, ParseException {double expx=Math.exp(ρ);double expmx=1/expx;b[0]=(expx-expmx)/ρ; for(int i=1 ;i<k+1;i++) {b[i+1-1]=( i*b[i-1] +Math.pow((-1) , i )*expx-expmx ) /ρ;System.out.println( b[i+1-1] +" ** "+i );}}

驗算

BINTGS120( 2, 3, b?? )= -1.2450533100532868

wolframalpha

Integrate[ x**3*exp(-2*x),{x,-1,1}]=

-1.245053310053286421652052233131526434740664367854223633997768300?

?

但這個方法在ρ很小時有精度損失

BINTGS120( 0.6, 13, b?? )= -0.08471914601714074

Wolframalpha

Integrate[ x**13*exp(-0.6*x),{x,-1,1}]=

-0.084304036045748143556492453678859949065231785751955121973510012...

因此BINTGS120程序僅在ρ較大的時候適用。

總結

以上是生活随笔為你收集整理的CNDO-INTGRL-SS-BINTGS-斯莱特轨道指数---递推方法的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。