First we write the system in differential operator notation:
(D−4)x+D2y(D+1)x+Dy=t2=0. (10)
Then, by eliminating x, we obtain
[(D+1)D2−(D−4)D]y=(D+1)t2−(D−4)0
or (D3+4D)y=t2+2t.
Since the roots of the auxiliary equation m(m2+4)=0 are m1=0,m2=2i, and m3=−2i, the complementary function is
yc=c1+c2cos2t+c3sin2t.
To determine the particular solution yp we use undetermined coefficients by assuming yp=At3+Bt2+Ct. Therefore
yp′=3At2+2Bt+C,yp′′=6At+2B,yp′′′=6Ayp′′′+4yp′=12At2+8Bt+6A+4C=t2+2t.
The last equality implies 12A=1,8B=2,6A+4C=0, and hence A=121,B=41,C=−81. Thus
y=yc+yp=c1+c2cos2t+c3sin2t+121t3+41t2−81t. (11)
Eliminating y from the system (9) leads to
[(D−4)−D(D+1)]x=t2 or (D2+4)x=−t2.
It should be obvious that
xc=c4cos2t+c5sin2t
and that undetermined coefficients can be applied to obtain a particular solution of the form xp=At2+Bt+C. In this case the usual differentiations and algebra yield xp=−41t2+81, and so
x=xc+xp=c4cos2t+c5sin2t−41t2+81. (12)
Now c4 and c5 can be expressed in terms of c2 and c3 by substituting (11) and (12) into either equation of (9). By using the second equation, we find, after combining terms,
(c5−2c4−2c2)sin2t+(2c5+c4+2c3)cos2t=0
so that c5−2c4−2c2=0 and 2c5+c4+2c3=0. Solving for c4 and c5 in terms of c2 and c3 gives c4=−51(4c2+2c3) and c5=51(2c2−4c3). Finally, a solution of (9) is found to be
x(t)=−51(4c2+2c3)cos2t+51(2c2−4c3)sin2t−41t2+81, y(t)=c1+c2cos2t+c3sin2t+121t3+41t2−81t.