Here \(E(z)= H(z)/H_0\) . If we have the functional form of \(E(z)\)  then we can easily integrate Eq. REF  by using any numerical integration method. However, in our case we have only \(H(z)\)  estimates at certain redshifts. In our analysis, to obtain continuous smooth values of \(H(z)\)  we use the Gaussian Process (GP), a well known hyper-parametric regression method [1]}. GP has been widely used in the literature to reconstruct the shapes of physical functions and is very useful for such functional reconstructions due to flexibility and simplicity.
