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

littlestupid

同一个版本7.11.0.584 (R2010b)，安装在一个32位台式机上，和一个64位笔记本上，操作系统都是win7。
同一个.m文件，在32位台式机上运行的好好的，但是在64位笔记本上就出错了！报错如下：
/ f4 K' R9 p/ E6 A3 F5 |??? Error using==> mtimes
! J. e8 U3 W9 a; [5 ]& K( r+ NInner matrixdimensions must agree.
4 T% W% I/ }! R4 _9 ~6 h: cError in ==>evalincaller at 14
r =evalin('caller',a);
+ G' Q0 Q1 j! m( J# O8 {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)));
经过检查，问题出在：
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));
上。于是，我做了如下尝试：
(1)     我改成：
! A& a7 ^7 _4 c; H3 y: K6 A' e; ZSS(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));
7 G- b' e  i) E: }在64位笔记本上就不出错了。但真正的公式里面，的确需要两个bessel公式相乘。
(2)     我又把两个bessel之间的.*改成*：
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));
+ u7 b) U6 w( I4 w64位笔记本还出错！报的错误仍然是“Inner matrix dimensionsmust agree.”
8 M, l8 h' e3 m" t(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);
: v, a) O  R* L& _7 iSS(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)
7 H& a1 n. F0 q* v/ J& H64位笔记本继续出错，报的错误又成了：“Inputs must be a scalarand a square matrix.”

出问题的程序段如下：
syms r z;
SS = zeros(2*n+1,1);
& G$ {- ^6 t2 L7 f/ ZSS = 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));
' G: [, z1 }9 e9 Cend
%#########################谐波功率#########################;
power=zeros(2*n+1,1);
8 r5 b4 ^. k( r- P" {" STP = 0;
( V& P0 O# ?( `6 i" b( {( Rfor kn=1:1:2*n+1
$ t8 N% N! m3 p+ \' Y+ ^    clear sp;
    sp = @(z,r)real(besselj(1,T(kn,1).*r).*besselj(1,T(kn,1).*r).*r);
5 q7 a. ?* j8 J7 I1 D; p    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));
6 E, q$ k3 v6 r  e    TP = TP + power(kn,1);
8 O9 T7 P* Y6 o# u6 Y5 F+ send
power = power/TP;
save power;
xlswrite('5-谐波功率.xls',power);
%#########################功率密度#########################;
2 l, S9 U7 v- H! a+ z" {2 Trx = 0:0.01：(ro-d)*1e3;
. j3 x* Z. U; [8 k- o. A- _$ er = rx*1e-3;
z = 0.5*p;
GLMD = rx';
; K& a, v3 R5 t1 p* @/ Hfigure(100);
for kn=1:1:2*n+1
    clear GLMDn;
    plot(rx,eval(SS(kn,1))/TP,'k','LineWidth',sqrt(abs(kn)));
    hold on;
9 T" I  }0 M7 d    grid on;
    box on;
* b, T9 a- g0 @$ l4 Y! u) K; D    xlabel('r-axis (mm)');
; V2 a+ u0 f# b: O0 \# I    ylabel('Power flux density (W/m2)');
    GLMDn = eval(SS(kn,1))/TP;
    GLMDn = GLMDn';
: Q, |+ a- o+ z! h8 V+ T- t    GLMD = [GLMD,GLMDn];
3 D- Q9 i/ s, m' _  pend
+ n. e" T! t# p1 Wsave GLMD;
xlswrite('5-Poynting.xls',GLMD);
3 }" z1 J* a. n1 t7 L

ExxNEN

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