“…Relevant numerical methods are based on the direct and inverse fast Fourier transforms (FFTs) for the imaginary-and real-time propagation and the Crank-Nicolson discretization [88][89][90][91]. Apart from the numerical approaches, the variational approximation (VA) is a powerful method for the semi-analytical study of the dynamics of BEC, such as investigating oscillation modes, frequency shift, stability, Faraday patterns, coupled modes, parametric and geometric resonances, temporal modulation of the s-wave scattering length, and localization in OL potentials [22,59,60,73,[92][93][94][95][96][97][98][99][100]. In particular, the VA was used to predict stable solitons with the strongly anisotropic structure in the plane including the direction of the polarizing magnetic field, which is combined with a trapping HO potential acting in the direction perpendicular to the plane [73,75].…”