1460 TJ(K,L)=TJ(K,L)+XP(K)*XP(L)*ND