在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