dsolve(x^2*diff(diff(y(x),x),x)-2*x*diff(y(x),x)-4*y(x)-x*sin(x)-(a*x^2+12*a+4)*cos(x)=0,y(x))