We study an asymptotic behavior of the return probability for the critical random matrix ensemble in the regime of strong multifractality. The return probability is expected to show critical scaling in the limit of large time or large system size. Using the supersymmetric virial expansion we confirm the scaling law and find analytical expressions for the fractal dimension of the wave functions $d_2$ and the dynamical scaling exponent $mu$. By comparing them we verify the validity of the Chalkers ansatz for dynamical scaling.