Abstract. A general methodology for the stability analysis of discrete approximations of nonstationary PDEs is applied to solve the Kuramoto-Tsuzuki equation, including also the Schrödinger problem. Stability regions are constructed for the explicit, backward and symmetrical Euler schemes. The obtained results are applied to solve the Kuramoto-Tsuzuki problem with a non-local integral boundary condition. Results of computational experiments are provided.