'Converts a singular initial condition to ordinary if possible.'

def _singularics_to_ord(self):

    a = list(self.y0)[0]
    b = self.y0[a]
    if ((len(self.y0) == 1) and (a == int(a)) and (a > 0)):
        y0 = []
        a = int(a)
        for i in range(a):
            y0.append(S(0))
        y0 += [(j * factorial((a + i))) for (i, j) in enumerate(b)]
        return HolonomicFunction(self.annihilator, self.x, self.x0, y0)
