1610 SD(3)=SD(3)+(TJ(M,N)-TJ(M,NA)-TJ(MA,N)+PY)^2