We present a novel representation for generalized hypergeometric functions of type p+1 F p which is in fact defined by an infinite series in nonnegative integer powers of its argument. We first construct a first order vector differential equation such that the unknown vector's coefficient is the sum of a constant matrix and a matrix premultiplied by the reciprocal of the independent variable whereas its first order derivative has unit matrix coefficient. An infinite process of factor extractions and power annihilations is employed yielding finally a vector differential equation that can be easily and analytically solved. Truncation of this scheme can be used to get approximations to hypergeometric functions of type p+1 F p . These functions have regular singularities at 0 and 1 values of the independent variable together with another regular singularity at infinity. Hence the factors are chosen to reflect the expected behavior of the functions at the singular point in a descending contribution order. Factorization is realized also for regular points. A simple, yet meaningful, implementation seems to give quite promising results.