conway said:
Whenever you "solve" a differential equation, the "solutions" are a set of eigenfunctions. You are always allowed to take linear combinations of these eigenfunctions to create new functions which are also solutions of the original differential equation.
That's why you can take the sinusoidal eigenfunctions of a guitar string and combine them to show what a plucked string with a triangular profile does.
But then the OP's first statement is really trivial, because
any function can be written as a linear combination in a complete basis.
I also wonder if you are correct about the linear combinations being solutions of the original differential equation, at least in the general case ...
Take the following example:
\frac{d^{2}f}{dx^{2}} + cf(x) = 0
two solutions of this are:
f_{1}=cos(ax), which solves it for c=a^{2}, and
f_{2}=sin(bx), which solves it for c=b^{2}
However, if we construct the arbitrary linear combination g=\alpha f_{1} + \beta f_{2}, as far as I can tell, there are
no choices of non-zero values for both \alpha, \beta and c for which g is a solution of the original diff. eq., except for in the degenerate case where a=b, in which case any values of \alpha and \beta will work.
More generally, the whole point of finding the eigenfunctions is that they are a complete set of linearly independent functions. That statement by itself precludes the possibility of finding other solutions that can be written as linear combinations of the eigenfunctions, since then it would not be a complete set, and any such solution would be linearly dependent.