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)); }