Methods of numerical solution of the Vlasov equation for a charged particle beam are concerned. These methods are based on the method of macroparticles and require a great number of computations. As a result of the investigation, we find optimal combinations of parameters, which allow to increase computational efficiency. Accuracy of the methods was determined by comparing of a numerical solution and a corresponding analytical solution of the Vlasov equation.