Returns haversine of x (angle in radians).
x
hav(x) == (1 - cos(x)) / 2 == sin(x / 2)².
static double hav(double x) { final sine = sin(x * 0.5); return sine * sine; }