support more special functions

Author: zhujch1

ext/TaylorDiffSFExt.jl

@@ -8,7 +8,8 @@ using ChainRules, ChainRulesCore
 
 dummy = (NoTangent(), 1)
 @variables z
-for func in (erf,)
+# logerfc, logerfcx, erfinv, gamma, digamma, trigamma
+for func in (erf, erfc, erfcinv, erfcx, erfi)
     F = typeof(func)
     # base case
     @eval function (op::$F)(t::TaylorScalar{T, 2}) where {T}