We present an analytical waveform family describing gravitational waves (GWs) from the inspiral, merger, and ringdown of nonspinning black-hole binaries including the effect of several nonquadrupole modes [ðl ¼ 2; m ¼ AE1Þ; ðl ¼ 3; m ¼ AE3Þ; ðl ¼ 4; m ¼ AE4Þ apart from ðl ¼ 2; m ¼ AE2Þ]. We first construct spin-weighted spherical harmonics modes of hybrid waveforms by matching numerical-relativity simulations (with mass ratio 1-10) describing the late inspiral, merger, and ringdown of the binary with post-Newtonian/effective-one-body waveforms describing the early inspiral. An analytical waveform family is constructed in frequency domain by modeling the Fourier transform of the hybrid waveforms making use of analytical functions inspired by perturbative calculations. The resulting highly accurate, ready-to-use waveforms are highly faithful (unfaithfulness ≃10 −4 -10 −2 ) for observation of GWs from nonspinning black-hole binaries and are extremely inexpensive to generate.