dsolve(diff(diff(y(x),x),x)+(a*x+b)*diff(y(x),x)+(a1*x^2+b1*x+c1)*y(x)=0,y(x))