“…where the subscripts denote the partial derivatives, x is the scaled horizontal coordinate, y denotes the scaled space coordinate perpendicular to x, t is the scaled time, the complex function u(x, y, t) is the amplitude of a surface wave packet, the real function v(x, y, t) is the velocity potential of the mean flow interacting with the surface wave, the parameter ε = 1 characterizes the focusing case, and the parameter ε = −1 characterizes the defocusing case [24]. When the media are inhomogeneous or the boundaries are nonuniform, variablecoefficient models are able to describe various situations more realistically than their constant-coefficient counterparts [26]. In this paper, we focus our interest on the variablecoefficient Davey-Stewartson (vcDS) equation for ocean waves, ultra-relativistic degenerate dense plasmas, and Bose-Einstein condensates [24]:…”