In the process of converting some R code to Julia, i found the matrix Eigen decomposition produce different result compared to R. the following code is tested in Julia 0.6, and R 3.3.3 through RCall 0.8.1

```
using RCall
a=rand(1:5,5,5)
a=a'*a
```

the R eigenvalues are in descending order, but they are the same as in Julia. however, some of the corresponding eigenvectors are different.