A definite integral of the hyperbolic secant function in terms of the Polygamma function given by: ​​
π/2
∫
0
2
Sech[μTan[x]]
x=
1
π
PolyGamma1,
1
2
+
μ
π

​​where Re(μ)>0.
In[]:=
Plot3D[
2
Sech[μTan[x]]
,{μ,-3,3},{x,-π/2,π/2},AxesLabelAutomatic]
Out[]=
In[]:=
fun1[μ_]:=NIntegrate[
2
Sech[μTan[x]]
,{x,0,π/2}],
PolyGamma1,
1
2
+
μ
π

π
//N
In[]:=
fun1closed[μ_]:=Integrate[
2
Sech[μTan[x]]
,{x,0,π/2}],
PolyGamma1,
1
2
+
μ
π

π

In[]:=
fun1closed[μ]
Out[]=

π
2
∫
0
2
Sech[μTan[x]]
x,
PolyGamma1,
1
2
+
μ
π

π

In[]:=
fun1[π+]
Out[]=
{0.27626-0.0783859,0.27626-0.0783859}
In[]:=
fun1[π]
Out[]=
{0.297557,0.297557}
In[]:=
fun1closed[π]
Out[]=

π
2
∫
0
2
Sech[πTan[x]]
x,
-4+
2
π
2
π

In[]:=
fun1[π/4]
Out[]=
{0.809105,0.809105,0.809105}
In[]:=
fun1closed[π/4]//FunctionExpand//FullSimplify
Out[]=

π
2
∫
0
2
Sech
1
4
πTan[x]
x,-
8Catalan
π
+π
We can also take a derivative with respect to μ to get Zeta values
In[]:=
D
2
Sech[μTan[x]]
,
PolyGamma1,
1
2
+
μ
π

π
,μ/.{μπ/4}//FunctionExpand//FullSimplify
In[]:=
NIntegrate
2
Sech
1
4
πTan[x]
Tan[x]Tanh
1
4
πTan[x],{x,0,π/2},-π+
56Zeta[3]
2
2
π
//N
Out[]=
{0.268635,0.268635}
In[]:=
Integrate
2
Sech
1
4
πTan[x]
Tan[x]Tanh
1
4
πTan[x],{x,0,π/2},-π+
56Zeta[3]
2
2
π

Out[]=

π
2
∫
0
2
Sech
1
4
πTan[x]
Tan[x]Tanh
1
4
πTan[x]x,-π+
28Zeta[3]
2
π

We can take the infinite sum over μ∈[1,∞) to get a constant value in terms of a definite integral and sum of the polygamma function given by:
In[]:=
∞
∑
μ=1

2
Sech[x]
2
x
+
2
μ
,
PolyGamma1,
1
2
+
μ
π

μπ

In[]:=

∞
∫
0
(-1+πxCoth[πx])
2
Sech[x]
2
2
x
x,
∞
∑
μ=1
PolyGamma1,
1
2
+
μ
π

πμ
//N
Out[]=
{1.30636,1.30636}