Bound electronic states calculated by expanding their eigenstates in terms of the eigenfunctions of a particle in a box. Bound states of single and multiple symmetric or non-symmetric wells are calculated within the single-band effective mass approximation. The method may be very useful to design narrow band quantum cascade photodetectors to work without applied bias in a photovoltaic mode. The method can be applied to more complex QW potentials.