急啊！高人帮我看下程序错在哪。

[xulin]

function MultiAR(P)
clear
p=input;
Np=input('number of the points in the wind field = ')
Td=input;
% 第一列为点编号，第二至四列分别为X,Y,Z三个方向的坐标
load 'coor.txt'
XYZ=coor;
X=coor(:,2);
Y=coor(:,3);
Z=coor(:,4);
%%%%%%%
fg=2.5;
kg=0.60;
s0=0.6;
A=0.736;
a=0.147;
vapp=2.5;
RR=[];
syms f
for s=0
for n=1:Np
for m=1:n
d=sqrt((X(n)-X(m)).^2+(Y(n)-Y(m)).^2+(Z(n)-Z(m)).^2);
D=sqrt((X(n)-X(m)).^2+(Y(n)-Y(m)).^2);
q=3300/(sqrt(1+(4*f/3).^2));
F=@(f)(s0*(1+(2*kg*f./fg).^2)*(A*exp((-2)*d*(1-A+a*A)./(a*q))+(1-A)*exp((-2).*d.*(1-A+a*A)./q))*cos(2*pi*f*(s.*dt-D./vapp)))./((1-(f./fg).^2).^2+(2*kg*f./fg).^2);
Q=quadl(F,0,10);
Q(n,m)=Q(m,n);
end
end
RR=[RR Q];
end


高人请指教啊 ，程序函数F=@(f)(...)怎么改都报错
“??? Error using ==> sym.mtimes
Inner matrix dimensions must agree.”

现在很急，还请大家帮帮我啊，在下感激涕零！！！