Brakuje zamykającego )
dla pierwszego sin()
( 3959 * acos( cos( radians(18.204540500000) )
* cos( radians( latitude ) )
* cos( radians( longitude )
- radians(-66.450958500000) )
+ sin( radians(18.204540500000 ) ) /* <--- here */
* sin( radians( latitude ) )
)
) AS distance
Chociaż trudno to zauważyć wizualnie, znalazłem to, kopiując twój kod do edytora tekstu, który obsługuje dopasowywanie nawiasów klamrowych. Zdecydowanie zaleca się użycie jednego, jeśli nie do tworzenia i testowania zapytań, to przynajmniej do debugowania.