The nonlinear Schrödinger equation is used to show how numerical methods can be used to solve mathematical problems present in nonlinear analysis. The Lanzos-Chevbychev Pseudospectral method is shown to be effective, flexible, and economical to meet various demands in practical applications of mathematical simulations using nonlinear differential equations. The electromagnetic wave propagation through an inhomogeneous, anisotropic, and complex space is used as an example to show how successful mathematical modeling could be used to explain the complex phenomenon of astronomical redshift that is the central issue in the widely debated Hubble tension.