#! /usr/bin/gawk -f
# Last edited on 2004-01-20 02:46:14 by stolfi

BEGIN{ L2 = log(2); }

/^[ ]*[0-9]/ {
  amx = $1; amy = $2; bmx = $3; bmy = $4;
  printf "(%2d,%2d,%2d,%2d) = ", amx, amy, bmx, bmy;
  for (bias = 1; bias <= 8; bias*=2)
    { zax = log(amx+bias)/L2;
      zay = log(amy+bias)/L2;
      zbx = log(bmx+bias)/L2;
      zby = log(bmy+bias)/L2;
      D = zax - zay - zbx + zby;
      S = D*D;
      printf " %7.2f", S;
    }
  printf "\n";
}