I am attempting to calculate the bearing from an object with geodesic coords given by latitude/longitude and a heading to another object at different latitude/longitude coordinates.
A rough sketch of the problem is
where the triangle represents object one, the arrow represents object one's heading, and the square represents object two.
I am aware this can be done using the haversine formula, as in , but how would one do this while simultaneously taking spheroid earth into consideration.