Quaternionic eigenvalue problem