Peut-être s'agit-il d'un rectangle sphérique limité par deux parallèles et deux méridiens.
Attendons...
mais pourquoi n'y a t'il pas les même résultat que le tableau de @Modeste
Sub SurfaceRectangleSphérique()
Dim Latmax, Latmin, Longmax, Longmin, n&, Ldeg#, pas#, lat#, s#
Latmax = 52.714667 'degrés
Latmin = 52.310158
Longmax = 13.756256
Longmin = 12.826538
n = 1000 'nombre d'itérations
Ldeg = 40000 / 360
pas = (Latmax - Latmin) / n
For n = 1 To n
lat = Latmin + (n - 1 / 2) * pas
s = s + Ldeg * pas * (Longmax - Longmin) * Ldeg * Cos(lat * Application.Pi / 180)
Next
MsgBox "Surface " & Format(s, "0.00 \k\m2") & " ou " & Format(s * 100, "0 \h\a")
End Sub
=(2*(ATAN(COS((ABS(S-N)-COS(N*PI()/180)*ABS(O-E))*PI()/360)/COS((ABS(S-N)+COS(N*PI()/180)*ABS(O-E))*PI()/360))
+ATAN(COS((ABS(S-N)-COS(S*PI()/180)*ABS(O-E))*PI()/360)/COS((ABS(S-N)+COS(S*PI()/180)*ABS(O-E))*PI()/360)))
-PI())*RT^2