Está faltando um
)
de fechamento para o primeiro 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
Embora seja difícil identificar visualmente, descobri isso copiando seu código em um editor de texto que suporta correspondência de chaves. É altamente recomendável usar um, se não para desenvolvimento e teste de consulta, pelo menos para depuração.