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]