In this paper, we propose a method for numerical modeling of the nuclear matter properties within the framework of relativistic heavy-ion collisions using a holographic equation of state. Model's free parameters are calibrated using lattice results for quark masses approximating physical values and adjusted to match the Regge spectra of $\rho$ mesons. Numerical simulations are performed using the iEBE-MUSIC framework, which incorporates the MUSIC relativistic hydrodynamics solver. We modify the code by implementing a tabulated holographic equation of state, enabling simulations of quark-gluon plasma evolution with dynamically generated initial conditions via the 3D Monte Carlo Glauber Model. Finally, the spectra of produced hadrons are computed using a hybrid iSS+UrQMD approach at the freeze-out stage.