A program is created to compute recursively the Moshinsky brackets. It is very fast and provides highly accurate results. In the case of the double precision computations with a single-processor consumer notebook, the computing time per bracket at any not small oscillator excitations is on the scale of 10 −8 s and the accuracy is very good for the total number of quanta up to 80. The program is easy to handle.