static double distance(double ax, double ay, double bx, double by) {
  ret sqrt(sqr(bx-ax)+sqr(by-ay));
}