1770 SXY=SXY+(A(I,J)-PA)*(A(IS,JB)-PB)