Introduction au Calcul Scientifique
ESI CS 1402 - Printemps 2010

Analyse de la matrice de pagerank

>	with(LinearAlgebra):

>	M:=map(evalf,Matrix([[0,1/5,0,0,0,0,0,0],[0,0,1/3,0,0,0,0],[0,1/5,0,1/4,0,0,0,0],[0,1/5,0,0,0,0,0,0],[0,1/5,0,1/4,0,0,0,0],[1/2,1/5,1/3,1/4,1/2,0,0,1],[1/2,0,0,0,0,1,0,0],[0,0,1/3,1/4,1/2,0,0,0]]));

>	r:=map(evalf,<<1/8,1/8,1/8,1/8,1/8,1/8,1/8,1/8>>);

>	M . r / MatrixNorm(M.r);

>	for i from 1 to 1000 do r := M . r / MatrixNorm(M.r); od:

>

r;