Abstract
Abstract Central to much of science, engineering, and society today is the building of mathematical models to represent complex processes. Recently, the non-relativistic limit of nonlinear Klein-Gordon equations in de Sitter spacetime has been used to derive Schrödinger equations with weighted nonlinear terms In this paper, numerical simulations are constructed to clarify the behavior of the solution in both one- and two-dimensions. These simulations are constructed based on the Crank-Nicolson scheme in the time direction and the Galerkin finite element in the spatial direction. The nonlinear system of algebraic equations resulting from the constructed scheme is solved using Newton’s method. It is also demonstrated that the numerical approximation converges to the exact one.