Te falta un ) de cierre para el primer 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
Aunque es difícil de detectar visualmente, encontré esto al copiar su código en un editor de texto que admita la combinación de llaves. Se recomienda encarecidamente utilizar uno, si no para el desarrollo y la prueba de consultas, al menos para la depuración.
