Работает вариант:
alpha:= arctan(((x1*y2) - (y1*x2))/( (x1*x2) + (y1*y2))) ;

Вероятно, это получено путем решения матрицы 2 х 2.