I can't figure out how to formulate this in a single expression as the number of options for each additional digit depends on the previous last digit's value, but each previous step represents the number of combinations, not a specific number.
Also, I would like to know the correct way to compute this directly, and not recursively.
Thanks!
So I'd like to implement my comment as an answer.
Let f(number of even digits taken,number of odd digits taken,is last digit even) be the recursive function with initial conditions $f(1,0,1)=4$,$f(0,1,0)=5$,$f(1,0,0)=0$,$f(0,1,1)=0$.
Now, for $m>0,n>0$ we get
$f(m,n,1)=\left(f(m-1,n,1)+f(m-1,n,0)\right)\cdot(4-(m-1))$,
$f(m,n,0)=f(m,n-1,1)\cdot(5-(n-1))$.
And $g(x)=\sum\limits_{m+n=x\\0\le m\le4\\0\le n\le5} (f(m,n,1)+f(m,n,0))$ will give the desired result.
As we have only $6$ values for $m$ and only $5$ values for $n$, the $f$ table will consist of $6\cdot 5\cdot 2=60$ values, not too much to do by hand.
If done right, the table will be similar to this:
f(m,n,1) | f(m,n,0) m=0 m=1 m=2 m=3 m=4| m=0 m=1 m=2 m=3 m=4 n=0 1 4 12 24 24| 1 0 0 0 0 n=1 0 20 120 360 480| 5 20 60 120 120 n=2 0 0 240 1440 2880| 0 80 480 1440 1920 n=3 0 0 0 1440 5760| 0 0 720 4320 8640 n=4 0 0 0 0 2880| 0 0 0 2880 11520 n=5 0 0 0 0 0| 0 0 0 0 2880(source code)