An iterative method for calculating several of the extreme eigensolutions of large real non-symmetric matrices