“…When α = 1, because the analytic solutions still remain unknown and the Newton's procedure relies heavily on its initial values, we instead employ the trust-regiondogleg algorithm built into Matlab to improve the convergence of iteration. At first, taking τ = 2.0 × 10 −3 , M = 100, β = 2, and Ω = [−10, 10], the mean square errors at t = 0.1 with the initial condition (38) for various α are reported in Table 5, where the solutions computed by using the coefficients (13) on a very fine time-space lattice, i.e., τ = 2.5 × 10 −4 , M = 400, are adopted as reference solutions (α = 1). As seen from 20,20], we display the evolution of the amplitude of the mobile soliton created by (23) for α = 0.98 and 1.0 in Fig.…”