Bivariate delta-evolution equations and convolution polynomials: Computing polynomial expansions of solutions

This paper describes an application of Rota and collaborator's ideas, about the foundation on combinatorial theory, to the computing of solutions of some linear functional partial differential equations. We give a dynamical interpretation of the convolution families of polynomials. Concretely, we interpret them as entries in the matrix representation of the exponentials of certain contractive linear operators in the ring of formal power series. This is the starting point to get symbolic solutions for some functional-partial differential equations. We introduce the bivariate convolution product of convolution families to obtain symbolic solutions for natural extensions of functional-evolution equations related to delta-operators. We put some examples to show how these symbolic methods allow us to get closed formulas for solutions of genuine partial differential equations. We create an adequate framework to base theoretically some of the performed constructions and to get some existence and uniqueness results.