  In this talk, we discuss the algorithm for the construction of analytical
coefficients of higher order epsilon expansion of some Horn type hypergeometric
functions of two variables around rational values of parameters.
