三维空间任意点绕任意轴旋转

编一段完整的程序求点p(10,20,10)绕直线段AB旋转45°后的坐标。直线段的两个端点坐标分别为（2，1，1）和（18，25，21）