Du saknar en avslutande )
för den första 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
Även om det är svårt att upptäcka visuellt, hittade jag detta genom att kopiera din kod till en textredigerare som stöder parentesmatchning. Det rekommenderas starkt att använda en, om inte för frågeutveckling och testning, så åtminstone för felsökning.