The XMM-LSS, XMM-CDFS, and XMM-COSMOS are three surveys with complementary properties in terms of luminosity and redshift coverage. We used these three surveys to derive Bayesian estimates of the unabsorbed luminosity function (LF) of AGN in the 2-10keV band. The LF estimates are presented as a set of samples from the posterior probability distribution of the LF parameters. The LF is parameterised as a double power-law, with either the luminosity and density evolution (LADE) model, or the luminosity-dependent density evolution (LDDE) model. The double power-law is described by Eq.(10) in the paper. The LADE and LDDE models are described by Eqs.(11-14) and Eqs.(15-17), respectively. A Fortran 2008 implementation of these models can be found in file src2/lumf_funcs.f90 of the LFTools package, in the classes doublepowerlaw, ladevol, and lddevol (see the paper). <P />(8 data files). <P />