% Palade, Catalin-Iulian, 312CA

function R = Iterative(nume, d, eps)
    fid = fopen('graf1','r');
    N = fgetl(fid);
    N = str2double(N);
    mat_adj = zeros(N);
    for i = 1 : N
        v = fgetl(fid);
        n = str2double(v(1));
        for j = 3:2:length(v)
            vj = str2double(v(j));
            mat_adj(n, vj) = 1;
        endfor
    endfor
    R_old = 1/N * ones(N, 1);
    suma = zeros(1, N);
    for i = 1 : N
        for j = 1 : N
            elem = mat_adj(i,j);
            suma(i) = suma(i) + elem;
        endfor
    endfor
    M = zeros(N);
    for i = 1 : N
        M(i,:) = mat_adj(i, :) ./ suma(i);
    endfor
    M = M'
    R = d*M*R_old + (1-d)/N * ones(N, 1)
    while abs(R - R_old) > eps
        R_old = R;
        R = d*M*R_old + (1-d)/N * ones(N, 1)
    endwhile
endfunction