月(lat1, lon1)と上の2つのポイントの緯度/経度を指定して、haversineの式と同じ結果を与える式を(lat2, lon2)使用して、2つのポイント間の距離をキロメートルで計算します。

入力

• lat1, lon1, lat2, lon2度（角度）の4つの整数値または
• ϕ1, λ1, ϕ2, λ2ラジアン単位の4つの10進数値。

出力

2つのポイント間の距離（キロメートル単位、精度または丸められた整数）。

ハーバーシン式

どこ

• r は球体の半径です（月の半径が1737 kmであると仮定）、
• ϕ1 ラジアン単位のポイント1の緯度
• ϕ2 ラジアン単位のポイント2の緯度
• λ1 ラジアン単位のポイント1の経度
• λ2 ラジアン単位のポイント2の経度
• d 2点間の循環距離

（ソース：https : //en.wikipedia.org/wiki/Haversine_formula）

他の可能な式

• d = r * acos(sin ϕ1 sin ϕ2 + cos ϕ1 cos ϕ2 cos(λ2 - λ1)) @milesの式。
• d = r * acos(cos(ϕ1 - ϕ2) + cos ϕ1 cos ϕ2 (cos(λ2 - λ1) - 1)) @Neilの式。

入力が度であり、出力が丸められた整数である例

42, 9, 50, 2  --> 284
50, 2, 42, 9  --> 284
4, -2, -2, 1  --> 203
77, 8, 77, 8  --> 0
10, 2, 88, 9  --> 2365

ルール

• 入力と出力は、任意の便利な形式で指定できます。
• 入力が度単位かラジアン単位かを回答で指定します。
• 無効な緯度/経度値を処理する必要はありません
• 完全なプログラムまたは機能のいずれかが受け入れられます。関数の場合、出力する代わりに出力を返すことができます。
• 可能であれば、他の人があなたのコードを試すことができるように、オンラインテスト環境へのリンクを含めてください！
• 標準的な抜け穴は禁止されています。
• これはコードゴルフなので、通常のゴルフルールがすべて適用され、最短のコード（バイト単位）が勝ちます。

Wolfram言語（Mathematica）、48バイト

