real analysis - Is $intlimits_0^inftyfrac{sin y}{y^{s+1}}dy=-Gamma(-s)sin(frac{pi s}{2})$ for $operatorname{Re}(s)in (-1,0)$ obvious?

This is a part of computation in Titchmash, Theories of Zeta Functions which I do not find obvious but there is no explanation. I did figure out the computation.

$$\int_0^\infty\frac{\sin(y)}{y^{s+1}}dy=-\Gamma(-s)\sin\left(\frac{\pi s}{2}\right)$$

Q: There is no explanation in the book for this step. Why is this obvious without explanation? My recipe goes as the following. It suffices to restrict to real axis part with $s\in (-1,0)$ region. Now integral is real valued in this region. Here I need $\Gamma(-s)=\frac{\Gamma(-s+1)}{s}$ extension to obtain real valuedness. Consider the integral as the imaginary part of $\int_0^{i\infty} \frac{e^{z}}{i^s z^{s+1}}dz$ where I have already rotated axis by $i$ multiplication. Now to obtain $\Gamma$ function, close contour from $(+\infty,0)$ axis portion and connect to $(0,i\infty)$ portion. Then close the contour by arc. The arc contour contribution is $0$ via exponential suppresion. Then apply residue theorem easily as the whole thing is holomorphic by $s\in (-1,0)$ region. Hence equality follows. This is not $1-2$ line naive computation though not hard. However, it did take me a while to figure out.

Answer

I realised this question has been asked before as you can see here. Anyway I will write down my solution here again. First of all consider Ramanuajan's Master Theorem.

Ramanujan's Master Theorem

Let $f(x)$ be an analytic function with a MacLaurin Expansion of the form

$$f(x)=\sum_{k=0}^{\infty}\frac{\phi(k)}{k!}(-x)^k$$ then the Mellin Transform of this function is given by

$$\int_0^{\infty}x^{p-1}f(x)dx=\Gamma(p)\phi(-p)$$

In order to use this Theorem we may expand the sine function as a series followed by the substitution $y^2=t$ which yields to

$$\begin{align*} \mathfrak{I}=\int_0^{\infty}y^{-s-1}\sin(y)dy&=\int_0^{\infty}y^{-s-1}\sum_{n=0}^\infty (-1)^n \frac{y^{2n+1}}{(2n+1)!}dy\\ &=\frac12\int_0^{\infty}y^{-s-1}\sum_{n=0}^\infty (-1)^n \frac{n!/(2n+1)!}{n!}(-y^2)^n[2ydy]\\ &=\frac12\int_0^\infty t^{-(s+1)/2}\sum_{n=0}^\infty\frac{n!/(2n+1)!}{n!}(-t)^ndt \end{align*}$$

Now we can use Ramanuajan's Master Theorem by setting $p=-\frac{s-1}2$ and $\phi(n)=\frac{n!}{(2n+1)!}=\frac{\Gamma(n+1)}{\Gamma(2(n+1))}$ and so we get

$$\begin{align*} \mathfrak{I}=\frac12\int_0^\infty t^{-(s+1)/2}\sum_{n=0}^\infty\frac{n!/(2n+1)!}{n!}(-t)^ndt&=\frac12\Gamma\left(-\frac{s-1}2\right)\frac{\Gamma\left(1+\frac{s-1}2\right)}{\Gamma\left(2\left(\frac{s-1}2+1\right)\right)}\\ &=\frac1{2\Gamma(s+1)}\Gamma\left(\frac{s+1}2\right)\Gamma\left(-\frac{s-1}2\right)\tag1\\ &=\frac1{2\Gamma(s+1)}\frac{\pi}{\sin\left(\pi\frac{s+1}2\right)}\\ &=\frac1{2\Gamma(s+1)}\frac{\pi}{\cos\left(\frac{\pi s}2\right)}\\ &=\frac{\pi}{\Gamma(s+1)}\frac{\sin\left(\frac{\pi s}2\right)}{2\sin\left(\frac{\pi s}2\right)\cos\left(\frac{\pi s}2\right)}\\ &=-\sin\left(\frac{\pi s}2\right)\frac{\pi}{\Gamma(s+1)\sin(\pi(s+1))}\tag2\\ &=-\sin\left(\frac{\pi s}2\right)\Gamma(-s) \end{align*}$$

$$\therefore~\mathfrak{I}=\int_0^{\infty}y^{-s-1}\sin(y)dy~=~-\Gamma(-s)\sin\left(\frac{\pi s}2\right)$$

For the simplification of the final solution we excessively used Euler's Reflection Formula which is a key property of the Gamma Function. Within line $(1)$ we applied the formula for $z=\frac{s+1}2$ and within line $(2)$ for $z=s+1$. The trigonometric reshaping utilized the double-angle formula as well as the periodic property of the sine function.