matlab程序总是正忙，显示不出结果，求教大神

Question

最终显示Bx，急求
syms f w
a=0.04;
R=0.02;
p=0.04;
Bx=0;
z=0;
I=100;
u=4*pi*10^-7;
x=p.*cos(w);
y=p.*sin(w);
b=a-R.*cos(f);
for k=0:2
for n=1:4
t=pi/4+pi*k/3+pi*(n-1)/2;
Xkn=(a-R.*cos(f)).*cos(f+t);
Ykn=(a-R.*cos(f)).*sin(f+t);
Zkn=R.*sin(f);
d=sqrt(((x-Xkn)^2+(y-Ykn)^2+(z-Zkn)^2)^3);
S1=((z-R*sin(f)).*(a*cos(f+t)-R*cos(2*f+t))-(y-b*sin(f+t)).*R*cos(f))/d;
S2=u*I/4*pi;
S3=int(S1,f,-pi/4,pi/4);
Bx=S2.*S3+Bx
end
end
这是程序

Answer 1

因为S1表达式含有两个f、w变量，所以用int()函数求S1(f,w)是错误的。Matlab分不清是求以f为变量的积分，还是以w为变量的积分。电脑始终在瞎转。。。
请你把要求的积分式贴出来，以便于帮你解决。

追问

