In terms of Hurwitz zeta function, Dirichlet beta function is given by $$\beta(s)=\frac1{4^s}\left(\zeta(s,\frac14)+\zeta(s,\frac34)\right).$$ Following the links, by means of analytic continuation of Hurwitz zeta function, I found that $$\beta(s)=\frac{\Gamma(1-s)}{4\pi i}\int_\gamma\frac{z^{s-1}}{\cosh z} dz$$ where $\gamma$ is the contour
And by using the functional equation of Dirichlet beta function, I found $\beta(2k+1)$, as shown by $(7)$, here.
I don't know much about these topics. I am following the links and formulas. Is my analytic continuation correct? Can it be useful for computation of $\beta(2k)$ which are related to Euler's zigzag numbers?
I would appreciate any help/reference in these topics. Thanks in advance.

The contour integral representation is correct, and it's very similar to the representation $$\zeta(s) = \frac{\Gamma(1-s)}{2 \pi i } \int_{\gamma} \frac{z^{s-1}}{e^{\color{red}{-}z}-1} \, \mathrm dz. $$
Using the principal branch of $z^{s}$ and assuming that $\Re(s) >0$, integrate the function $$f(z) = \frac{z^{s-1}}{\cosh (z)} $$ counterclockwise on a contour $(\mathcal{C})$ that consists of a ray just below the branch cut on the negative real axis, a small semicircle about the origin, and a ray just above the branch cut on the negative real axis.
We get $$\begin{align} \int_{\mathcal{C}} \frac{z^{s-1}}{\cosh(z)} \, \mathrm dz &= \int_{\infty}^{\epsilon} \frac{(te^{-i \pi})^{s-1}}{\cosh(te^{- i \pi})} \, e^{- i \pi} \, \mathrm dt + \int_{|r|= \epsilon} f(z) \, \mathrm dz + \int_{\epsilon}^{\infty} \frac{(te^{i \pi})^{s-1}}{\cosh (te^{ i \pi})} \, e^{i \pi }\mathrm \, \mathrm dt \\ &=-e^{- i \pi s}\int^{\infty}_{\epsilon} \frac{t^{s-1}}{\cosh(t)} \, \mathrm dt + \int_{|r|= \epsilon} f(z) \, \mathrm dz + e^{i \pi s} \int_{\epsilon}^{\infty} \frac{t^{s-1}}{\cosh (t)} \, \mathrm dt. \end{align}$$
Since $\Re(z) >0$, the second integral on the right vanishes as $\epsilon \to 0$, and we get
$$ \begin{align} \lim_{\epsilon \to 0}\int_{\mathcal{C}} \frac{z^{s-1}}{\cosh(z)} \, \mathrm dz &= \left(-e^{- i \pi s} + e^{i \pi s} \right) \int_{0}^{\infty} \frac{t^{s-1}}{\cosh (t)} \, \mathrm dt \\ &= 2i \sin \left(\pi s \right) \, 2 \, \Gamma(s) \, \beta(s) \\ &= \frac{4 \pi i \, \beta(s)}{\Gamma(1-s)} . \end{align}$$
However, since the magnitude of $\frac{1}{\cosh (z)}$ decays exponentially to zero as $\Re(z) \to - \infty$, the integral $$\int_{\gamma} \frac{z^{s-1}}{\cosh (z)} \, \mathrm dz$$ has the same value as long as it doesn't curve around any of the poles of $\frac{1}{\cosh(z)}$ on the imaginary axis.
(See this related question that I asked a long time ago.)
Therefore, we indeed have $$\beta(s) = \frac{\Gamma(1-s)}{4 \pi i} \int_{\mathcal{\gamma}} \frac{z^{s-1}}{\cosh(z)} \, \mathrm dz$$ for all $s \in \mathbb{C}. $
The function $\Gamma(1-s)$ has simple poles at the positive integers, but the limit of the expression on the right side is finite because the integral vanishes when $s$ is a positive integer.
To derive the functional equation, let's close the contour with a square with vertices at $z= \pm 2N \pi \pm i 2N \pi$, where $N$ is some large positive integer.
The function $\frac{1}{\cosh(z)}$ has simple poles at the imaginary half- integers.
But on the square, $\frac{1}{\cosh (z)}$ is uniformly bounded.
So by the estimation lemma, the integral on the square vanishes as $N \to \infty$ through the positive integers if $\operatorname{Re}(s) <0.$
(A more subtle argument shows that it actually vanishes for $\Re(s) <1$. It has to do with the fact that the magnitude of $\operatorname{sech}(z)$ is going to $0$ exponentially fast as $\Re(z) \to \pm \infty$.)
Therefore, for $\Re(s) <0$, we have $$ \begin{align} \beta(s) &= - \frac{\Gamma(1-s)}{4 \pi i } \, 2 \pi i \left( \sum_{k=0}^{\infty} \operatorname{Res} \left[ f(z), \frac{(2k+1) \pi e^{-i \pi/2}}{2} \right] + \sum_{k=0}^{\infty} \operatorname{Res} \left[ f(z), \frac{(2k+1) \pi e^{i \pi/2}}{2} \right] \right) \\ &= - \frac{\Gamma(1-s)}{2} \left( \sum_{k=0}^{\infty} \frac{\left(\tfrac{(2k+1) \pi e^{-i \pi/2}}{2} \right)^{s-1}}{\sinh \left(\tfrac{(2k+1)\pi e^{- i \pi/2}}{2} \right)} + \sum_{k=0}^{\infty} \frac{\left(\tfrac{(2k+1) \pi e^{i \pi/2}}{2} \right)^{s-1}}{\sinh \left(\tfrac{(2k+1)\pi e^{i \pi/2}}{2} \right)}\right) \\ &= -\frac{\Gamma(1-s)}{2} \left(\frac{\pi}{2} \right)^{s-1} \, \frac{1}{i} \left(-e^{-i \pi(s-1)/2} +e^{i \pi(s-1)/2} \right)\sum_{k=0}^{\infty} \frac{(-1)^{k}}{(2k+1)^{1-s}} \\ &= -\Gamma(1-s) \left(\frac{\pi}{2} \right)^{s-1} \sin \left(\frac{\pi(s-1)}{2} \right) \beta(1-s). \end{align}$$
Replacing $s$ with $1-s$, we have
$$ \begin{align} \beta(1-s) &= - \Gamma(s) \left(\frac{\pi}{2} \right)^{-s} \sin \left(- \frac{\pi s}{2} \right) \beta(s) \\ &= \left(\frac{\pi}{2} \right)^{-s}\sin \left( \frac{\pi s}{2} \right) \Gamma(s) \beta(s). \end{align} $$