十四面体matlab开发程序-Kelvin模型

function Kelvin（dxdydz）
sq=sqrt（2）
face01=[0 sq -2*sq;sq 0 -2*sq;0 -sq -2*sq;-sq 0 -2*sq]
face02=[0 sq 2*sq;sq 0 2*sq;0 -sq 2*sq;-sq 0 2*sq];
face03=[0 2*sq -sq;sq 2*sq 0;0 2*sq sq;-sq 2*sq 0];
face04=[0 -2*sq -sq;-sq -2*sq 0;0 -2*sq sq;sq -2*sq 0];
face05=[2*sq 0 -sq;2*sq -sq 0;2*sq 0 sq;2*sq sq 0];
face06=[-2*sq 0 -sq;-2*sq sq 0;-2*sq 0 sq;-2*sq -sq 0];
face07=[0 sq -2*sq;sq 0 -2*sq;2*sq 0 -sq;2*sq sq 0;sq 2*sq 0;0 2*sq -sq];
face08=[0 sq 2*sq;sq 0 2*sq;2*sq 0 sq;2*sq sq 0;sq 2*sq 0;0 2*sq sq];
face09=[sq 0 -2*sq;0 -sq -2*sq;0 -2*sq -sq;sq -2*sq 0;2*sq -sq 0;2*sq 0 -sq];
face10=[sq 0 2*sq;0 -sq 2*sq;0 -2*sq sq;sq -2*sq 0;2*sq -sq 0;2*sq 0 sq];
face11=[0 -sq -2*sq;-sq 0 -2*sq;-2*sq 0 -sq;-2*sq -sq 0;-sq -2*sq 0;0 -2*sq -sq];
face12=[0 -sq 2*sq;-sq 0 2*sq;-2*sq 0 sq;-2*sq -sq 0;-sq -2*sq 0;0 -2*sq sq];
face13=[-sq 0 -2*sq;0 sq -2*sq;0 2*sq -sq;-sq 2*sq 0;-2*sq sq 0;-2*sq 0 -sq];
face14=[-sq 0 2*sq;0 sq 2*sq;0 2*sq sq;-sq 2*sq 0;-2*sq sq 0;-2*sq 0 sq];
face01（:1）=face01（:1）+dx;face01（:2）=face01（:2）+dy;face01（:3）=face01（:3）+dz;
face02（:1）=face02（:1）+dx;face02（:2）=face02（:2）+dy;face02（:3）=face02（:3）+dz;
face03（:1）=face03（:1）+dx;face03（:2）=face03（:2）+dy;face03（:3）=face03（:3）+dz;
face04（:1）=face04（:1）+dx;face04（:2）=face04（:2）+dy;face04（:3）=face04（:3）+dz;
face05（:1）=face05（:1）+dx;face05（:2）=face05（:2）+dy;face05（:3）=face05（:3）+dz;
face06（:1）=face06（:1）+dx;face06（:2）=fac