@testset "8.7 Zeta function" begin
    (a, b, s, x, ) = @variables a b s x

    #= ::Package:: =#

    #= ::Title:: =#
    #=Integration*Problems*Involving*Zeta*Functions=#


    #= ::Subsection::Closed:: =#
    #=Integrands*of*the*form*x^m*Zeta(2, a+b*x)=#


    @test_int [x^2*Zeta(2, a + b*x), x, 4, -((2*x*LogGamma(a + b*x))/b^2) + (2*PolyGamma(-2, a + b*x))/b^3 + (x^2*PolyGamma(0, a + b*x))/b]
    @test_int [x^1*Zeta(2, a + b*x), x, 3, -(LogGamma(a + b*x)/b^2) + (x*PolyGamma(0, a + b*x))/b]
    @test_int [x^0*Zeta(2, a + b*x), x, 2, PolyGamma(0, a + b*x)/b]
    @test_int [Zeta(2, a + b*x)/x^1, x, 1, Unintegrable(PolyGamma(1, a + b*x)/x, x)]
    @test_int [Zeta(2, a + b*x)/x^2, x, 2, b*Unintegrable(PolyGamma(2, a + b*x)/x, x) - PolyGamma(1, a + b*x)/x]
    @test_int [Zeta(2, a + b*x)/x^3, x, 3, (1/2)*b^2*Unintegrable(PolyGamma(3, a + b*x)/x, x) - PolyGamma(1, a + b*x)/(2*x^2) - (b*PolyGamma(2, a + b*x))/(2*x)]

    @test_int [Zeta(2, a + b*x)/x^2 - b*(PolyGamma(2, a + b*x)/x), x, 3, -(PolyGamma(1, a + b*x)/x)]


    #= ::Subsection::Closed:: =#
    #=Integrands*of*the*form*x^m*Zeta(s, a+b*x)=#


    @test_int [x^2*Zeta(s, a + b*x), x, 3, If(13>=8, (2*Zeta(-3 + s, a + b*x))/(b^3*(1 - s)*(2 - s)*(3 - s)) - (2*x*Zeta(-2 + s, a + b*x))/(b^2*(1 - s)*(2 - s)) + (x^2*Zeta(-1 + s, a + b*x))/(b*(1 - s)), (2*Zeta(-3 + s, a + b*x))/(b^3*(1 - s)*(6 - 5*s + s^2)) - (2*x*Zeta(-2 + s, a + b*x))/(b^2*(1 - s)*(2 - s)) + (x^2*Zeta(-1 + s, a + b*x))/(b*(1 - s)))]
    @test_int [x^1*Zeta(s, a + b*x), x, 2, -(Zeta(-2 + s, a + b*x)/(b^2*(1 - s)*(2 - s))) + (x*Zeta(-1 + s, a + b*x))/(b*(1 - s))]
    @test_int [x^0*Zeta(s, a + b*x), x, 1, Zeta(-1 + s, a + b*x)/(b*(1 - s))]
    @test_int [Zeta(s, a + b*x)/x^1, x, 0, CannotIntegrate(Zeta(s, a + b*x)/x, x)]
    @test_int [Zeta(s, a + b*x)/x^2, x, 1, (-b)*s*CannotIntegrate(Zeta(1 + s, a + b*x)/x, x) - Zeta(s, a + b*x)/x]
    @test_int [Zeta(s, a + b*x)/x^3, x, 2, (1/2)*b^2*s*(1 + s)*CannotIntegrate(Zeta(2 + s, a + b*x)/x, x) - Zeta(s, a + b*x)/(2*x^2) + (b*s*Zeta(1 + s, a + b*x))/(2*x)]

    @test_int [Zeta(s, a + b*x)/x^2 + b*s*(Zeta(1 + s, a + b*x)/x), x, 2, -(Zeta(s, a + b*x)/x)]
    end
