Вычислить угол образованный медианой и биссектрисой треугольника

Function Ugol(a, b, c, Optional grad)
 
'вычисляет угол между биссектрисой и медианой треугольника АВС, опущенными из вершины А
'a,b,c - координаты вершин, одномерные массивы с одинаковыми верхними и нижними индексами
'grad - необязательный параметр. Если True, то значение угла в градусах, иначе в радианах.
'код для VBA Excel
 
Const RAD2GRAD = 57.2957795130823 'градусов в радиане, 180/atn(1)/4
Dim lb&, ub&, i&, labl#, lacl#
  lb = LBound(a): ub = UBound(a)
  ReDim ab#(lb To ub), ac#(lb To ub), bi#(lb To ub), md#(lb To ub)
  For i = lb To ub
    ab(i) = b(i) - a(i): ac(i) = c(i) - a(i)
  Next
  With WorksheetFunction
    labl = Sqr(.SumSq(ab))
    lacl = Sqr(.SumSq(ac))
    For i = lb To ub
      bi(i) = ab(i) / labl + ac(i) / lacl: md(i) = (ab(i) + ac(i)) / 2
    Next
    Ugol = .Acos(.SumProduct(bi, md) / Sqr(.SumSq(bi) * .SumSq(md)))
  End With
  If grad Then Ugol = Ugol * RAD2GRAD
End Function
 
Sub test()
Dim a(), b(), c(), u
  a = Array(1, 0)
  b = Array(3, 0)
  c = Array(3, -5)
  u = Ugol(a, b, c, True)
  Debug.Print u
End Sub