We present and expand the simplest possible quantum cosmological bouncing model already discussed in previous works: the trajectory formulation of quantum mechanics applied to cosmology (through the Wheeler–De Witt equation) in the Friedmann–Lemaître–Robertson–Walker (FLRW) minisuperspace without spatial curvature. The initial conditions that were previously assumed were such that the wave function would not change its functional form but instead provide a dynamics to its parameters. Here, we consider a more general situation, in practice consisting of modified Gaussian wave functions, aiming at obtaining a nonsingular bounce from a contracting phase. Whereas previous works consistently obtain very symmetric bounces, we find that it is possible to produce highly non-symmetric solutions, and even cases for which multiple bounces naturally occur. We also introduce a means of treating the shear in this category of models by quantizing in the Bianchi I minisuperspace.