Y "+ 9y = xsin3x to solve differential equations

Y "+ 9y = xsin3x to solve differential equations

(x*sin(3*x))/36 - (x^2*cos(3*x))/12 + C1*cos(3*x) + C2*sin(3*x)
Matlab solution, the process of reference books step by step