matlab编程实现pagerank公式

G=[0 1 0 0 1 0 0 0 0 0 0 0 0 0;
   1 0 1 0 1 0 0 0 0 0 0 0 0 0;
   0 1 0 1 0 0 0 0 0 0 0 0 0 0; 
   0 0 1 0 1 0 1 0 0 0 0 0 0 0;
   1 1 0 1 0 1 0 0 0 0 0 0 0 0;
   0 0 0 0 1 0 0 0 0 0 1 1 0 0;
   0 0 0 1 0 0 0 1 1 0 0 0 0 0;
   0 0 0 0 0 0 1 0 0 0 0 0 0 0;
   0 0 0 0 0 0 1 0 0 1 0 0 0 0;
   0 0 0 0 0 0 0 0 1 0 1 0 0 1;
   0 0 0 0 0 1 0 0 0 1 0 0 0 0;
   0 0 0 0 0 1 0 0 0 0 0 0 1 0;
   0 0 0 0 0 0 0 0 0 0 0 1 0 1;
   0 0 0 0 0 0 0 0 0 1 0 0 1 0];
[nn]=size（G）;
p=0.85;
delta=（1-p）/n;
sn=sum（G1）;     % 按列求矩阵 G 各列的列和
D=diag（1./sn）;
A=p*G*D + delta;
% 幂迭代法
x=ones（n1）/n; % 迭代初始向量
z=zeros（n1）;
cnt=0;  % 用于记录迭步数
while max（abs（x-z）） > 0.0001 
   z = x;
   x = A*x;
   cnt=cnt+1;
end
[x1index]=sort（x）;
x1=flipud（x1）;
index=flipud（index）;
% 输出结果
out=[1:n; x1‘; ind