A Maple code is presented for algebraic collective model (ACM) calculations. The ACM is an algebraic version of the Bohr model of the atomic nucleus, in which all required matrix elements are derived by exploiting the model's SU(1,1)×SO(5) dynamical group. This paper reviews the mathematical formulation of the ACM, and serves as a manual for the code.The code enables a wide range of model Hamiltonians to be analysed. This range includes essentially all Hamiltonians that are rational functions of the model's quadrupole momentsqM and are at most quadratic in the corresponding conjugate momentaπN (−2 ≤ M, N ≤ 2). The code makes use of expressions for matrix elements derived elsewhere and newly derived matrix elements of the operators [π ⊗q ⊗π]0 and [π ⊗π]LM . The code is made efficient by use of an analytical expression for the needed SO(5)-reduced matrix elements, and use of SO(5) ⊃ SO(3) Clebsch-Gordan coefficients obtained from precomputed data files provided with the code.The parameter a in the basis {|(a, λ) ν } for L 2 (R + , β 4 dβ) is a useful scale parameter that implicitly defines the U(1) ⊂ SU(1, 1) subgroup (see Section III). The group SO(5) in the chain (5) is the group of linear transformations of the five quadrupole moments {q M } that leave β 2 invariant, and SO(3) is the rotational subgroup that transforms the quadrupole moments as a basis for the 5-dimensional L = 2 irrep. An extra 'missing label' α in the range 1 ≤ α ≤ d vL is needed to distinguish the multiplicity, d vL , of SO(3) irreps of the same angular momentum L in an SO(5) irrep of seniority v (seniority is the SO(5) analogue of angular momentum). This multiplicity is given [18,24] by