An algorithm to compute Liouvillian solutions of prime order linear difference-differential equations