CNDO-INTGRL-SS-AINTGS-斯莱特轨道指数

在Approximate Molecular Orbital Theory by Pople John A., Beveridge David L. 的182页出现了AINTGS函数，这个函数只有一个作用就是计算积分

 

 

因此这个函数的积分程序是

static double []a3=new double[17];

public static  void AINTGS( double ρ, double k  ) throws IOException, ParseException {
				
				a3[1-1]=Math.exp(-ρ)/ρ;
				
				for(int n=1 ;n<k+1 ;n++) {
					
					a3[n+1-1]=( a3[n-1]*(n)+Math.exp(-ρ) )/ρ;
					
					//System.out.println(  a3[n+1-1] +" ***  "+n );
				}
				
			}

验算

AINTGS(  2 , 3   )= 0.32142129768695515

print ( N(sympy.integrate(  x**3*exp(-2*x)  ,( x,1,oo   ) ) )  )= 0.321421297686955