Magneto-responsive slender bodies are used in a range of promising applications, such as artificial cilia, magnetic fiber networks and cellular actuators. To accurately describe the magneto-elastic deformations, both the demagnetization field as well as the resulting magnetic loads on the body should be properly accounted for. The calculation of the demagnetization field for a general sample shape is very challenging, which has hampered the experimental characterization of the intrinsic magnetic and magnetoelastic properties. Here, a methodology is developed to accurately calculate the demagnetization field for slender bodies,i.e, for long beams having a rectangular or circular cross-section. We propose two different expressions for the mag...