%pageRank

M=[1/2, 1/2, 0;  
   1/2, 0,   0; 
   0,   1/2, 1];
N=3;
R=ones(3,3)/N
r=[1/N;1/N;1/N];
interval=1:30;
beta=0.8;
for i=interval
    x(i)=r(1);
    y(i)=r(2);
    z(i)=r(3);
    r=(beta*M+(1-beta)*R)*r
    
end;
plot(interval, x, '+-', interval,y,'-', interval,z,'*-');
xlabel('iterations');
ylabel('ranks');
legend('y','a','m');
