ClrHome
rand180-90→A
rand360-180→B
Disp "LAT",A
Disp "LON",B
3→dim(L₁)
For(C,1,3)
rand60+15→H
H→L₁(C)
End
3→dim(L₂)
rand360→L₂(1)

L₂(1)+rand30+90→L₂(2)
If L₂(2)>360
L₂(2)-360→L₂(2)

L₂(2)+rand30+90→L₂(3)
If L₂(3)>360
L₂(3)-360→L₂(3)
3→dim(L₃)
3→dim(L₄)
For(E,1,3)
sin⁻¹(sin(A)*sin(L₁(E))+cos(L₂(E))*cos(L₁(E))*cos(A))→D
D→L₃(E)
cos⁻¹(((sin(L₁(E))-sin(L₃(E))*sin(A))/cos(L₃(E))*cos(A)))→G
G→L₄(E)
End
3→dim(L₅)
For(F,1,3)
If L₂(F)>90
L₄(F)*(⁻1)→L₄(F)
B+L₄(F)→L₅(F)
End
{3,3}→dim([J])
List►matr(L₅,L₃,L₁,[J])
Disp [J]