Use the method of variation of arbitrary constants. The solution to the homogeneous equation is:
[tex]
y_{0}(t) = C_{1} \, \cos{t} + C_{2} \, \sin{t}[/tex]
Then, assume [itex]C_{i} \rightarrow C_{i}(t), i = 1, 2[/itex] where these functions satisfy the following conditions:
[tex]
\left[\begin{array}{cc}<br />
\cos{t} & \sin{t} \\<br />
<br />
-\sin{t} & \cos{t}<br />
\end{array}\right] \cdot \left[\begin{array}{c}<br />
C'_{1}(t) \\<br />
<br />
C'_{2}(t)<br />
\end{array}\right] = \left[\begin{array}{c}<br />
0 \\<br />
g(t)<br />
\end{array}\right][/tex]
The solution for [itex]C'_{i}(t)[/itex] is:
[tex]
\left[\begin{array}{c}<br />
C'_{1}(t) \\<br />
<br />
C'_{2}(t)<br />
\end{array}\right] = \left[\begin{array}{cc}<br />
\cos{t} & -\sin{t} \\<br />
<br />
\sin{t} & \cos{t}<br />
\end{array}\right] \cdot \left[\begin{array}{c}<br />
0 \\<br />
g(t)<br />
\end{array}\right] = \left[\begin{array}{c}<br />
-\sin{t} \, g(t) \\<br />
<br />
\cos{t} \, g(t)<br />
\end{array}\right][/tex]
One integration gives the following:
[tex]
C_{1}(t) = C_{1} - \int_{t_{0}}^{t}{g(t') \, \sin{t'} \, dt'}[/tex]
[tex]
C_{2}(t) = C_{2} + \int_{t_{0}}^{t}{g(t') \, \cos{t'} \, dt'}[/tex]
Substituting this into the expression for the general solution, one gets:
[tex]
y(t) = C_{1} \, \cos{t} + C_{2} \, \sin{t} + \int_{t_{0}}^{t}{g(t') \, \left(\sin{t} \, \cos{t'} - \cos{t} \, \sin{t'}\right) \, dt' \right)}[/tex]
where, the integrating constants [itex]C_{1/2}[/itex], are determined from the initial conditions:
[tex]
y(t_{0}) = C_{1}[/tex]
[tex]
y'(t_{0}) = C_{2}[/tex]
Using the addition theorem for the sine function:
[tex]
\sin{(t - t')} = \sin{t} \, \cos{t'} - \cos{t} \, \sin{t'}[/tex]
we see that the above integral can be written as:
[tex]
\int_{t_{0}}^{t}{g(t') \, \sin{(t - t')} \, dt'}[/tex]
Take [itex]t_{0} = 0[/itex] and compare with the definition for convolution, you will get your desired result.