Okie, you know that (and how) factorial n! is defined on the integers. We define Gamma as:
[tex]\Gamma(x)=\int_{0}^{\infty}t^{x-1}e^{-t}dt[/tex]
for all x>0 real numbers. Then you, having completed calc 2, can prove the following:
-the integral in the definition of Gamma does indeed converge when x>0
-for all x>0, [tex]\Gamma(x+1)=x\Gamma(x)[/tex], the "Functional equation" of gamma
-if n is a non-negative integer, then [tex]n!=\Gamma(n+1)[/tex], which justifies calling Gamma an extension of factorial.
Some calculators will use Gamma to make sense of arguments of x! that are not non-negative integers, and are essentially defining [tex]x!=\Gamma(x+1)[/tex] when x is a real number greater than -1. I find it odd that yours does half integral values but not other decimals. I believe you, it just seems like an odd thing to do.
To prove your formula, you can use the identity involving sine CRGreathouse supplied to find [tex]\Gamma(1/2)[/tex], then use the functional equation above to find Gamma at 3/2, 5/2, ... n/2 (using induction) then translate back to factorial to get the equation in your first post. Alternatively, look directly at the integral definition of [tex]\Gamma(1/2)[/tex] and try to relate it back to the Gaussian probability integral.