dsolve(diff(diff(y(x),x),x)+4*x*diff(y(x),x)+(4*x^2+2)*y(x)=0,y(x))