“…55]. Computing Liouvillian solutions is already solved in [9,11,2,8] where this problem is reduced to computing hypergeometric solutions of some other operatorL (or system, in [8]) which has n times more finite singularities (defined in [10]) than the original operator L. The problem is that computing a hypergeometric solution is done with a combinatorial algorithm [12,7] where the number of combinations depends exponentially on the number of singularities. We give a more direct approach, based on Theorems 2 and 3, that avoids introducingL and the corresponding increase in the number of singularities.…”