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.