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))