The two-dimensional finite difference time domain algorithm is used to numerically reconstruct the electron density profile in O-mode ultrashort pulse reflectometry. A Gaussian pulse is employed as the source of a probing electromagnetic wave. The Gaussian pulse duration is chosen in such a manner as to have its frequency spectrum cover the whole range of the plasma frequency. By using a number of numerical band-pass filters, it is possible to compute the time delays of the frequency components of the reflected signal from the plasma. The electron density profile is reconstructed by substituting the time delays into the Abel integral equation. As a result of simulation, the reconstructed electron density profile agrees well with the assumed profile.