利用fsolve函数求解六点确定三正交平面法向量的问题

该程序为利用fsolve求解六点确定三正交面法向量的程序，已知量为六点坐标，其中点1，点2，点3位于同一平面上，点4，点5位于另一平面上，点6位于第三平面上，三平面两两垂直，待求量为三个平面的法向量 

主程序为：

fun=@liudiansanmian
x0=[0,0,1,0,1,0,1,0,0];
x=fsolve(fun,x0);
vector1=[x(1) x(2) x(3)];%平面法向量
vector2=[x(4) x(5) x(6)];%平面法向量
vector3=[x(7) x(8) x(9)];%平面法向量
norm1=norm(vector1);%向量的模
norm2=norm(vector2);%向量的模
norm3=norm(vector3);%向量的模
normalized_vector1=[x(1)/norm1 x(2)/norm1 x(3)/norm1];%平面法向量归一化
normalized_vector2=[x(4)/norm2 x(5)/norm2 x(6)/norm2];%平面法向量归一化
normalized_vector3=[x(7)/norm3 x(8)/norm3 x(9)/norm3];%平面法向量归一化

函数程序为： 

function f=liudiansanmian(x)
f1=x(1)*(-124)+x(2)*51+x(3)*265+1;
f2=x(1)*(-39)+x(2)*24+x(3)*295+1;
f3=x(1)*(-93)+x(2)*(-56)+x(3)*260+1;
f4=x(4)*(-110)+x(5)*(-68)+x(6)*224+1;
f5=x(4)*(-72)+x(5)*(-55)+x(6)*136+1;
f6=x(1)*x(4)+x(2)*x(5)+x(3)*x(6);
f7=x(7)*(-135)+x(8)*10+x(9)*186;
f8=x(1)*x(7)+x(2)*x(8)+x(3)*x(9);
f9=x(4)*x(7)+x(5)*x(8)+x(6)*x(9);
f=[f1;f2;f3;f4;f5;f6;f7;f8;f9];

计算结果为：

normalized_vector1 =

    0.3712    0.1504   -0.9163

normalized_vector2 =

    0.4921    0.8050    0.3314

normalized_vector3 =

   -0.5281   -0.7775   -0.3415