Is there a way to evaluate the derivative of $x$! without using Gamma function?

If subtracting $1$ successively from a number $x$ eventually gives $1$, then $x$ must be an integer, too, and so the definition $$x! := x (x - 1) \cdots (2) (1)$$ is only defined for positive integers $x$. (We can extend this by defining $0!$ to be the value $1$ of the empty product; this is a good choice in the sense that the identity $x! = x (x - 1)!$ then holds for $x = 1$, too.)

In particular, so defined the factorial is a function $\Bbb N \to \Bbb N$, and in particular, its derivative doesn't exist at any point: By definition, if a (real) function $f$ is differentiable at a point $a$ in its domain, $f$ must be defined on some open interval $(a - \epsilon, a + \epsilon)$ (sometimes, and depending on application, we'll allow a half-open interval with $a$ at its endpoint, which is enough to compute a left- or right-hand limit of the difference quotient).

In short, if we want to make sense of $ \frac{d}{dx}(x!$) we must extend $x!$ to some (differentiable) function, like $\Gamma(x + 1)$. (This extension is not arbitrary, by the way, it is the unique extension of $x!$ to $(-1, \infty)$ that satisfies certain natural properties.)


According to your approach regarding $N(x)$ it seems you have a function $x!:\mathbb{R}\rightarrow\mathbb{R}$

\begin{align*} x!:=\prod_{k=0}^{n-1}(x-k) \end{align*}

consisting of $n$ factors $x-k$ in mind. If so, we can write the function using the Pochhammer symbol $$(x)_n=x(x-1)(x-2)\cdots (x-n+1)$$ which can be written as polynomial in $x$ using the Stirling Numbers of the first kind $\begin{bmatrix}n\\k \end{bmatrix}$ \begin{align*} (x)_n=\prod_{k=0}^{n-1}(x-k)=\sum_{k=0}^n\begin{bmatrix}n\\k \end{bmatrix}x^k \end{align*} We can so find a representation for the derivative \begin{align*} \frac{d}{dx}&x(x-1)(x-2)\cdots (x-n+1)=\frac{d}{dx}(x)_n\\ &=\frac{d}{dx}\sum_{k=0}^n\begin{bmatrix}n\\k \end{bmatrix}x^k\\ &=\sum_{k=1}^nk\begin{bmatrix}n\\k \end{bmatrix}x^{k-1}\\ \end{align*}

Alternatively, if you consider generalisations of the factorial function different to the Gamma function you should consider according to @Travis answer which properties a generalisation should have and then analysing the properties regarding derivatives.

A nice page providing interesting alternatives to the standard definition of $x!$ is Is the Gamma function mis-defined by P. Luschny.


Expanding on Lucian's answer (or perhaps making it more clear)

$$f(x):=\frac{d}{dx}\ln(x!)=\frac{\left(\frac{d}{dx}x!\right)}{x!}\tag0$$

$$\frac{d}{dx}\ln(x!)=\frac{d}{dx}\ln(1\cdot2\cdot3\dots x)=\frac{d}{dx}\ln(1)+\ln(2)+\ln(3)+\dots\ln(x)$$

$$f(x)=\frac{d}{dx}\sum_{n=1}^x\ln(n)$$

A discrete representation nonetheless, but we will make do with it.

$$F(x)=\int_1^xf(t)dt+c=\sum_{n=1}^x\ln(n)$$

$$\implies F(x)=\ln(x)+F(x-1)\tag1$$

$$\implies\frac{d}{dx}F(x)=\frac{d}{dx}\ln(x)+F(x-1)$$

$$f(x)=\frac1x+f(x-1)$$

$$f(x)-C=\sum_{n=1}^x\frac1n\tag{backwards $(1)$}$$

$$f(x)=C+\sum_{n=1}^x\frac1n$$

$$\frac{d}{dx}x!=x!\left(C+\sum_{n=1}^x\frac1n\right)$$

where $c$ and $C$ are constants. If we sneak in with the Gamma function, we find the exact value of the constant:

$$\frac{d}{dx}x!=x!\left(-\gamma+\sum_{n=1}^x\frac1n\right)$$

where $\gamma$ is the Euler-Mascheroni constant.


$(0)$ - chain rule/logarithmic differentiation.

$(1)$ - if $F(x):=\sum_{n=1}^xg(n)$, then it comes naturally that $F(x)=g(x)+F(x-1)$ from it's definition. The converse is not true, however, and that is why we get extra constants in $(\text{backwards }(1))$.