If two points are close enough to each other and different from the poles, you can use a few simple triggers:
float dy = lat2 - lat1;
float dx = cosf(M_PI/180*lat1)*(long2 - long1);
float angle = atan2f(dy, dx);
EDIT: , latN longN - dx dy - , . angle, , . , , 180/M_PI.