Use element-wise operations for the multiplications and divisions in the anonymous functions. Also exp must have an argument.
Try this —
F = @(x,y)exp(x.*y);
ymin = @(x)log(x)./x;
ymax = @(x)1./x;
q = integral2(F,1,exp(1),ymin,ymax)
Make appropriate changes to get the desired result.
.
