Here is a set of difference equations for y[n]
62022240 + 545995032 n + 2056791388 n^2 + 4333244560 n^3 + 5587600700 n^4 + 4517982000 n^5 + 2238010000 n^6 + 621200000 n^7 + 74000000 n^8 + (-19027008 - 158454120 n - 566231672 n^2 - 1135130960 n^3 - 1397526400 n^4 - 1082880000 n^5 - 516080000 n^6 - 138400000 n^7 - 16000000 n^8) y[n] + 3 (2 + 5 n) (3 + 5 n)^2 (4 + 5 n) (6 + 5 n) (5 + 6 n) (7 + 6 n) (31 + 20 n) y[1 + n] == 0, y[1] == 158/31
or equivalently
4 (3 + 2 n) (5168520 + n (42053906 + 5 n (28672669 + 5 n (10621126 + 5 n (2308917 + 10 n (147271 + 20 n (2551 + 370 n))))))) + (-8 (1 + n) (3 + 2 n) (8 + 5 n) (7 + 10 n) (9 + 10 n) (11 + 10 n) (13 + 10 n) (11 + 20 n)) y[n] + 3 (2 + 5 n) (3 + 5 n)^2 (4 + 5 n) (6 + 5 n) (5 + 6 n) (7 + 6 n) (31 + 20 n) y[1 + n] == 0, y[1] == 158/31
I would appreciate it if a reader could provide the Maple solution to these equations--as I have been analyzing it with Mathematica, and obtaining a "hypergeometric-rich" solution of large complexity.
Maple doesn't return any hypergeometrics. Its solution has a lot of Gamma functions and a finite summation: