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