用matlab写段代码，已知三维空间三个点坐标，求其构成的法向量，不使用叉积的方法

在MATLAB中，如果你想要计算通过三维空间中三个不共线点P1(x1, y1, z1)，P2(x2, y2, z2)和P3(x3, y3, z3)的平面的法向量，可以利用向量减法和向量内积的性质来求解。这里假设这三个点不是同一直线上的一点。

首先，创建两个向量，分别代表两个差分向量：
v1 = P2 - P1;
v2 = P3 - P1;

然后，你可以通过这两个向量的向量内积（也称为点积或标量积）来得到法向量n。因为向量n是垂直于这两条直线的，所以它将与这两个向量都垂直。法向量n的公式是：
n = v1 × v2; // 这里省略了实际的交叉乘积运算

但是题目要求不使用叉积直接计算，我们可以利用向量的坐标表示来实现。由于向量n = a × b 等价于 n_x = (a_y * b_z - a_z * b_y), n_y = (a_z * b_x - a_x * b_z), n_z = (a_x * b_y - a_y * b_x)，你可以分别计算出每个分量：

[x1, y1, z1] = [p1(1), p1(2), p1(3)];
[x2, y2, z2] = [p2(1), p2(2), p2(3)];
[x3, y3, z3] = [p3(1), p3(2), p3(3)];

dx = x2 - x1;
dy = y2 - y1;
dz = z2 - z1;

dxa = y2 - y3;
dyb = z2 - z3;
dzb = x2 - x3;

n_x = dy * dzb - dz * dyb;
n_y = dz * dx - dx * dz;
n_z = dx * dyb - dy * dxa;

% 计算后的法向量就是[n_x, n_y, n_z]
normal_vector = [n_x, n_y, n_z];