Esta es una respuesta a la solicitud de comentarios a mostrar cómo llegar a la fórmula para las entradas de $L^h$
Aquí está un ejemplo de cómo podemos encontrar las entradas de $L^h$ simbólicamente. Yo lo hago con el $4 \times 4 $ triangular de la matriz de $L$
$$ M(h) = \exp( h \cdot \log (L)) = L^h $$
y obtener por $M$
$$ M(h) = \left[ \begin{array} {rrrr}
1 & . & . & . \\
1h & 1 & . & . \\
\frac 12( 1 h^2+ 1 h) & 1h & 1 & . \\
\frac 16( 1h^3+ 3 h^2+ 2 h) & \frac 12 (1h^2+ 1h) & 1h & 1
\end{array} \right]$$
Aquí un ojo entrenado reconoce los números de Stirling de 1ª clase como los coeficientes de las potencias de $h$ y debido a que la estructura de la matriz tiene esta constante de las diagonales es fácil hacer la fórmula para la transferencia:
$$ M(h)\cdot S_k = S_{k+h} $$
Un paso más se muestra, que la evaluación de los polinomios en las entradas conduce a la binomial números, que es un conocido de la propiedad de los números de Stirling de primera especie (la vandermonde de la matriz de LDU-se descompone en las matrices de Stirling-los números 2º clase y de los coeficientes binomiales y por lo tanto reduce por la multiplicación de la matriz de los números de Stirling 1'st tipo (que es la inversa de la 2ª clase de la matriz) para binomios)
Me he divertido para continuar un poco. Factorización de la smbolic entradas, suponiendo que el hypothese que tenemos siempre los números de Stirling de 1ª clase y la fracción de cofactores de la recíproca fatorials dar
$$ M(h) = \left[ \begin{array} {llll}
1 & . & . & . \\
1(h) & 1 & . & . \\
\frac 12( h(h+1)) & 1(h) & 1 & . \\
\frac 16( h(h+1)(h+2)) & \frac 12 ( h(h+1)) & 1(h) & 1
\end{array} \right]$$
y esto da inmediatamente los coeficientes binomiales
$$ M(h) = \left[ \begin{array} {cccc}
1 & . & . & . \\
\binom{h}{1} & 1 & . & . \\
\binom{h+1}{2} & \binom{h}{1} & 1 & . \\
\binom{h+2}{3} & \binom{h+1}{2} & \binom{h}{1} & 1
\end{array} \right]$$
and a routine could solve the problem given the vector S of dimension n in the following way:
T = 0*S ; \\ initialize an empty array of size of S
b = 1; \\ contains the current binomial
for(j=1,n,
for(k=j,n, T[k]+=S[k+1-j]*b);
b *= (h+j-1)/j;
);
return(T);
So we have $n^2/2$ operaciones por el bucle.