A numerical solution of Wheeler-De Witt equation for a quantum cosmological model simulating boson and fermion creation in the early Universe evolution is presented. This solution is based on a Wheeler-DeWitt equation obtained by Krechet, Fil'chenkov, and Shikin, in the framework of quantum geometrodynamics for a Bianchi-I metric.