图的最小生成树(普里姆Prim算法),MATLAB

% Prim算法
% 输带权邻接矩阵 a  
a = []
n = size(a,1);
T = []; c = 0; v = 1; sb = 2:n;   % 1是第一个红色点，sb 是白点集
% 构造初始候选边集
for j = 2:n
   b(1,j-1) = 1;
   b(2,j-1) = j;
   b(3,j-1) = a(1,j);
end  
while size(T,2) < n-1
   [Y,i] = min(b(3,:));   % 在候选边集中找最短边
   T(:,size(T,2)+1) = b(:,i);
   c = c+b(3,i);
   v = b(2,i);   % v是新红点
   temp = find(sb==b(2,i));
   sb(temp) = [];
b(:,i) = [];
% 调整候选边集
    for j = 1:length(sb)    % length(sb) 返回 sb 的行数和列数中最大的
       d = a(v,b(2,j));
        if  d < b(3,j)
           b(1,j) = v;
           b(3,j) = d;
        end
    end  
end  
T, c