The production of , and mesons decaying into the dimuon final state is studied with the LHCb detector using a data sample corresponding to an integrated luminosity of collected in proton-proton collisions at a centre-of-mass energy of TeV. The differential production cross-sections times dimuon branching fractions are measured as functions of the transverse momentum and rapidity, over the ranges GeV/ and . The total cross-sections in this kinematic region, assuming unpolarised production, are measured to be sigma (pp -> Upsilon(1S)X) x B(Upsilon(1S) -> mu(+)mu(-)) = 1.111 +/- 0.043 +/- 0.044 nb, sigma (pp -> Upsilon(2S)X) x B(Upsilon(2S) -> mu(+)mu(-)) = 0.264 +/- 0.023 +/- 0.011 nb, sigma (pp -> Upsilon(3S)X) x B(Upsilon(3S) -> mu(+)mu(-))s = 0.159 +/- 0.020 +/- 0.007 nb, where the first uncertainty is statistical and the second systematic.