我有三个空间点坐标m(a,b,c),n(d,e,f),p(g,h,i),我想求zxz表示法下mnp平面的欧拉角，怎么写这个Python函数

可以使用以下代码来计算平面mnp的ZXZ欧拉角：

import numpy as np

def zxz_euler_angles(m, n, p):
    # 计算mnp平面法向量
    normal = np.cross(n-m, p-m)
    normal = normal / np.linalg.norm(normal)
    
    # 计算旋转轴的方向
    z_axis = np.array([0, 0, 1])
    axis = np.cross(z_axis, normal)
    axis = axis / np.linalg.norm(axis)
    
    # 计算第一个旋转角
    cos1 = np.dot(normal, z_axis)
    sin1 = np.linalg.norm(axis)
    angle1 = np.arctan2(sin1, cos1)
    
    # 计算第二个旋转角
    cos2 = np.dot(np.dot(np.transpose(axis), np.transpose(normal)), np.transpose(z_axis))
    sin2 = np.dot(np.cross(normal, z_axis), axis)
    angle2 = np.arctan2(sin2, cos2)
    
    # 计算第三个旋转角
    cos3 = np.dot(np.transpose(normal), np.transpose(axis))
    sin3 = np.linalg.norm(normal)
    angle3 = np.arctan2(sin3, cos3)
    
    return (angle1, angle2, angle3)

这个函数使用numpy库来进行向量计算，首先计算平面mnp的法向量，然后计算旋转轴的方向，最后按照ZXZ欧拉角的顺序计算三个旋转角。函数返回一个包含三个角度值的元组。