1985 DIAG=DIAG+2*TI(KI,KJ)*TI(KI,KJ)