如何矢量积分计算？

Question

我有一个涉及4个输入的函数：fun=@(T,x1,x2,x3)。对于每个T，我想计算[integral3(fun(1,x1,x2,x3),0,1,0,x2,0,x3),integral3(fun(2,x1,x2,x3),0,2,0,x2,0,x3)...,integral3(fun(100,x1,x2,x3),0,100,0,x2,0,x3)]

当然，我可以写一个for循环来解决这个问题。但是，我想对它进行矢量化而不是使用for循环。可以在matlab中做到吗？