1580 SXY=SXY+(A(I,J)-PX)*(A(IS,J)-PX)