In chapters 13 and 14 we focused on perturbation theories for time-independent and time-dependent perturbations. For time-independent perturbations we obtained first-order and second-order corrections to energy eigenvalues and eigenfunctions. For time-dependent perturbations we obtained the probability for a system to make a transition from an ith initial state to fth final state. We have also discussed how to study the effect of adiabatic and sudden perturbations. In the present chapter we consider the variational method. The variational method can be used to construct an approximate ground state eigenfunction. In some cases the method can be used for finding approximately the energy of first few discrete states of a bound system. The perturbation theory uses a set of unperturbed states and the perturbation Hamiltonian H(1) to construct the eigenfunctions and eigenvalues through highly formalized approximations. In contrast, the variational method is based on a guess for the state-a trial state which is used to determine the corresponding energy. Any well behaved trial state introduced will yield an appropriate energy, although some trial states will be better than other, but whatever the guess, we can always make systematic improvements. We can construct a reasonably accurate eigenfunction often by a physical insight and ingenuity.