1810 DIAG=DIAG+TI(I,I)*TI(I,I)