1685 SDT=SDT+ND*(A(I,JE)-PY)^2