Abstract In the present paper the magnetic flux penetration dynamics of type-II superconductors in the flux creep regime is studied by analytically solving the nonlinear diffusion equation for the magnetic flux induction, assuming that an applied field parallel to the surface of the sample and using a power-law dependence of the differential resistivity on the magnetic field induction. An exact solution of nonlinear diffusion equation for the magnetic induction B(r, t) is obtained by using a well-known self-similar technique. We study the problem in the framework of a macroscopic approach, in which all length scales are larger than the flux-line spacing; thus, the superconductor is considered as a uniform medium.