function Length=GetLength(iEnd,Nodes) node_j=iEnd(1); node_k=iEnd(2); x_j=Nodes{node_j}.x; x_k=Nodes{node_k}.x; y_j=Nodes{node_j}.y; y_k=Nodes{node_k}.y; L=(x_k-x_j)^2+(y_k-y_j)^2; Length=sqrt(L); end