Problem:
Formally speaking lets say we have,
$ F_n = a_{k}F_{n-1} + a_{k-1}F_{n-2} + a_{k-2}F_{n-3} + \cdots + a_1F_{n-k} $
Base case being,
$F_0 = b_0\\F_1 = b_1\\ \cdots\\F_{k-1} = b_{k-1}$
We are to find the generating function of this recurrence relation in terms of the coefficients $a_i$ and $b_j$.
My attempts:
From symmetry the solution I get is $f(x)(1 - a_k - a_{k-1}x - a_{k-2}x^2 \cdots - a_{k-(k-1)}x^{k-1}) = b_0 + (b_0+b_1)x + (b_0+b_1+b_2)x^2 + \cdots + (b_0+b_1+b_2+ \cdots b_{k-1})x^{k-1}$
$\implies f(x)= \frac {b_0 + (b_0+b_1)x + (b_0+b_1+b_2)x^2 + \cdots + (b_0+b_1+b_2+ \cdots b_{k-1})x^{k-1}}{1-a_k-a_{k-1}x-a_{k-2}x^2- \cdots -a_{k-(k-1)}x^{k-1}} $
where $f(x)$ is the generating function. But I am not sure I can prove it.
Also I would like to express the $t_{th}$ term ie coefficient of $x^t$ in the above. There should be a way to express in $kxk-matrix$ form $M$ such that $M^t[1][1]$ might have the computed coefficients of $x^t$ .
PS:
I have seen many recurrence relation questions across math.stack exchange like this or this or youtube tutorials like this or this but no one seems to be making any effort to generate a generalised formula for it.
This seems to be the closest where they seem to have done the reverse of my question.
This should be well know as I found there is even a wolframalpha tool doing it