Calculation of distance on the ellipsoid

Use for small distances. Error in 10 km distance, between latitudes 50S and 50N, can be up to 2 mm.

Help section
Input: Lon1, Lat1, Lon2, Lat2 (decimal degrees)
Output: dS-lon (m), dS-lat (m), dS (m), Azimuth (deg)