Few-body systems are the simplest systems where predictions of realistic models can be accurately computed and the results compared to complete experimental measurements. As an example, I will discuss the construction of a three-nucleon interaction, called the Tucson-Melbourne interaction in the literature, which is based on the interchange of virtual pions. This interaction is subjected to the constraints of the chiral symmetry of the pion-nucleon (strong) interaction and of pion-nucleon scattering data. This three-body force is meant to be used in the non-relativistic many-nucleon Schroedinger equation. I will describe contemporary mathematical techniques to solve the three and four-nucleon bound states (nuclei) with realistic two- and three-body forces and conclude with recent bound state (and some continuum) results.