Jasper Yu says to YSITD cos(lat1/pi/180) * cos(lat2/pi/180) * sin((lon2-lon1)/pi/360)^2 + sin((lat2-lat1)/pi/360)^2