hello Mr. Holoborodko. for example this function I can create:

>>myfun = @(x)arrayfun( @(x) integral(@(t) t,eps, x, 'RelTol', eps), x);
>>myfun([1 2])
ans =
   0.500000000000000   2.000000000000000

but with your toolbox:

>>mufun = @(x)arrayfun( @(x) quadgk(@(t) t,100*eps('mp'), x, 'RelTol', 100*eps('mp')), x);
>>mufun(mp([1 2]))
ans = 
    [1x1 mp]    [1x1 mp]

however when I try to give scalar without arrayfun, it works:

mufun = @(x) quadgk(@(t) t,100*eps('mp'), x, 'RelTol', 100*eps('mp'));
mufun(mp(2))
ans = 
    2.000000000000000000000000000000001

what does it mean when it displays [1x1 mp] ?