Numerical simulations of Eq. (REF ) have been carried out in real space. A straightforward finite-difference scheme in space and an Euler scheme in time have been employed [1]}, using \(\Delta t=0.1\)  and \(\Delta x=1\) . A homogeneous initial condition \(m(i,(j),k)=-1\)  for all \(i,j,k=1,2,...,L\)  has been considered, while the opposite spin value has been considered at the Dirichlet boundary corresponding to \(i=0\) .
