python代码如下:
import math
### 第一类反常积分(无穷区间)数值分析
### y = 1/( x^2 )
### 积分区间[1,+inf)
def Func(x):...return 1/ pow(x,2)
def Improp1(Func, a, inf, eps = 1e-6):
### [a,+inf) (inf > 0时) 或者 (-inf,a](...*math.fabs(s) or math.fabs(Func(x)) > eps):
s1 = AdaptiveTrapzCtrl(Func,x,x+math.fabs(h),eps)...*math.fabs(t)): break
if (math.fabs(t) <= eps and math.fabs...return t
s = Improp1(Func, 1, 10, eps = 1e-6)
print(s)
输出结果
?