Contents:
great_circle_normal
Returns the unit normal(s) to the great circle(s)
n-vector(s) [no unit] of position A and B, decomposed in E.