Find all of the solution of the following differential equation s.t. $\lim_{x\to +\infty}\frac{y(x)}{x^2}=0$ $$y^{(33)}(x)-y^{(30)}(x)=e^{-x}$$
I have immediately dismissed the idea of using the variation of parameters method, since I need 33 linearly independent solutions of the homogeneous; however, putting it into a system is leading nowhere. Any hint?
I would try this: $$y^{(33)}(x)-y^{(30)}(x)=e^{-x}$$ $$y^{(31)}(x)-y^{(28)}(x)=e^{-x}+C_1x+C_2$$
$$y^{(5)}(x)-y^{(2)}(x)=e^{-x}+\sum_{i=0}^{27}C_{i}x^{i}$$ $$y^{(3)}(x)-y(x)=e^{-x}+\sum_{i=0}^{29}C_{i}x^{i}$$ The differential equation is of order 33 so you need 33 constants.
Since the differential equation is linear you can solve this apart: $$y^{(3)}(x)-y(x)=e^{-x}$$ The characteristic polynomial is: $$\implies r^3-1=0 \implies (r-1)(r^2+r+1)=0$$ $$\implies r=1, r=-\dfrac 12 \pm i\dfrac {\sqrt 3}2$$ The solution of the homogeneous equation is: $$y(x)=e^{-x/2} \left (A \cos (\dfrac {\sqrt 3 x}2)+B\sin(\dfrac {\sqrt 3 x}2) \right)+Ce^x$$ For the particular solution try: $$y_p(x)=Fe^{-x}$$ $$-Fe^{-x}-Fe^{-x}=e^{-x} \implies F=-\frac 12$$ So the general solution of the DE is: $$y(x)=e^{-x/2} \left (A \cos (\dfrac {\sqrt 3 x}2)+B\sin(\dfrac {\sqrt 3 x}2) \right)+Ce^x-\frac 12e^{-x}$$
In absence of initial conditions the solution is: $$y(x)=e^{-x/2} \left (A \cos (\dfrac {\sqrt 3 x}2)+B\sin(\dfrac {\sqrt 3 x}2) \right)+Ce^x-\frac 12e^{-x}+\sum_{n=0}^{29}a_{n}x^{n}$$