使用普通python计算3D空间中2个线段之间的角度

Question

我在3D空间中有两个由其起点和终点坐标定义的线段。

我需要在没有numpy或其他花哨的lib的python中计算这两个段之间的角度。

我写了以下代码，对我来说似乎是正确的：

def CalcAngle(beg1, end1, beg2, end2):

  vect1 = [None, None, None]
  vect2 = [None, None, None]

  for cc in range(3):
    vect1[cc] = end1[cc] - beg1[cc]
    vect2[cc] = end2[cc] - beg2[cc]

  length1 = math.sqrt(vect1[0]**2+vect1[1]**2+vect1[2]**2)
  length2 = math.sqrt(vect2[0]**2+vect2[1]**2+vect2[2]**2)

  dotProd = 0.0
  for cc in range(3):
    dotProd += (vect1[cc] * vect2[cc])

  angleRad = math.acos( dotProd / (length1 * length2) )
  angleDeg = math.degrees(angleRad)

  return angleRad, angleDeg

# 1st segment
beg1 = [47.737, 31.8, 21.762]
end1 = [43.55, 38.956, 42.209]

# 2nd segment 
beg2 = [35.993, 18.033, 56.988]
end2 = [31.497, 18.783, 68.881]

angleRad, angleDeg = CalcAngle(beg1, end1, beg2, end2)

# angleRad is 0.312973304876
# angleDeg is 17.9320494697

我在这里错过了什么吗？ 有什么情况下我的代码无法正常工作吗？

谢谢。