Numerical solution of a fully nonlinear one dimensional in space and three dimensional in velocity space electron kinetic equation is presented. Direct Simulation Monte Carlo (DSMC) method used for the nonlinear Landau-Fokker-Planck (LFP) collision operator is combined with Particle-in-Cell (PiC) simulations. An assumption of a self-consistent ambipolar electric field is used. The illustrative simulation results for the relaxation of the initial temperature perturbation are compared with the antecedent analytical and numerical results.