We report on the inclusive production cross sections of [Formula: see text], [Formula: see text], [Formula: see text](1S), [Formula: see text](2S) and [Formula: see text](3S), measured at forward rapidity with the ALICE detector in [Formula: see text] collisions at a center-of-mass energy [Formula: see text] TeV. The analysis is based on data collected at the LHC and corresponds to an integrated luminosity of 1.23 pb[Formula: see text]. Quarkonia are reconstructed in the dimuon-decay channel. The differential production cross sections are measured as a function of the transverse momentum [Formula: see text] and rapidity y, over the [Formula: see text] ranges [Formula: see text] GeV/c for [Formula: see text], [Formula: see text] GeV/c for all other resonances, and for [Formula: see text]. The cross sections, integrated over [Formula: see text] and y, and assuming unpolarized quarkonia, are [Formula: see text] [Formula: see text]b, [Formula: see text] [Formula: see text]b, [Formula: see text] nb, [Formula: see text] nb and [Formula: see text] nb, where the first uncertainty is statistical and the second one is systematic. These values agree, within at most [Formula: see text], with measurements performed by the LHCb collaboration in the same rapidity range.