eg.
matlab.
There are libs for major programming languages and math scripts which provide all these methods.
In the QR method - the eigenvectors are the product of the orthogonal transformation in each iteration. Which is what the Olver (example code post #3) paper does.