Numerical simulation of short laser pulse relativistic self-focusing in underdense plasma