by ncforest

对矩阵积分可使用int和quadv，但当矩阵稍大时，效率便会偏低，使用循环反而更快。
int是精确的解析解，因此速度很慢，可以使用数值解quadv提高速度，使用前对矩阵先inline

>> syms x
>> fun=[sin(x),3*x;4,(cos(x)+1)]
fun =

[   sin(x),      3*x]
[        4, cos(x)+1]
>> intf=int(fun,x,0,1)

intf =

[ -cos(1)+1,       3/2]
[         4,  1+sin(1)]

>> intf2=quadv(inline(fun),0,1)

intf2 =

    0.4597    1.5000
    4.0000    1.8415

也可以使用arrayfun函数对quadl等积分函数进行矢量化运算来达到相同的目录，参见关于arrayfun的另一篇blog。