$2 = H * W
while($1 < $2)
  $3 = ($1 / W) + 1
  $4 = ($1 % W) + 1
  $5 = (A[$3,$4] - B[$3,$4])*(A[$3,$4] - B[$3,$4])+1
  $6 = $5 / ($5 * $5)
  $5 = 1 - $6
  $7 = $7 + ($3 * $5)
  $8 = $8 + ($4 * $5)
  $1 = $1 + 1
end
return($7,$8)