code: static double sinSumFromHav(double x, double y) {
          double a = sqrt(x * (1 - x));
          double b = sqrt(y * (1 - y));
          return 2 * (a + b - 2 * (a * y + b * x));
      }
docs: Returns sin(arcHav(x) + arcHav(y)).