The present day spectrum of the extragalactic background light ( EBL ) in UV , optical and IR wavelengths is the integral result of multiple astrophysical processes going on throughout the evolution of the Universe .
The relevant processes include star formation , stellar evolution , light absorption and emission by the cosmic dust .
The properties of these processes are known with uncertainties which contribute to the EBL spectrum precision .
In the present paper we develop a numerical model of the EBL spectrum while maintaining the explicit dependence on the astrophysical parameters involved .
We constructed a Markov Chain in the parameter space by using the likelihood function built with the up-to-date upper and lower bounds on the EBL intensity .
The posterior distributions built with the Markov Chain Monte Carlo method are used to determine an allowed range of the individual parameters of the model .
Consequently , the star formation rate multiplication factor is constrained in the range 1.01 < C _ { \mbox { sfr } } < 1.69 at 68 \% C.L .
The method also results in the bounds on the lifetime , radius , dust particle density and opacity of the molecular clouds that have large ambiguity otherwise .
It is shown that there is a reasonable agreement between the model and the intensity bounds while the astrophysical parameters of the best fit model are close to their estimates from literature .