同一个版本的matlab、同一个.m文件，为何一个顺利执行、另一个出错？

littlestupid

同一个版本7.11.0.584 (R2010b)，安装在一个32位台式机上，和一个64位笔记本上，操作系统都是win7。
. |) H* r+ y8 O: A4 W1 R0 |3 x% P同一个.m文件，在32位台式机上运行的好好的，但是在64位笔记本上就出错了！报错如下：
. U1 C" R( M( Z??? Error using==> mtimes
Inner matrixdimensions must agree.
Error in ==>evalincaller at 14
r =evalin('caller',a);
Error in ==>sym.eval at 33
r =evalin('caller',['evalincaller(''' sym2str(a) ''')']);
Error in ==>ColdHRBWO7 at 247
   plot(rx,eval(SS(kn,1))/TP,'k','LineWidth',sqrt(abs(kn)));
% E4 F: @, h, `; U4 o+ Z经过检查，问题出在：
# V/ s7 n, ^0 Y  W$ c# c) VSS(kn,1) =real(omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*besselj(1,T(kn,1).*r).*besselj(1,T(kn,1).*r));
6 `4 K% ^6 Q& k* c上。于是，我做了如下尝试：
(1)     我改成：
SS(kn,1) =real(omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*besselj(1,T(kn,1).*r));
在64位笔记本上就不出错了。但真正的公式里面，的确需要两个bessel公式相乘。
(2)     我又把两个bessel之间的.*改成*：
: C" G  F" m( M/ aSS(kn,1) =real(omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*besselj(1,T(kn,1).*r)*besselj(1,T(kn,1).*r));
64位笔记本还出错！报的错误仍然是“Inner matrix dimensionsmust agree.”
(3)     我又改成如下两种形式：
SS(kn,1) =real(omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*(besselj(1,T(kn,1).*r))^2);
SS(kn,1) =real(omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*(besselj(1,T(kn,1).*r)).^2)
4 x! h! f( \6 D3 p) N64位笔记本继续出错，报的错误又成了：“Inputs must be a scalarand a square matrix.”

出问题的程序段如下：
4 C& m0 L9 Y4 t* |! x/ zsyms r z;
4 }9 v' J" {1 D% K* ~. W5 }SS = zeros(2*n+1,1);
SS = sym(SS);
for kn=1:1:2*n+1
    SS(kn,1) =real(omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*besselj(1,T(kn,1).*r).*besselj(1,T(kn,1).*r));
end
2 K0 X0 a5 B: x7 H%#########################谐波功率#########################;
/ w5 @& l$ c) n; ~8 U; a+ Bpower=zeros(2*n+1,1);
TP = 0;
# i" J; |0 e3 ufor kn=1:1:2*n+1
    clear sp;
0 c! J0 W( D- T9 j. q' J( e$ K    sp = @(z,r)real(besselj(1,T(kn,1).*r).*besselj(1,T(kn,1).*r).*r);
7 O- j9 ]6 k1 R% u1 Y1 g    power(kn,1) =real(2*pi/p*omg*epsl*T(kn,1)*T(kn,1)*bt(kn,1)*A(kn,1)*A(kn,1)*quad2d(sp,0,0.5*p,0,Rw,'RelTol',1e-5,'MaxFunEvals',9000));
    TP = TP + power(kn,1);
end
power = power/TP;
save power;
xlswrite('5-谐波功率.xls',power);
" @6 B$ A* O* k; ~7 T0 {%#########################功率密度#########################;
rx = 0:0.01：(ro-d)*1e3;
" r& i+ ?/ z& z. Y6 g; r% ]' Dr = rx*1e-3;
z = 0.5*p;
GLMD = rx';
$ C7 W/ x4 }. t: bfigure(100);
+ c7 R" ^/ x' {7 W; o; @for kn=1:1:2*n+1
  j9 p8 z0 K# B4 l5 L    clear GLMDn;
- V5 j0 a0 C6 j    plot(rx,eval(SS(kn,1))/TP,'k','LineWidth',sqrt(abs(kn)));
    hold on;
0 [) \0 c0 w- P9 C3 [    grid on;
9 V9 V6 S/ c  t/ C' W7 I$ y    box on;
, S& s8 r" j) V    xlabel('r-axis (mm)');
1 g  E! E3 L% I( h3 R/ C: y    ylabel('Power flux density (W/m2)');
    GLMDn = eval(SS(kn,1))/TP;
4 U' t, B0 [+ ~    GLMDn = GLMDn';
% v* L+ {$ k, c( z& {    GLMD = [GLMD,GLMDn];
end
save GLMD;
# ^  v, v# r  o0 k9 r1 l) X, h1 |xlswrite('5-Poynting.xls',GLMD);

ExxNEN

我的天，把两个文件放在同一目录下，在一个M文件中就可以直接使用另一个文件的函数了啊