% matrice du web
c1 = [0 1 0 0 1 0 0 0 0 0];
c2 = [0 0 1 1 0 0 0 0 0 0];
c3 = [0 0 0 1 0 0 0 0 0 1];
c4 = [0 1 0 0 0 0 1 0 0 0];
c5 = [0 0 0 0 0 1 0 0 0 0];
c6 = [0 0 0 0 1 0 1 0 0 0];
c7 = [0 0 0 0 0 0 0 1 0 0];
c8 = [0 0 0 0 0 0 0 0 1 0];
c9 = [0 0 0 0 0 0 1 1 0 0];
c10= [0 0 0 0 0 0 0 0 0 0];

sprintf('la matrice C de connectivite du web') 
C = [c1;c2;c3;c4;c5;c6;c7;c8;c9;c10]

pause

sprintf('la matrice Q ')
epsi=1.E-20;
for j=1:10    
 Q(:,j) = C(:,j)/(epsi+sum(C(:,j)));
end

Q

pause

sprintf('le spectre de la matrice Q ') 

eig(Q)

sprintf('construction de la matrice P ')

e = ones(10,1);
d = zeros(1,10);
for i=1:10
    if (sum(C(:,i))== 0)
       d(i) = 1;
    end
end

P = Q + 0.1*e*d

pause

sprintf('le spectre de la matrice P ') 
eig(P)
%sprintf('le conditionnement de la matrice P ') 
%cond(P-eye(10,10),2)

sprintf('La methode de la puissance pour P') 
x0=e;
epsi= 1.E-4;
x1 = puissance(epsi,P,x0);
sprintf('Le score de pertinence est ')
x1

pause

alfa = 0.85
A = alfa*P + 0.1*(1-alfa)*e*e'

pause

sprintf('le spectre de la matrice A ') 
eig(A)
sprintf('La methode de la puissance pour A') 


x0=e;
x1 = puissance(epsi,A,x0);
sprintf('Le score de pertinence est ')
x1

pause

for i=1:10
 alfa = i/10.;   
 sprintf('le conditionnement de la matrice alfa*P -  I ') 
 cond(eye(10,10)-alfa*P,2)
end