以下是我程序的一部分，我尽量把变量表达的清晰点，每次运行之后得到的矩阵J（M,M）正对角线下方的数据全部为零，我不明白是为什么
public A(M,N)
public J(M,M)
public R(N)
public QQ(N)
local d,e,f
J=0
for d=1 to M
 for e=1 to M
  for f=1 to N
   if A(d,f)!=0 .and. A(e,f)!=0
    J(d,e)=J(d,e)+2*(A(d,f)*A(e,f)*R(f)*ABS(QQ(f)))
   endif
  endfor
 endfor
endfor

程序对不对因为没有判断标准就不好说。你所说的问题，也许是数据产生的正常现象吧。