We review our calculation method, Gaussian expansion method (GEM), and its applications to various few-body (3-to 5-body) systems such as 1) few-nucleon systems, 2) few-body structure of hypernuclei, 3) clustering structure of light nuclei and unstable nuclei, 4) exotic atoms/molecules, 5) cold atoms, 6) nuclear astrophysics and 7) structure of exotic hadrons. Showing examples in our published papers, we explain i) high accuracy of GEM calculations and its reason, ii) wide applicability of GEM, iii) successful predictions by GEM calculations before measurements. GEM was proposed 30 years ago and has been applied to a variety of subjects. To solve few-body Schrödinger equations accurately, use is made of the Rayleigh-Ritz variational method for bound states, the complex-scaling method for resonant states and the Kohn-type variational principle to S-matrix for scattering states. The total wave function is expanded in terms of few-body Gaussian basis functions spanned over all the sets of rearrangement Jacobi coordinates. Gaussians with ranges in geometric progression work very well both for short-range and long-range behavior of the few-body wave functions. Use of Gaussians with complex ranges gives much more accurate solution especially when the wave function has many oscillations. * Electronic address: hiyama@phys.kyushu-u.ac.jp † Electronic address: mkamimura@riken.jp C. Prediction of spin-orbit splitting in hypernuclei 15 D. Prediction of neutron-rich hypernuclei 16 E. Prediction of hypernuclear states with strangeness S = −2 17 F. Strategy of studying hypernuclei and Y N and Y Y interactions 18