Compute the coordinates of the end-point of a displacement on a sphere. lat e lon are the coordinates of the starting point, range is the covered distance of the displacements along a great circle and azimuth is the direction of the displacement relative to the North.
This function can also be used to define a spherical coordinate system with rotated poles.
[lato, lono] = reckon(lat, lon, range, azimuth)
[lato, lono] = reckon(lat, lon, range, azimuth, units)
The units of all input and output parameters can be either degrees (default) or radians.