dsolve(x^4*diff(diff(diff(diff(y(x),x),x),x),x)+6*x^3*diff(diff(diff(y(x),x),x),x)+(4*x^4+(-rho^2-sigma^2+7)*x^2)*diff(diff(y(x),x),x)+(16*x^3+(-rho^2-sigma^2+1)*x)*diff(y(x),x)+(rho^2*sigma^2+8*x^2)*y(x)=0,y(x))