Numerical simulations of self-focusing of ultrafast laser pulses