可以通过 for 循环和 ty_integral2 函数配合使用实现。

示例： f=1/（n+x^2+y^2），n 取 1~10 的情况下，x 的范围是 0 到 y， y 积分,y 的范围是 0 到 1。

1.运行下面代码：

using TyMath
for n = 1:10
    global q,e
    fun(y, x) = @. 1 / n + x^2 + y^2
    xmax(y) = y
    q,e = ty_integral2(fun, 0, 1, 0, xmax)
end

2.在命令行中输入 q、e，可以看到实现了二重积分和向量化积分的叠加运算：