I know how to do this on a plane, the trouble I am having is in, as you say, adapting the solution to a sphere.
On the plane I used simple trigonometry with the coordinates of an unknown point C(X,Y) being
X= a cos(alpha) + b sin(alpha)
Y= -a cos(alpha) + b sin (alpha)
with the coordinates of A...