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