错误使用 integralCalc/finalInputChecks 输入函数必须返回 ‘double‘ 或 ‘single‘ 值。报错

a=0.5;
p=0.5;
r=0.03;
q=0.03;
s=0.3;
t=0;
T=1;
v=2;
tau=0.0056;
pai=3.14;
z=0.5;
x=1;
syms rho n;
u=@(rho)double(real(exp(-(v^2+rho.^2).*tau./2).*(2*x)^(1/p+(1-v)/2).*exp(1/(4*x))*pai./sin(1i.*rho.*pai).*z.^((1-rho.*1i)./2).*exp(-z/2).*symsum(igamma((v-rho.*1i)./2+n,0)./igamma((v-rho.*1i)./2,0)./igamma(1-rho.*1i+n,0).*igamma(1-rho.*1i,0)*z^n/factorial(n),n,0,15)./igamma((v+rho.*1i)./2,0)-z.^((1+rho.*1i)./2).*exp(-z/2).*symsum(igamma((v-rho.*1i)./2+n,0)./igamma((v+rho.*1i)./2,0)./igamma(1+rho.*1i+n,0).*igamma(1+rho.*1i,0)*z^n/factorial(n),n,0,15)./igamma((v-rho.*1i)./2,0).*(real(igamma((v+rho.*1i)./2,0)).^2+imag(igamma((v+rho.*1i)./2,0).^2).*(real(igamma(1/p-(v+rho.*1i)./2,0)).^2+imag(igamma(1/p-(v+rho.*1i)./2,0)).^2).*sin(a*(rho.*pai).^(-p)))));
%%%%%%积分的被积式子
%fplot(u,[0,15]);
re=integral(u,0,10);

在计算下图的较为复杂的积分时：

总会出现如下报错：

错误使用 integralCalc/finalInputChecks (第 522 行)
输入函数必须返回 'double' 或 'single' 值。找到 'sym'。

出错 integralCalc/iterateScalarValued (第 315 行)
                finalInputChecks(x,fx);

出错 integralCalc/vadapt (第 132 行)
            [q,errbnd] = iterateScalarValued(u,tinterval,pathlen);

出错 integralCalc (第 75 行)
        [q,errbnd] = vadapt(@AtoBInvTransform,interval);

出错 integral (第 87 行)
Q = integralCalc(fun,a,b,opstruct);

出错 test (第 17 行)
re=integral(u,0,10);

对目标函数画出的图像如下图：

也会显示有数组符号的问题

目前已经做过很多修改尝试，均未成功实现积分（包括int函数、trapz函数）

希望可以求出积分的值，谢谢大家！！！