Why am I getting the "Index exceeds matrix dimensions" error when trying to perform triple integral of a function as indicated below? Thanks

h=@(x,y,z) cos(x)/sqrt((y.^2.)+(z.^2.)-2.*y.*z.*cos(x))

thirdans=triplequad(h,0,pi,1,3,1,3)

??? Index exceeds matrix dimensions.

Error in ==> quad at 85

if ~isfinite(y(7))

Answer 1

[Matt Fig]

You forgot one dot:

h=@(x,y,z) cos(x)./sqrt((y.^2.)+(z.^2.)-2.*y.*z.*cos(x));