Referring to this functon: $ J_n(x)= \sum_{k=0}^{\infty}\frac{\substack{(-1)^k}} {\substack{k!\Gamma(k+n+1)}}(\frac{\substack{x}} {\substack{2}})^{2k+n} , n\geq0.$
When changing the $n$ subscript value (i.e. determining the order of Bessel function), I know that when graphing order zero, Bessel function peaks about $x=0$ and graphing order one makes it intercept the origin at $(0,0)$. However, I would like to know the reason for this.


Just look at the first few terms. $$J_0(x) = 1 - \frac{1}{4}x^2 + \dots$$ and $$J_1(x) = \frac{1}{2}x - \frac{1}{16} x^3+\dots$$ As with all convergent Maclaurin series, these tell us the behavior near $x=0$.
Perhaps the $\Gamma$ function in there is making it difficult. When $n$ is an integer, those gamma functions can be written as factorials instead.