Five-body Moshinsky brackets that relate harmonic oscillator wavefunctions in two different sets of Jacobi coordinates make it straightforward to calculate some matrix elements in the variational calculations of five-body systems. The analytical expression of these transformation coefficients and the computer code written in the Mathematica language are presented here for accurate calculations. Program summary Program title: FBMB Catalogue identifier: AEZT_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEZT_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 15635 No. of bytes in distributed program, including test data, etc.: 1498331 Distribution format: tar.gz Programming language: Mathematica. Computer: All computers running Mathematica. Operating system: Operating systems running Mathematica. RAM: 20265384 bytes Classification: 17.17. Nature of problem: The calculation of five-body Moshinsky brackets Solution method: The method is based on the explicit expressions of five-body Moshinsky brackets presented in [1] Unusual features: The possibility of calculating five-body Moshinsky brackets Running time: Less than one second for a single five-body Moshinsky bracket References: [1] Shuyuan Xiao, Xueli Mu, Zhixuan Deng, Hong Chen, Journal of Mathematical Physics 56, 042102 (2015) (C) 2016 Elsevier B.V. All rights reserved.