metodo de la potencia

function [vector, lambda] = powerade(A, zv, iter)
    // zv es un vector inicial, iter cantidad de iteraciones
    // metodo de la potencia, aproxima el autovalor de mayor
    // valor absoluto, y su autovector asociado
    n = size(A, 1);
    for rep=1:iter
        wn = A*zv;
        zn = wn / norm(wn, 'inf');
        for k=1:n
            if(abs(zn(k)) > 1e-12)
                lambda = wn(k) / zv(k);
                break;
            end
        end
        zv = zn;
    end
    vector = zn;
endfunction