Let $a>0,x\geq 0$, the lower regularized incomplete gamma function is defined as : $$P(a,x)=\frac{\gamma(a,x)}{\Gamma(a)} = \int_0^x \frac{e^{-t}t^{a-1}}{\Gamma(a)}dt.$$
I have read in the paper of Gautschi "The incomplete gamma functions since Tricomi" that $P(.,x): a\mapsto \frac{\gamma(a,x)}{\Gamma(a)}$ is decreasing on $\mathbb{R}_+^*$
But how to prove this result ?