The postulates of quantum mechanics associate the time-evolution of a system with its time-dependent Schrödinger equation. We start by examining different solutions to this equation for a particle, represented in terms of a Gaussian wave packet, in different scenarios: free, scattered from a potential energy barrier, or trapped in a potential energy well. In some cases, we encounter stationary solutions, in which the probability density does not change in time (a standing wave). These solutions are identified as eigenfunctions of a system Hamiltonian. The properties of the Hamiltonian as a Hermitian operator are introduced, and particularly, the fact that its proper eigenfunctions can compose an orthonormal set, and that the corresponding eigenvalues are real-valued. Learning that all operators that relate to measurables are Hermitian, and that their eigenvalues relate to the measured values, we conclude that the eigenvalues of the energy operator (Hamiltonian) are the energy levels of the quantum system.