Numerical solution of the Smoluchowski kinetic equation and asymptotics of the distribution function

We obtain the numerical solution of the Smoluchowski kinetic equation for model kernels U(M1, M2) varies as (M1+M2)lambda and U(M1, M2) varies as (M1M2)lambda 2/, 0< lambda <or=2. We show that the behaviour of the solution for the kernel U varies as (M1+M2)lambda at 0< lambda <or=1 and U varies as (M1M2)lambda/2 at 0< lambda <or=2 becomes self-similar after some time. The shape of the scaling function is analysed; in particular, a simple approximate expression for it at U varies as (M1+M2)lambda is found. An interesting result is obtained for U varies as (M1M2)lambda 2/, 0< lambda <1: the asymptotic behaviour of the scaling function proved to be non-power. We develop the procedure for determining tcr, the critical indices and the exponent of the power-law asymptotics of the Smoluchowski equation solution. The concrete values of these quantities for the model kernels are obtained. The stability conditions of the algorithm for numerical solving are analysed. The possibilities afforded by numerical solution for investigation of the Smoluchowski equation are discussed.