请教MATLAB的for循环向量化问题

[zcc8607]

sym t;
for x=-0.2:0.02:0.2
for y=-0.3:0.02:0.3
for z=-0.4:0.02:0.4
Po=[x,y,z];
dl=[0,1,0];%积分元向量
Pl=[vz/2,t,-wz];%积分元向量坐标,t为积分变量
R=Po-Pl;
S=cross(dl,R)/((R*R')^(3/2));
K=double(int(S,t,-0.1,0.1));
end
end
end

请问这个问题可以怎么转化？谢谢！