Feature recurse roots

There are still possible improvements to the code: 
1) Faster upper-Hessenberg eigvals
2) Removing duplicate colleague matrix code
3) More wide ranging testing
but I will make a pull request now because these are 
not too pressing issues.