Numerical experiment for focus of rotating and propagating LIB in plasma. I: Quasi-neutral approximation