​​​​θ=Pi/2;a=0.3;r=0.2;​​F[r_,c_]:=
aSin
θ
2

3
Cosh[r]
Coth[r]PolyLog-
1
2
,
2
Tanh[r]
+
1-
2
a
-
2
c
1+
2
a
2
c
Sin
θ
2
+cCos
θ
2
+
1-
2
a
-
2
c
-
1-
2
a
-
2
c
1+
2
a
2
c
Cos
θ
2

3
Cosh[r]
Coth[r]PolyLog-
1
2
,
2
Tanh[r]
;​​​​
M
r
=Maximize[F[r,c],c]​​​​​​​​
Maximize0.0417288+
c
2
+
0.91-
2
c
1+
0.09
2
c
2
+0.139096
0.91-
2
c
-
0.91-
2
c
1+
0.09
2
c
,c