求积分的数值解遇到错误

求积分，但是显式积分求不出来，我打算用数值解，如下
syms a z x
y=@(x)x.((x+a).^2+z.^2).^(-5/2).(4xa) %a,z都是常数
quad(y,0,1)

遇到如下错误,请问这两个错误怎么解决
??? Error using ==> sym.mtimes
Inner matrix dimensions must agree.

Error in ==> @(x)x*((x+a).^2+z.^2).^(1/2)(4x*a)

Error in ==> quad at 63
y = f(x, varargin{:});

a和z都是常数的话就不要定义为符号变量，比如a和z都等于1，你的代码可以改写如下：

y = @(x,a,z) x .* ((x + a).^2 + z.^2).^(1/2) .* (4 .* x .* a);
q = integral(@(x) y(x,1,1), 0, 1);

运行结果：
q =

2.6908

码字不易，有用麻烦点一下采纳哦

syms a z x
y = @(x) x * ((x + a).^2 + z.^2).^(1/2) * (4 * x * a);
integral(y, 0, 1);