小弟想复现论文中的一个方程，求解有限元软件输出的节点的刚体位移和转角，公式在图里
x、y、z和xi、yi、zi是变形前后的节点坐标是已知向量，abcθxθyθz是我要求解出来的刚体位移转角
Q这个公式是个解析式，除了已知向量x、y、z和xi、yi、zi，剩下的都是符号，我用循环求和∑数据量太大会卡死
偏微分方程也不好解，请问各位大佬有什么编码思路

将12式带入11式，这个可以手算。

得到方程组后，将数据带入求解方程组。

n取值不要太大，不然计算很慢，内存不够也容易卡死，可以先取小一点的值，慢慢尝试