# Constraining the cosmological equation of state at extremely high redshift

We construct a fully Bayesian methodology tuned to the inference of the dark energy equation of state as a function of redshift in a completely non-parametric way, without imposing any a priori w(z) functional form. As a check of the method, we test the complete scheme through the extensive use of synthetic data sets produced from a known cosmological model, and having the same relative errors and redshift distribution as the real data. This allows to confirm the accurate w(z) reconstruction capabilities of the method proposed. Through an updated form of the L_{X}-T_{a} correlation, we are able to use a large GRB sample from the Swift satellite as accurate tracers of the expansion history of the Universe out to z=9. Within the assumption of a GR universe with an arbitrary dark energy term, the maximum likelihood solution combining supernovae type Ia, our GRB sample and BAO observations is extremely close to a constant w=-1 solution out to z=9, although variations of +- 0.3 can not be ruled out. If one imposes the restriction of a constant w, we obtain w=-1^{+0.06}_{-0.04}, highly consistent with a cosmological constant, firmly determined out to z=9. As additional constraints in the low redshift range, the method yields highly accurate maximum likelihood inferences for the present day parameters :H_{0}=69.5, km s^{-1} Mpc^{-1} and Omega_Lambda_0 = 0.72.