Non-Singular cosmological solutions are found in a simple model of Einstein gravity coupled non-minimally to a self-interacting real scalar field. The solutions are either full bouncing solutions or partial bounces with a minimal finite size and
a well-behaved future. The characteristics of the solutions are presented and analyzed with an effective potential as a central tool. Bouncing solutions exist for the Higgs-like self interaction potential which is bounded from below, contrary to previous claims in favor of unbounded potentials. Appearance of higher power terms in the potential, gives rise to bouncing solutions even if the self-interaction potential has a symmetric minimum only. The bouncing solutions exist in a finite domain in parameter space with nonzero measure in the space of initial conditions. In other regions of parameter space, solutions with a Big Bang behavior exist, as well as others. While full bounce solutions are unstable under anisotropic or inhomogeneous perturbations for the spatially flat case (k=0), their spatially closed (k>0) counterparts can be stable. | 
