ERROR: MethodError: no method matching iterate(::Type{HypergeometricFunctions.ℕ₀})

[Li-shiyue]

Is there anyway toI put TaylorN in the first or second argument position of a hypergeometric function？
I've tried a lot of methods but none of them solved it

x, y = set_variables("x y", order=8)
Base.float(t::HomogeneousPolynomial) = HomogeneousPolynomial(float.(t.coeffs))
Base.float(t::TaylorN) = TaylorN(float.(t.coeffs), t.order)
TaylorSeries.TaylorN{T}(a::S) where {T,S} = (R=promote_type(T,S); TaylorN(R(a), get_order()))
Base.eps(::Type{TaylorN{T}}) where {T} = eps(one(T))
Base.eps(a::TaylorN) = eps(constant_term(a))
HypergeometricFunctions.pFq((x,),(2,),x+y)pFq((2,),(4,),x+0.5y)

ERROR: MethodError: no method matching iterate(::Type{HypergeometricFunctions.ℕ₀})
Closest candidates are:
iterate(!Matched::Union{LinRange, StepRangeLen})
@ Base range.jl:880
iterate(!Matched::Union{LinRange, StepRangeLen}, !Matched::Integer)
@ Base range.jl:880
iterate(!Matched::T) where T<:Union{Base.KeySet{<:Any, <:Dict}, Base.ValueIterator{<:Dict}}
@ Base dict.jl:698
...