package main
import (
    "fmt"
    "math"
)
func main() {
    lat1 := 29.490295
    lng1 := 106.486654
    lat2 := 29.615467
    lng2 := 106.581515
    fmt.Println(EarthDistance(lat1, lng1, lat2, lng2))
}
// 返回值的單位為米
func EarthDistance(lat1, lng1, lat2, lng2 float64) float64 {
    radius := float64(6371000) // 6378137
    rad := math.Pi/180.0
    lat1 = lat1 * rad
    lng1 = lng1 * rad
    lat2 = lat2 * rad
    lng2 = lng2 * rad
    theta := lng2 - lng1
    dist := math.Acos(math.Sin(lat1) * math.Sin(lat2) + math.Cos(lat1) * math.Cos(lat2) * math.Cos(theta))
    return dist * radius
}

<?php// 返回值的單位為米function pc_sphere_distance($lat1, $lon1, $lat2, $lon2, $radius = 6371000) {  $rad = doubleval(M_PI/180.0);  $lat1 = doubleval($lat1) * $rad;  $lon1 = doubleval($lon1) * $rad;  $lat2 = doubleval($lat2) * $rad;  $lon2 = doubleval($lon2) * $rad;  $theta = $lon2 - $lon1;  $dist = acos(sin($lat1) * sin($lat2) + cos($lat1) * cos($lat2) * cos($theta));  return $dist * $radius * 1000;}$lat1 = 29.490295;$lon1 = 106.486654;$lat2 = 29.615467;$lon2 = 106.581515;echo pc_sphere_distance($lat1, $lon1, $lat2, $lon2);