%¼ÆËã¸Ë³¤º¯Êý
function [BarLength]=BarsLength(x,y,z,ele)
   BarLength=zeros(size(ele,1),1);
   for iEle =1: size(ele,1)
      BarLength(iEle,1)=((x(ele(iEle,3))-x(ele(iEle,2)))^2+(y(ele(iEle,3))-...
      y(ele(iEle,2)))^2+(z(ele(iEle,3))-z(ele(iEle,2)))^2)^0.5;
   end
end