function?LST（a）
%a是坐標矩陣
[mn]=size（a）;
A=zeros（size（a1））;
for?i=1:m
????for?j=1:i
??????A（ij）=norm（a（j:）-a（i:））;%坐標生成鄰接矩陣
????end
end
A=A+A‘;
A（A==0）=inf;
E=[];
p=1;%設置生成樹的起始頂點
b=2:length（A）;?%設置生成樹以外頂點
while?size（E2）~=length（A）-1
	t=A（pb）;%表示在已經形成的樹每個點到下一個沒標記的點的距離，以后的程序要比較以尋求最小
	t=t（:）;%取出與p關聯的所有邊
	d=min（t）;%取上述邊中的最小邊
	[jk]=f