matlab编程mindlin解,MATLAB对mindlin解进行三重积分

利用MATLAB2017b对mindlin解进行三重积分。采用int进行多次积分时间长，结果出现复数，不是具体数值。采用integral3和其他函数又出现其余问题。作为一个初学者，看不懂出错的问题，不知如何解决。请求各位前辈们帮忙，非常感谢！

x=20.6;a=0.2;syms t y z;

r3=sqrt(x.^2+y.^2+(z-t).^2);

r4=sqrt(x.^2+y.^2+(z+t).^2);

f=@(t,y,z)t*(1-2*a)./r3.^3-t*(1-2*a)*(5-4*a)./r4.^3+3*t*x.^2./r3.^5+3*(3-4*a)*x.^2./r4.^5+4*t*(1-a)*(1-2*a)*(3-x.^2*(3*r4+z+t)./(r4.^2*(r4+z+t)))./(r4*(r4+z+t).^2)-6*t*t*(3*t-(3-2/a)*(z+t)+5*x.^2./(r4.^2))./r4.^5;

q=integral3(f,0,1,0,1,0,1)

>> Untitled

错误使用 integral2Calc>integral2t/tensor (line 231)

输入函数必须返回 'double' 或 'single' 值。找到 'sym'。

出错 integral2Calc>integral2t (line 55)

[Qsub,esub] = tensor(thetaL,thetaR,phiB,phiT);

出错 integral2Calc (line 9)

[q,errbnd] = integral