示例代码：
def quad_custom(func, a, b, N=1000):
dx = (b - a) / N
integral = 0
for i in range(N):
x1 = a + i*dx
x2 = a + (i+1)*dx
integral += func(x1) * dx
return integral
# 定义要积分的函数
def f(x):
return x**2
# 使用自定义的积分策略进行积分
result = quad_custom(f, 0, 1)
print(result)
在上面的代码中，quad_custom 函数接受一个函数 func、积分区间 [a, b] 和可选的积分精度 N，并返回对函数 func 在 [a, b] 区间上的积分值。这个函数使用矩形法计算积分，将积分区间分成 N 个小矩形，并计算每个小矩形的面积，然后将它们相加得到总的积分值。
在示例中，我们定义了一个函数 f(x) = x**2，然后使用自定义的积分策略 quad_custom 对这个函数在区间 [0, 1] 上进行积分。最后打印出积分结果。
