Secondary fluorescence effects are important sources of characteristic X-ray emissions, especially for materials with complicated geometries. Currently, three approaches are used to calculate fluorescence X-ray intensities. One is using Monte Carlo simulations, which are accurate but have drawbacks such as long computation times. The second one is to use analytical models, which are computationally efficient, but limited to specific geometries. The last approach is a hybrid model, which combines Monte Carlo simulations and analytical calculations. In this article, a program is developed by combining Monte Carlo simulations for X-ray depth distributions and an analytical model to calculate the secondary fluorescence. The X-ray depth distribution curves of both the characteristic and bremsstrahlung X-rays obtained from Monte Carlo program MC X-ray allow us to quickly calculate the total fluorescence X-ray intensities. The fluorescence correction program can be applied to both bulk and multilayer materials. Examples for both cases are shown. Simulated results of our program are compared with both experimental data from the literature and simulation data from PENEPMA and DTSA-II. The practical application of the hybrid model is presented by comparing with the complete Monte Carlo program.