“…Further, high-precision solutions of eigenvalues and eigenfunctions play important roles in research and engineering. For example, the representative dynamic analysis of structures based on high-precision solutions of natural frequency (eigenvalue) or vibration mode (eigenfunction) (Chang and Kim, 2016;Wang et al, 2018) enables the accurate identification of the number, size and location of the damages in cracked structures (Fan and Qiao, 2011;Yang and Yang, 2018). However, solutions that meet user-specified error tolerances are challenging to provide, especially for issues involving the coefficients of variable matrices (Huang et al, 2013), coincident and adjacent approximate eigenvalues (Shen and Shieh, 1999), continuous orders of eigenpairs (Lee, 2000), varying boundary conditions (Jin et al, 2013) and variable cross-sections or curvatures (Raveendranath et al, 2000;Tseng et al, 1997;Zhou and Cheung, 2000).…”