In [7]:

    

function evalue(A)
    n = size(A)[1]
    v = ones(n)
    
    i = 0
    v_next = A*v / norm(A*v)
    while norm(v_next - v) > 1e-5 && i < 10000
        i += 1
        v = v_next
        v_next = A*v / norm(A*v)
    end
    
    if i == 10000
        return false
    end

    v_next
end


    

    
Out[7]:





evalue (generic function with 1 method)

In [10]:

    

a = randn(3, 3)
eig(a)


    

    
Out[10]:





([-0.713190862923451,0.7343674157670204,0.9178607072846022],
3x3 Array{Float64,2}:
 -0.690279  0.947295  -0.0191414
  0.292139  0.245601  -0.182269 
  0.661943  0.205697   0.983062 )

In [11]:

    

println(evalue(a))


    

    




[0.019183593054240104,0.18227737265349772,-0.9830599926637577]