The simulations indicated that DIRA is capable of determining elemental composition of tissues, which are needed in brachytherapy with low energy (< 50 keV) photons and proton therapy. The algorithm provided quantitative monoenergetic images with beam hardening artifacts removed. Its convergence was fast, image sharpness expressed via the modulation transfer function was maintained, and image noise did not increase with the number of iterations.