Example

> B = rand(3,3)
 B =
     0.57      0.831      0.462  
    0.683     0.0735      0.435  
    0.396      0.464       0.63  
> C = symm(B)
 C =
     0.57      0.757      0.429  
    0.757     0.0735       0.45  
    0.429       0.45       0.63  
> eigenS = eigs(C)
 eigenS =
   val          vec          
> C* eigenS.vec[;1] - eigenS.val[1] * eigenS.vec[;1]
 2.78e-16  
 5.55e-17  
 2.22e-16  
> C* eigenS.vec[;2] - eigenS.val[2] * eigenS.vec[;2]
-1.39e-16  
-6.25e-17  
-1.67e-16  
> C* eigenS.vec[;3] - eigenS.val[3] * eigenS.vec[;3]
-4.44e-16  
 2.22e-16  
        0