Matlab积分Mindlin解

Matlab积分Mindlin解
我在用Matlab二重积分Mindlin解的时候,出现无法计算的问题.
syms a c;
for i=1:n
for j=1:n
z=d1*i-d1/2;
R1=sqrt((d*cos(a))^2+(z-c)^2);
R2=sqrt((d*cos(a))^2+(z+c)^2);
A=(1+u)/(8*pi*(1-u))*((3-4*u)/R1+(z-c)^2/R1^3+(8*u^2-12*u+5)/R2+((3-4*u)*(z+c)^2-2*c*z)/R2^3+6*c*z*(z+c)^2/R2^5);
I1(i,j)=2*int(int(A,c,(j-1)*d1,j*d1),a,0,pi/2);
end
end
其中n,d1,d,u都是可以给出的,但是为什么二重积分没办法计算呢?