Matrix reduction and approximation to principal axes