我想为不同数量的T计算以下积分int[0,T] int[0,x] F(T,x,y)dy dx，比如T=1:0.25:11。我假设int[0,x] F(T,x,y)dy是函数g(x)并应用trapz来帮助评估迭代积分，以下是我的代码：
M=T.'*linspace(0,1,100);
fun=@(T,x)integral(@(y)F(T,x,y),0,x);
Points=bsxfun(fun,T',M);
Ans=1/100*trapz(Points,2);
但是，第三行给出了一条错误消息：A and B must be floating point scalars.我理解错误消息，但在我看来，我没有发现代码中的错误。任何人都可以解决这个问题吗？