After a short presentation of the Boltzmann-Fokker-Planck (BFP) equation, which was derived in a previous work, two numerical approaches to solve this equation are investigated-the multigroup method and a diamond scheme applied in a consistent way to space and energy variables. Because of the parabolic nature of the Fokker-Planck operator, it is shown that the standard neutron transport codes cannot solve such an equation. With the one-dimensional time-dependent BFP-1 code, many numerical results have been produced. All deal with the transport of charged particles in dense plasmas because such a problem is very severe from a numerical point of view. Other applications can be imagined since the BFP formalism is quite general.