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

littlestupid

同一个版本7.11.0.584 (R2010b)，安装在一个32位台式机上，和一个64位笔记本上，操作系统都是win7。
* l# B: Z8 r6 \( [+ y7 t  _$ j同一个.m文件，在32位台式机上运行的好好的，但是在64位笔记本上就出错了！报错如下：
??? Error using==> mtimes
+ @+ ?) j+ {7 R) `/ Y/ sInner matrixdimensions must agree.
Error in ==>evalincaller at 14
r =evalin('caller',a);
Error in ==>sym.eval at 33
  r* h- ^$ p  x" G9 J9 S  a' dr =evalin('caller',['evalincaller(''' sym2str(a) ''')']);
Error in ==>ColdHRBWO7 at 247
   plot(rx,eval(SS(kn,1))/TP,'k','LineWidth',sqrt(abs(kn)));
. z# V' Y; r9 Z8 G经过检查，问题出在：
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));
0 v+ E( z- ^+ |5 R* G上。于是，我做了如下尝试：
: S6 {: E& [& W2 O4 f7 K: z8 ?(1)     我改成：
0 s& Q* J# u9 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));
在64位笔记本上就不出错了。但真正的公式里面，的确需要两个bessel公式相乘。
  w, d# M: y  ^: o1 B- Y' r6 Q! N(2)     我又把两个bessel之间的.*改成*：
3 `" g6 E) @: A( x9 tSS(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));
. C0 d9 S3 o" U( u6 g2 s; Q64位笔记本还出错！报的错误仍然是“Inner matrix dimensionsmust agree.”
5 u3 K: A0 B6 {+ d, d(3)     我又改成如下两种形式：
3 r7 T$ e- Z7 B0 d8 U$ ?5 h2 `2 n6 qSS(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$ S6 z% j6 k3 P8 ?, C7 e) 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)).^2)
64位笔记本继续出错，报的错误又成了：“Inputs must be a scalarand a square matrix.”

出问题的程序段如下：
2 R( P( q8 T/ ?syms r z;
/ H6 [( ?6 a' N: ?' V7 n  ~SS = zeros(2*n+1,1);
! x/ u' y8 F+ \2 rSS = sym(SS);
, e. q: G: K8 E" ofor 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));
- f; x% d) g* Y# W7 E1 eend
%#########################谐波功率#########################;
& B' {/ y/ i2 q' B' l+ @power=zeros(2*n+1,1);
TP = 0;
for kn=1:1:2*n+1
    clear sp;
) v* q" U6 [9 O# B, N6 ~3 I    sp = @(z,r)real(besselj(1,T(kn,1).*r).*besselj(1,T(kn,1).*r).*r);
    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));
" o; [3 C$ s% T; v0 r    TP = TP + power(kn,1);
end
power = power/TP;
3 H- n8 e& p' Dsave power;
xlswrite('5-谐波功率.xls',power);
%#########################功率密度#########################;
rx = 0:0.01：(ro-d)*1e3;
' n8 l4 V( p* `  M4 er = rx*1e-3;
, w2 L! z" L# y% b9 [" d5 Wz = 0.5*p;
GLMD = rx';
! n0 s9 u4 c5 K( ~. n% k( j9 w; mfigure(100);
, `) O2 e8 z* W* zfor kn=1:1:2*n+1
. Q% A2 J/ j$ l8 e% p    clear GLMDn;
7 o+ M- X6 z5 C, r    plot(rx,eval(SS(kn,1))/TP,'k','LineWidth',sqrt(abs(kn)));
" [1 K6 a( ^4 G& W2 w    hold on;
& _+ p: ~7 L1 Y' b8 }0 v9 |    grid on;
6 @: Z3 ^6 Q7 u8 L    box on;
    xlabel('r-axis (mm)');
% G7 |, ]4 w, j1 c4 u; T) s- Y    ylabel('Power flux density (W/m2)');
2 r3 X! ^4 I* s7 {6 t0 r6 Q3 y/ B    GLMDn = eval(SS(kn,1))/TP;
8 Y" d! |/ V3 E6 y% C    GLMDn = GLMDn';
    GLMD = [GLMD,GLMDn];
- x  |: V) x" M' O) N8 k" Bend
save GLMD;
$ z" K3 P3 }: ]; W/ l/ M+ nxlswrite('5-Poynting.xls',GLMD);

ExxNEN

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