I have to solve the system
$y'=-y-2z+2e^{-ax}$,
$z'=3y+4z+e^{-ax}$.
I expressed z from the first equation and plugged it into the second equation. The result is equation $y''-3y'+2y=4(2-a)e^{-ax}$. How to solve this equation?
Any help is welcome. Thanks in advance.
$$y''-3y'+2y=4(2-a)e^{-ax}$$ Solve the homogeneous part first: $$y''-3y'+2y=0$$ Polynomial characteristic is $$r^2-3r+2=(r-1)(r-2)=0 $$ $$ \implies S_r=\{1,2\}$$ $$y_h=c_1e^x+c_2e^{2x}$$ For the particular solution it depends on $a$ if $a \ne -1,-2$ try the particular solution: $$y_p=Ae^{-ax}$$ $a$ if $a = -1$ or $a=-2$ try the particular solution: $$y_p=Axe^{-ax}$$ Plug this in the original inhomogeneous DE and find the constant A then the solution is $$y=y_h+y_p$$
You can add both DE and solve a first order DE instead $$(y+z)'=2(y+z)+3e^{-ax}$$ $$((y+z)e^{-2x})'=3e^{-(a+2)x}$$ This is easy to solve. $$y=-z +3e^{2x}\int e^{-(a+2)x}dx$$ The integral depends on the value of a you have to discuss ($a \ne -2$ or $a=-2$). Plug y in this DE $$z'=3y+4z+e^{-ax}$$ Then you have a first DE in $z$ that's easy to integrate.