DSolve[Sin[x] - Cos[x]*y[x] - 3*Sin[x]*Derivative[1][y][x] + 3*(1 + Cos[x])*Derivative[2][y][x] + (x + Sin[x])*Derivative[3][y][x] == 0,y[x],x]
dsolve((sin(x)+x)*diff(diff(diff(y(x),x),x),x)+3*(cos(x)+1)*diff(diff(y(x),x),x)-3*sin(x)*diff(y(x),x)-y(x)*cos(x)+sin(x)=0,y(x))