This study is devoted to designing two hybrid computational algorithms to find approximate solutions for a class of singularly perturbed parabolic convection–diffusion–reaction problems with two small parameters. In our approaches, the time discretization is first performed by the well-known Rothe method and Taylor series procedures, which reduce the underlying model problem into a sequence of boundary value problems (BVPs). Hence, a matrix collocation technique based on novel shifted Delannoy functions (SDFs) is employed to solve each BVP at each time step. We show that our proposed hybrid approximate techniques are uniformly convergent in order