matlab生成一个三角形平面的法向量

%求三角形法向量
function [x,y,z]=fa_vector(nodes,x1,x2,x3)
%三个点的坐标
p1=[nodes(x1,2),nodes(x1,3),nodes(x1,4)];
p2=[nodes(x2,2),nodes(x2,3),nodes(x1,4)];
p3=[nodes(x3,2),nodes(x3,3),nodes(x3,4)];
%两个边向量
a=p2-p1;
b=p3-p1;
%求法向量
c=cross(b,a);
%归一化
 norm = sqrt(c(1,1)^2+c(1,2)^2+c(1,3)^2);
 x=c(1,1)/norm;
 y=c(1,2)/norm;
 z=c(1,3)/norm;
end