ArcCos[Sin@#*Sin@#3+Cos[#2-#4]Cos@#*Cos@#3]1737&

オンラインでお試しください！

式使用d = r * acos( sin ϕ1 sin ϕ2 + cos ϕ1 cos ϕ2 cos(λ2 - λ1) )場所をr = 1737

R + 地球圏、54 47バイト

function(p,q)geosphere::distHaversine(p,q,1737)

オンラインでお試しください！

入力をlongitude,latitude度単位の2要素ベクトルとして受け取ります。TIOにはgeosphereパッケージがありませんが、以下の関数と同じ結果を返すことを保証します。

7バイトを削ってくれたJonathan Allanに感謝します。

R、64バイト

function(p,l,q,k)1737*acos(sin(p)*sin(q)+cos(p)*cos(q)*cos(k-l))

オンラインでお試しください！

テストケースのように4つの入力を受け取りますが、度ではなくラジアン単位です。

JavaScript（Node.js）、65バイト

(a,b,c,d,C=Math.cos)=>1737*Math.acos(C(a-c)+C(a)*C(c)*(C(d-b)-1))

オンラインでお試しください！

Kevin Cruijssenの回答、Miles 'およびNeilのコメント、およびArnauldの要求に基づきます。

JavaScript（ES7）、90バイト

注：はるかに短いソリューションについては、@OlivierGrégoireの投稿を参照してください

TFeldの答えの直接の移植版。ラジアン単位で入力します。

(a,b,c,d,M=Math)=>3474*M.asin((M.sin((c-a)/2)**2+M.cos(c)*M.cos(a)*M.sin((d-b)/2)**2)**.5)

オンラインでお試しください！

悪名高いwith()85バイトを使用して

6バイトを節約してくれた@ l4m2に感謝

with(Math)f=(a,b,c,d)=>3474*asin((sin((c-a)/2)**2+cos(c)*cos(a)*sin((d-b)/2)**2)**.5)

オンラインでお試しください！

APL（Dyalog Unicode）、40 35 バイト^SBCS

匿名の暗黙関数。左引数として{ϕ₁、λ₁}を、右引数として{ϕ₂、λ₂}を取ります。

式2 ^r√（sin²（^{（ϕ₁-ϕ₂）} ⁄ ₂）+ cos ϕ₁ cos ϕ₂sin²（^{（λ₁–λ₂）} ⁄ ₂））を使用します。

3474ׯ1○.5*⍨1⊥(×⍨1○2÷⍨-)×1,2×.○∘⊃,¨

オンラインでお試しください！（r関数は度をラジアンに変換します）

,¨ 対応する要素を連結します。{{ϕ₁、ϕ₂}、{λ{、λ₂}}

⊃ 最初のものを選んでください。{ϕ₁、ϕ₂}

∘ それから

2×.○ それらの余弦の積; cos ϕ₁ cos ϕ₂
 ^{点灯。ドット「積」ですが、乗算の代わりにトリガー関数セレクター（2はコサイン）、プラスの代わりに時間}

1, その前に1を追加します。{1、cos ϕ₁ cos ϕ₂}

(… )× 次の関数を{ϕ₁、λϕ}および{ϕ₂、λ₂}に適用した結果を乗算します。

 - それらの違い; {ϕ₁-ϕ₂、λ₁-λ₂}

 2÷⍨ それを2で割ります。{ ^{（ϕ₁-ϕ₂）} ⁄ ₂、^{（λ₁-λ₂）} ⁄ ₂ }

 1○ そのサイン; {sin（^{（ϕ₁-ϕ₂）} ⁄ ₂）、sin（^{（λ₁-λ₂）} ⁄ ₂）}

 ×⍨ それを二乗する（点灯。自己乗算）。{sin²（^{（ϕ₁-ϕ₂）} ⁄ ₂）、sin²（^{（λ₁-λ₂）} ⁄ ₂）}

{sin²（^{（ϕ₁-ϕ₂）} ⁄ ₂）、cos ϕ₁ cos ϕ₂sin²（^{（λ₁-λ₂）} ⁄ ₂）}

1⊥ その合計（lit. base-1で評価）; sin²（^{（ϕ₁-ϕ₂）} ⁄ ₂）+ cos ϕ₁ cos ϕ₂sin²（^{（λ₁-λ₂）} ⁄ ₂）

.5*⍨ その平方根（それを半分のべき乗する）

 ¯1○ その逆正弦

 3474× これを掛けます

度単位の入力を許可する関数は次のとおりです。

○÷∘180

÷180 引数を180で割った値

○ πを掛ける

Python 2、95バイト

lambda a,b,c,d:3474*asin((sin((c-a)/2)**2+cos(c)*cos(a)*sin((d-b)/2)**2)**.5)
from math import*

オンラインでお試しください！

ラジアン単位で入力します。

I / Oがスラックされる前の古いバージョン：入力を整数度として受け取り、丸められたdistを返します

Python 2、135バイト

lambda a,b,c,d:int(round(3474*asin((sin((r(c)-r(a))/2)**2+cos(r(c))*cos(r(a))*sin((r(d)-r(b))/2)**2)**.5)))
from math import*
r=radians

オンラインでお試しください！

Java 8、113 92 88 82バイト

(a,b,c,d)->1737*Math.acos(Math.cos(a-c)+Math.cos(a)*Math.cos(c)*(Math.cos(d-b)-1))

入力a,b,c,dはϕ1,λ1,ϕ2,λ2ラジアン単位です。

@milesの短い式を使用して-21バイト。@OlivierGrégoreの
おかげで-4バイトになりました。なぜなら、とを直接使用する代わりに、すべてのas でまだ使用しているからです。@Neilの短い式 を使用して-6バイト。{Math m=null;return ...;}Math.m.returnMath

オンラインでお試しください。

説明：

(a,b,c,d)->                  // Method with four double parameters and double return-type
  1737*Math.acos(            //  Return 1737 multiplied with the acos of:
   Math.cos(a-c)             //   the cos of `a` minus `c`,
   +Math.cos(a)*Math.cos(c)  //   plus the cos of `a` multiplied with the cos of `c`
   *(Math.cos(d-b)-1))       //   multiplied with the cos of `d` minus `b` minus 1

Japt、55 50バイト

MsU *MsW +McU *McW *McX-V
ToMP1/7l¹ñ@McX aUÃv *#­7

必ずしも他の回答ほど正確ではありませんが、少年はこの回答を楽しんでいます。詳しく説明させてください。
ほとんどの言語では、この課題は非常に簡単ですが、Japtには、アークサインもアークコサインも組み込まれていないという不幸な性質があります。もちろん、JavascriptをJaptに埋め込むことができますが、それはFeng Shuiの反対です。

このマイナーな迷惑を克服するために私たちがしなければならないのは、おおよそのアークコサインであり、私たちは行ってもいいです！

最初の部分は、アークコサインに供給されるすべてのものです。

MsU *MsW +McU *McW *McX-V
MsU                        // Take the sine of the first input and
    *MsW...                // multiply by the cos of the second one etc.

結果は暗黙的に保存され、U後で使用されます。

その後、アークコサインの適切な近似値を見つける必要があります。私は怠け者で、数学が苦手なので、明らかにそれをブルートフォースするだけです。

ToMP1/7l¹ñ@McX aUÃv *#­7
T                       // Take 0
 o                      // and create a range from it
  MP                    // to π
    1/7l¹               // with resolution 1/7!.
         ñ@             // Sort this range so that
           McX          // the cosine of a given value
               aU       // is closest to U, e.g. the whole trig lot
                        // we want to take arccosine of.
                 Ã      // When that's done,
                  v     // get the first element
                    *#­7 // and multiply it by 1737, returning implicitly.

ジェネレーターの解像度に任意の数を使用できた可能性があり、手動テストで7!は十分に大きいが、かなり高速であることが示されました。

入力をラジアンとして受け取り、丸められていない数値を出力します。

おかげで5バイト削られた Oliverの。

オンラインでお試しください！

Haskell、68 66 52 51バイト

s=cos;(a!b)c d=1737*acos(s(a-c)+s a*s c*(s(d-b)-1))

オンラインでお試しください！

BMOのおかげで-1バイト

ルビー、87 70 69バイト

->a,b,c,d{extend Math;1737*acos(cos(a-c)+cos(a)*cos(c)*(cos(d-b)-1))}

オンラインでお試しください！

Kevin Cruijssenに感謝します。

ゼリー、 23 22  18 バイト

-4マイルのおかげで（の使用バイト{と}使用しながら、それらの式を。

;I}ÆẠP+ÆSP${ÆA×⁽£ġ

結果を（浮動小数点として）返すラジアン単位の[ϕ1, ϕ2,]左と[λ1, λ2]右を受け入れる2項関数。

オンラインでお試しください！

Mine ...（aを使用して1バイトも保存しました{）

,IÆẠCH;ÆẠ{Ḣ+PƊ½ÆṢ×⁽µṣ

オンラインで試す

Mapping Toolboxを備えたMATLAB、26バイト

@(x)distance(x{:})*9.65*pi

チャレンジで説明したのと同じ順序で、セル配列として4つの入力を受け取る匿名関数。

1737/180等しいため、これにより正確な結果が得られることに注意してください（月の半径が1737 kmであると仮定）9.65。

Matlab R2017bでの実行例：

Python 3、79バイト

from geopy import*
lambda a,b:distance.great_circle(a,b,radius=1737).kilometers

TIOにはgeopy.pyがありません

APL（Dyalog Unicode）、29 バイト^SBCS

完全なプログラム。stdinに{ϕ₁、st₂}を要求し、次に{λ₁、λ₂}を要求します。標準出力に出力します。

式r acos（sin ϕ₁ sin ϕ₂ + cos（λ₂–λ₁）cos ϕ₁ cos ϕ₂）を使用します。

1737ׯ2○+/(2○-/⎕)×@2×/1 2∘.○⎕

オンラインでお試しください！（r関数は度をラジアンに変換します）

⎕ {ϕ₁、ϕ₂}のプロンプト

1 2∘.○ デカルトトリガ関数アプリケーション。{{sin ϕ₁、sin ϕ₂}、{cos ϕ₁、cos ϕ₂}}

×/ 行ごとの製品。{sin ϕ₁ sin ϕ₂、cos ϕ₁ cos ϕ₂}

(… )×@2 2番目の要素で、次の要素を乗算します。

 ⎕ {λ₁、λ₂}のプロンプト

 -/ それらの違い; λ₁–λ₂

 2○ その余弦。cos（λ₁–λ₂）

{sin ϕ₁ sin ϕ₂、cos（λ₁–λ₂）cos ϕ₁ cos ϕ₂}

+/ 和; sin ϕ₁ sin ϕ₂ + cos（λ₁–λ₂）cos ϕ₁ cos ϕ₂

¯2○ その余弦。cos（sin ϕ₁ sin ϕ₂ + cos（λ₁–λ₂）cos ϕ₁ cos ϕ₂）

1737× それにrを掛けます。1737 cos（sin ϕ₁ sin ϕ₂ + cos（λ₁–λ₂）cos ϕ₁ cos ϕ₂）

度単位の入力を許可する関数は次のとおりです。

○÷∘180

÷180 引数を180で割った値

○ πを掛ける

C（gcc）、100 88 65 64バイト

88→65 @ マイルの式を使用
65→64 @ マイルの式を使用

#define d(w,x,y,z)1737*acos(cos(w-y)+cos(w)*cos(y)*(cos(z-x)-1))

オンラインでお試しください！

Excel、53バイト

=1737*ACOS(COS(A1-C1)+COS(A1)*COS(C1)*(COS(D1-B1)-1))

@Neilの式を使用します。ラジアンで入力します。

ロブスター、66バイト

def h(a,c,b,d):1737*radians arccos a.sin*b.sin+a.cos*b.cos*cos d-c

マイルの式を使用しますが、入力は度単位です。これにより、半径を乗算する前にラジアンに変換する追加のステップが追加されます。

Pythonの3、119の 103バイト

これは度を使用します。

from math import*
def f(L):a,o,A,O=map(radians,L);return 1737*acos(cos(a-A)+cos(a)*cos(A)*(cos(O-o)-1))

オンラインでお試しください！

PHP、88バイト

オリバー港の回答

function f($a,$b,$c,$d,$e=cos){return 1737*acos($e($a-$c)+$e($a)*$e($c)*($e($d-$b)-1));}

オンラインでお試しください！

SmileBASIC、60バイト

INPUT X,Y,S,T?1737*ACOS(COS(X-S)+COS(X)*COS(S)*(COS(T-Y)-1))