def _rewrite_saxena_1(fac, po, g, x):
    _, s = _get_coeff_exp(po, x)
    a, b = _get_coeff_exp(g.argument, x)
    period = g.get_period()
    a = _my_principal_branch(a, period)

    # We substitute t = x**b.
    C = fac/(abs(b)*a**((s + 1)/b - 1))
    # Absorb a factor of (at)**((1 + s)/b - 1).

    def tr(l):
        return [a + (1 + s)/b - 1 for a in l]
    return C, meijerg(tr(g.an), tr(g.aother), tr(g.bm), tr(g.bother),
                      a*x)
Rewrite the integral fac*po*g dx, from zero to infinity, as
dataset/ETHPy150Open sympy/sympy/sympy/integrals/meijerint.py/_rewrite_saxena_1