The intermediate integration is numerical in nature so the operation needs to be vectorized. Unfortunately the integral function does not handle vectors for integral limits. Try the following
F = @(t) cos(t);
A = @(t) arrayfun(@(s)-integral(F,-10,s),s);
Z = @(t1,t2) -integral(A,t1,t2);
Best Answer