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

xulinzero

function MultiAR(P)
clear
% 多维多变量AR法合成地震动场程序
% AR法计算模型 AR(p)
% p为AR模型的阶数
% Np为空间相关的点的数目
p=input('Order of AR p = ')
Np=input('number of the points in the wind field = ')
Td=input
% coor.txt为所模拟随机地震场的空间相关点的序号和坐标
% 第一列为空间点编号，第二至四列分别为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:p
    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.”

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

[ 本帖最后由 zhangnan3509 于 2008-11-20 08:53 编辑 ]