This article is concerned with the numerical solution of nonlinear hyperbolic Schr$\ddot{o}$dinger equations (NHSEs) via an efficient Haar wavelet collocation method (HWCM). The time derivative is approximated in the governing equations by the central difference scheme, while the space derivatives are replaced by finite Haar series, which transform it to full algebraic form. The experimental rate of convergence follows the theoretical statements of convergence and the conservation laws of energy and mass are also presented, which strengthens the proposed method to be convergent and conservative. The Haar wavelets based on numerical results for solitary wave shape of $\left|\mathsf{\phi }\right|$ are discussed in detail. The proposed approach provides a fast convergent approximation to the NHSEs. The reliability and efficiency of the method are illustrated by computing the maximum error norm and the experimental rate of convergence for different problems. Comparisons are performed with various existing methods in recent literature and better performance of the proposed method is shown in various tables and figures.

In this study, we implemented the well-known Crank–Nicolson scheme for the numerical solution of Schrödinger equation. The numerical results converge to the exact solution because the Crank–Nicolson scheme is unconditionally stable and accurate. We have compared the results for different parameters with analytical solution, and it is found that the Crank–Nicolson scheme is suitable for the numerical solution of Schrödinger equations. Three different problems are included to verify the accuracy, stability, and capability of the Crank–Nicolson scheme.